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A FORTRAN IV PROGRAM FOR PREDICTING THE UNCONTROLLED 

DYNAMIC RESPONSE CHARACTERISTICS OF A SPINNING, 

CABLE-CONNECTED, TWO-BODY SPACE STATION 

By William E. Thomas, Jr. 

Manned Spacecraft Center 

SUMMARY 


This paper presents a 12-degree-of-freedom digital computer program for deter- 
mining the nonlinear motion of two rigid bodies connected by massless cables and sub- 
jected to time -dependent external sinusoidal forces and torques. The equations used in 
the program and a sample problem are included. 


INTRODUCTION 


Long-duration space missions may tax the ability of man to withstand long periods 
of weightlessness. Some uncertainty is presently associated with this ability. Because 
of this uncertainty, the study of techniques to provide an artificial gravity environment 
is very desirable. One of the most attractive means of achieving a high ratio of artifi- 
cial gravity to spin velocity (equivalent to that of a large diameter space station) is to 
separate the manned vehicle from its adjacent booster stage by means of flexible cables 
and to spin the two bodies about their composite mass center. The problem of deter- 
mining the dynamic response characteristics of this type of configuration is dealt with 
in this paper by means of a 12-degree-of-freedom digital computer program. 

The disturbances presently allowed for in the program include time-dependent 
external sinusoidal forces and torques on both bodies. These disturbances can be used 
to determine the general dynamic characteristics of a given configuration and to deter- 
mine the effect on these characteristics of altering body and cable parameters and/or 
cable configurations. The equations defining these external disturbances are pro- 
gramed in an external force subroutine distinct from the rest of the program (appen- 
dix A). Equations representing other external disturbances (gravity gradient torques, 
reaction control system jets, associated control equations, et cetera) may be included 
in the program by replacing and/or adding the appropriate punched cards in the force 
subroutine. Perturbations resulting from .internal mass shifts and fixed internal rotat- 
ing masses may be included in the program by adding the appropriate terms to Euler's 
dynamical equations. Instantaneous cable length and forces are calculated in a subrou- 
tine which represents cable free length as a constant. This subroutine may be easily 



modified to study vehicle motion during cable extension or retraction by simply repre- 
senting cable free length as an arbitrary function of time and/or any of the dependent 
variables in the program. Appendix B provides data on the general input and output of 
the program. 


SYMBOLS 


AA, BB 

AF ,AF ,AF 
x, n’ y, n’ z, n 


AG , AG , AG 
x, n’ y, n’ z, n 


a, b, c 


CD 


n 


CG ,CG ,CG 
x, n’ y, n’ z, n 


CG . ,CG . , CG . 

x’,n’ y’,n’ z’,n 


CK 


n 


instantaneous angular velocity of body 2 about 
c - comp’ de S /sec 

defined by equations (18) and (19) 


components of amplitude of the impressed force acting 
on body n along i n ~, j -, and k n -axis, respec- 
tively, lb (n = 1, 2) 


components of amplitude of the impressed torque acting 
on body n about i^-, j^-, and ld-axis, respec- 
tively, in-lb (n = 1, 2) 

components of translational velocity of c. g. ^ relative 
to c *S- com p» directed parallel to the X*Y’ plane and 


perpendicular to the plane of directed radially 

outward from c. g. com p> and directed perpendicular 
to b and in the plane of in/sec 


equivalent viscous damping coefficient for cable n , 
Ib-sec/in. (n = 2, . . . , N + 1) 

components of torque acting through c. g. n about i -, 
j -, and k n ~axis, respectively, because of N cables, 
in-lb (n = 1, 2) 

components of torque acting through c. g. n about i^-, 
j^-, and k^-axis, respectively, because of N cables, 
in-lb (n = 1, 2) 

spring constant of cable n , lb/in. (n = 2, . . . , N + 1) 
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CL 

n 

C, 

f, max 
cable 

n 


c.g. 


comp 


c.g. 


n 


D 

FD 

n 

FK 

n 


unstretched length of cable n , in. (n = 2, . . . , N + 1) 

number of the highest stressed cable 

cable connected to points P 1 and P 0 (n = 2, . . . , 

Ij 11 cty n 

N + 1) 

center of gravity of composite configuration 

center of gravity of body n (n = 1, 2) 

vector sum of a and c, in/sec 

total damping force in cable n , lb (n = 2, . . . , N + 1) 

total spring force in cable n , lb (n = 2, . . , , N + 1) 


F force on the highest stressed cable, lb 

c, max ’ 


F 

x, n 


,F 


Y,n 


,F 


z, n 


components of impressed force acting on body n along 
i -, j and k n ~axis, respectively, lb (n = 1, 2) 


F ,C , F ,C , F ,C 
x, 1 n’ y, 1 n’ z, 1 n 


F x, 2 G n* F y, 2 C n* F z, 2 C n 


G , G , G 
x, n’ y, n’ z, n 

G x, 1, n’ G y, 1, n’ G z, 1, n 
G x, 2, n’ G y, 2, n’ G z, 2, n 


components of force (because of cable n ) acting on body 1 
parallel to ij-, j^-, and k^axis, respectively, lb 
(n = 2, . . . , N + 1) 

components of force (because of cable^) acting on body 2 
parallel to ig-, jg-, and kg-axis, respectively, lb 
(n = 2, . . . , N + 1) 

components of impressed torque acting on body n about 
iJi~, and k^-axis, respectively, in-lb (n = 1, 2) 

components of torque acting through c.g. ^ about ij-, 
j and k^-axis, respectively, because of cable n , 
in-lb (n = 2, . . . , N + 1) 

components of torque acting through c. g. g about ig-, 
j 2 ~, and kg -axis, respectively, because of cable n , 
in-lb (n = 2, . . . , N + 1) 
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I i , ,n ,I j , ,n’ I k’,n 

body n moments of inertia about i’ j’ and 

j n ’ J n ’ 


2 

k^-axis, respectively, lb- sec -in. (n = 1, 2) 

^B’ ^RB’ k RB 

orthogonal pseudorigid body axes fixed at c. g. com p 

i , j , k 
n’ J n’ n 

arbitrary orthogonal axes body fixed at c. g. (n = 1, 2) 


principal axes of inertia for body n (n = 1, 2) 

V^n 

unit vectors directed parallel to and in the positive 
direction of the i^-, j n ~, and k n -axis, respectively 


(n = 1, 2) 

r , p , & 

n ’ J n* n 

unit vectors directed parallel to and in the positive 
direction of the i^-, j^-, and k^-axis, respectively 


(n = 1, 2) 

KK , ZZ 
ir n 

functions defined by equations (13) and (14), (n = 1, 2) 

K) 

matrix of direction cosines for transforming vector 
components from the principal body axes of body n 
to the arbitrary body axes of body n (n = 1, 2) 

M 

n 

2 

mass of body n, lb- sec /in. (n = 1, 2) 

N 

number of cables 

P P 

l,n’ z, n 

points on bodies 1 and 2, respectively, for which rela- 
tive displacement vectors are determined (n = 1, 2, 

. . . , N + 1) (cable attachment points for n = 2, . . . , 
N + 1) 

P P 
*l,n 2,n 

magnitude of relative displacement vector from 

point P t to point P,, , in. (n = 1, 2, . . . , N + 1) 

l, n Ay n 

TF , TF , TF 
x, n’ y, n’ z, n 

components of total force acting on body n along i n ~, 
j -, and k n ~axis, respectively, lb (n = 1, 2) 

TG , TG , TG 
x, n’ y, n’ z, n 

components of total moment acting through c. g. n about 
i^-, j^-, and k^-axis, respectively, in-lb (n = 1, 2) 

t 

time, sec 
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u 


" v" w” 


n 7 n 


n 


X , Y ,Z 
n’ n’ n 


X P, l,n’ Y P, l,n’ Z P, 1, n 


X P, 2, n’ Y P, 2, n’ Z P, 2, n 


components of translational velocity vector of c. g. 
along i j -, and k n ~axis, respectively, in/sec 
(n= 1, 2) 

components of relative displacement vector from 

point Pj n to point Pg parallel to i^ j.-, and 

k^-axis, respectively, in. (n = 1, 2, . . . , N + 1) 


components of relative displacement vector from c. g. ^ 
to P l n along l r , j r> and kj-axls, respectively, 
in. (n = 1, 2, . . . , N + 1) 


components of relative displacement vector from c. g. ^ 
to Pg n along ig-, jg - * and kg-axis, respectively, 
in. (n = 1, 2, . . . , N + 1) 


X,Y,Z 


inertially fixed orthogonal axes 


X’ , Y' , Z’ 


X 


c.g. 




X 1 , Y’ , Z' 
c.g.’ c.g.’ c.g. 


X ,Y , 
n’ n’ 



a 1, n’ a 2, n’ °3, n 

[r] 


orthogonal axes parallel with X-, Y-, and Z-axis, 
respectively, with origin at c.g. 

» * v/v/iiiy 

components of c. g. displacement vector along 

LuIIly 

X-, Y-, and Z-axis, respectively, in. 

components of relative displacement vector from 

c.g- comp to c *S*i alon S X'-, Y'-, and Z'-axis, 
respectively, in. 

components of c.g- n displacement vector along X-, 

Y-, and Z-axis, respectively, in. 

spin-plane pointing error, deg 

functions defined by equations (37), (38), and (39) 

(n = 1, 2, . . . , N + 1) 

general coordinate transformation matrix, defined by 
equation (1) 
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[r] 



y 


orthogonal transformation matrix for transforming 
vector components from the arbitrary body axes of 
body 1 to the arbitrary body axes of body 2 

orthogonal transformation matrix for transforming 
vector components from the inertially fixed axis 
system to the arbitrary body axes of body n (n = 1, 2) 

orthogonal transformation matrix for transforming 
vector components from the inertially fixed axis 
system to the pseudorigid body axis system 

orthogonal transformation matrix for transforming 
vector components from the pseudorigid body axis 
system to the arbitrary body axes of body n (n = 1, 2) 

angle between a and D, deg 


<t> 




\L> ,6 , <p 

^n’ n’ n 


i// > 9 

r s,n’ 


s, n' s, n 


x, n’ y, n’ z, n 


x, n’ y, n’ z, n 


°x, 1L’ *y, 1L’ \ 1L 


general Euler angles 

Euler angles defining angular orientation of arbitrary 
body -fixed axes ig, j g > anc * ^ respect to 

arbitrary body -fixed axes i^, j^, and k^, deg 

Euler angles defining angular orientation of pseudorigid 
body axes i RR , j RB , and k RR with respect to axes 

X’, Y’, and Z\ deg 

Euler angles defining angular orientation of arbitrary 
body-fixed axes i n , j n , and k n with respect to 

inertially fixed axes X, Y, and Z, deg (n = 1, 2) 

Euler angles defining angular orientation of arbitrary 
body -fixed axes i n , j , and k n with respect to 

pseudorigid body axes i RB , j RR , and k RR , deg 

(n = 1, 2) 

body n angular velocity components about i -, j -, 
and k n ~axis, respectively, deg/sec (n = 1, 2) 

body n angular velocity components about i^-, j^-, 
and k^-axis, respectively, deg/sec (n = 1, 2) 

functions defined by equation (33) 
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frequency of the impressed force acting on body n, 

’ n rad/sec (n = 1, 2) 

o> T frequency of the impressed torque acting on body n, 

,n rad/sec (n = 1, 2) 


Subscripts: 

m defines matrix row 

p defines matrix column 

A dot over a symbol indicates differentiation with respect to t. 

ANALYSIS 


Axes Systems and General Vehicle Orientation 

The general body orientation of the spacecraft is shown in figure 1. Each body 
has two body -fixed orthogonal axes systems having origins at the body center of gravity. 
One of the systems must be coincident with the body principal axes. The rotational 
equations of motion for each body are written with respect to this system and thus will 
reduce to Euler’s dynamical equations. Angular orientation of the other axes system 
within the body is completely arbitrary; however, the system is usually located coin- 
cident with geometrically symmetric axes, if such axes exist. This arbitrary system 
is angularly located in an inertial frame by a set of Euler angles defined in figure 2. 

The translational equations of motion for each body are written with respect to this 
arbitrary system. The arbitrary systems in the two bodies are angularly related to 
each other by a set of relative Euler angles which reduce to pitch, yaw, and roll for 
small angles (fig. 3). The two axes systems within a given body are related by a set 
of direction cosines. 


Composite body motion is broken up into pseudorigid body motion and flexible 
body motion. The pseudorigid body is defined as a straight line connecting the two body 
centers of gravity. This pseudorigid body has a body-fixed orthogonal axes system with 


its origin at the composite center of gravity. An auxiliary set of reference axes X', 

Y’ , and Z’ (parallel to the fixed inertial axes and located at the composite center of 
gravity) is used to angularly orient the pseudorigid body in inertial space. Pseudorigid 
body orientation is shown in figure 4. Note that the i^ B axis is coincident with the 

pseudorigid body and directed toward body 2 and that the j BB axis is restricted to the 
X’Y’ plane (a line can have no roll displacement). The angles ^ RB and 0 RB define 


pseudorigid body inertial angular response, and the coordinates 




and 



g. 


define pseudorigid body inertial translational response. 


The flexible cables 
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connecting the two bodies will superimpose structural oscillations on the pseudorigid 
body motion. The structural response of each body is measured by a set of Euler an- 
gles relating the respective arbitrary body axes to the pseudorigid body axes (fig. 5) 
and by monitoring the instantaneous length of the pseudorigid body jP 2 y The 

program will also calculate the instantaneous angular velocity of body 2 about the com- 
posite center of gravity (fig. 6). 

A number of coordinate transformations are required by the program. To sim- 
plify the description of many of these transformations, the following general matrix 
will be required. 



r n 

r 

12 

r 

13 


cos 9 cos x // 

cos 9 sin 0/ 

-sin 9 

[ r ] = 

r 21 

r 22 

r 

23 

= 

-sin 0/ cos 0 + 
sin 0 sin 9 cos 0/ 

COS 0 cos 0/ + 
sin 0 sin 9 sin 0/ 

sin 0 cos 9 


_ r 31 

r 32 

r 

33 


sin 0/ sin 0 + 
cos 0 sin 9 cos 0/ 

-sin 0 cos 0/ + 
cos 0 sin 9 sin 0/ 

cos 0 cos 9 


The principal body axes are related to the arbitrary body axes (for a given body) as 
follows 


/. > 
l 


(V'l 

n 


n 

V n 

>-KI< 

>n > 

k 


k’ 

\ n ) 


K n ) 


where n = 1, 2 and 



is given by 



V 

n 

i • 
n 

5 

i • 
n 

k’ 

n 

V 

n 

r n’ 


^n ’ 

k T 

n 

r ; 

k • 
n 


k • 
n 

k’ 

n 


( 2 ) 


(3) 
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The inertial axes are related to the arbitrary body axes (for a given body) as follows 


/ \ 



/ = \ 

i 



X 

n 




l ” 

>=[ r n] 

< 

Y > 

k 

l n J 



I 



\ / 


( 4 ) 


where n = 1, 2 and 



is given by equation (1) after the following substitutions are 


made: [ T] = r L T =r , \j/ = xj/ , 0=0, and <p = <p where m = 1, 2, 3 
L J L n] mp mp, n’ r n’ n n ’ ’ 

and p = 1, 2, 3. The arbitrary body axes of body 1 are related to the arbitrary body 

axes of body 2 as follows 


(l\ 


V 

\ j 2 

> = [r]< 

h > 

k 9 

\ 2 J 




(5) 


where [ T ] is given_by equation (1) after the following substitutions are made: 

[ r] = [T], r = r , 0=0, = \ //, and 0 = 0 where m = 1 , 2, 3 and p = 1 , 

L J mp mp 

2, 3. The inertial axes are angularly related to the pseudorigid body axes as follows 


/ \ 
X RB 


/=\ 

X 

VRB 


[ 

, k RB, 


\ J 


( 6 ) 


where is given by equation (1) after the following substitutions are made: 

t ] H ■ [ r R ]’ % - r R, mp’ 6 ' 9 RB> * - * RB’ “ d *’ °- ' 0 Where “ ’ *■ 2 ’ 3 
and p = 1, 2, 3. The Euler angles in equation (6) are derived quantities obtained from 
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the following equations. 



The pseudorigid body axes are related to the arbitrary body axes (for a given body) as 
follows 


fO 

n 


(\ \ 


RB 

1 in 

>= [ r R,^ 

j RB / 

k 


k RB 
\ / 

l n J 



(9) 


where n = 1, 2 and is given by equation (1) after the following substitutions 

0 = 9- 0 = and 0 = 


are made: [ r ] = [T ~| 

it; nj 

where m = 1, 2, 3 and 


mp _ ^mp, R, n’ “ ” "s, n’ Y ~ Y s,n’ ' ' s, n 

p = 1, 2, 3. By comparing equations (4) and (6) to equa- 
tion (9), the structural Euler angles can be obtained as follows 


* = tan’ 1 ! J’ 11 > «- Fl2 - 22 

Sj n v -* * __ x . 


11, n* R, 11 + r i2, n r R, 12 + F 13, n F R, 13, 


( 10 ) 


,/KK 

*8, n = tan (zz 


n , 


(ID 
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9 = tan 

s,n 


-l/‘ r il,n r R, 13 ~ r !2, n r R, 32 ~ F 13, n F R, 33^ 




( 12 ) 


2 + ZZ 2 
n 


where KK n and ZZ are given by 

^n = r 21, n r R, 31 + r 22,n r R, 32 + r 23,n F R, 33 

and 


(13) 


ZZ n " r 31, n r R, 31 + F 32, n r R, 32 + r 33, n r R, 33 


(14) 


The instantaneous spin-plane variables (fig. 6) are derived quantities obtained from the 
following equations 


, - 1/ BB \ 

y = tan (sx; 


(15) 


■4 

a = cos -1 ^cos cos y'j 


(W 


and 


A 2 “ 


VbB 2 + AA 2 


P l, 1 P 2, 1 * V(^.g.) 2 + ( T ’c.g.) 2 + ( Z, c.g.) ! 


(17) 


where AA and BB are given by 


^ = (V *c.g.) C0S *RB ' ( X 2 - X c.g.) sin ^RB 


(18) 
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and 


BB = - 




+ 




cos \jj. 


r3 


sin 9 


RB 


(V^c.g .) 008 ^ 


(19) 


Equations of Motion 

The equations used in the program place no restrictions (within the limitations 
imposed by cable interference) on either angular or translational displacement of the 
two rigid bodies relative to an inertial frame and to each other. The bodies may also 
have completely general geometrical and inertial properties. 

The rotational equations of motion used in the program are 


I., & - a W (l. t - L t \ = TG 

i , n x, n y, n z^j^n k r ,ny x, 


n 


I., --a O’ (l, -I., \ = TG 

j',n y, n x,n z,n^k’,n l ',nj y,n 


( 20 ) 

( 21 ) 


and 


I,, ^ 

k’ , n z, 


n 


- 


a 

n y 


jn^iSn " V,n) TG z,n 


( 22 ) 


where n = 1, 2. Integration of these equations yields £2* , £2' , and £? . Equa- 

x ? ii y, n z, ii 

tions (20), (21), and (22) may be modified to include the effects of internal mass shifts 
(for example, crew movements) and fixed internal rotating masses. Body n angular 
velocity components about the arbitrary body axes can be obtained from 


(a ) 

I 

fs? ^ 

x, n 


x,n 

/ £2 

\ y,n 

ii 

r^i 

l2_j 

£2* \ 
y»n / 

£2 


£2' 

l Z ’ n J 


l z ’ n ) 


(23) 
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where n = 1, 2. The translational equations of motion for body n are 


and 


M u” 

+ M / 

r ft 

W M ■ 

- ft 

V” 

n n 

n( 


n n 

z, 

n n 

M v” 

+ M / 

'ft 

u M - 

• ft 

w” 

n n 

n\ 

< z > 

n n 

x, 

n n 

M w” 

+ M i 

(a 

v M ■ 

- ft 

u” 

n n 

n 

\ x 3 

,n n 

y, 

n n 


n 


n 


(24) 

(25) 


(26) 


where n = 1, 2. Integration of these equations yields u”, v” , and w” . The time 
rates of change of the inertial Euler angles for body n are given by 


9 = ft cos 0 - ft sin 0 
n y,n n z,n n 


0 = ft + tan 9 /ft sin 0 + ft cos 0 \ 
n x,n n\ y,n n z,n n) 


(27) 

(28) 


and 


ft sin 0 + ft cos 0 

^ _ y> n n z > n n 

Mi 


cos 9 


(29) 


n 


where n = 1, 2. Integration of these equations results in the Euler angles shown in 
figure 2. The time rates of change of the relative Euler angles are given by 


* = (°y, 2 - Vil) cos 5 ' ( B z, 2 ' n z, l L ) sin * 
t, 1 l) + tan »[( n y, 2 - °y, Xl) 81 " ? + 2 - \ 1l) cos * 


ft 

X, 


(30) 

(31) 
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and 



cos 0 


where 


(n 

x, 1L 


fft ^ 
x, 1 

* fi y, il 

>=m< 

ft , \ 

y» 1 ' 

il , 


ft , 

l z ’ 7 


(33) 


Integration of equations (30), (31), and (32) yields the relative Euler angles. The com- 
ponents (in the arbitrary axis system of body 1) of the time rate of change of P. P 9 


X = Y ft . - Z ft , - uV - Z„ - ft i + , ft 1 

n n z, 1 n y, 1 1 P, 1, n y, 1 P, 1, n z, 1 


+ r il a l,n +r 2l“2,n + r 31 a 3,n 


(34) 


Y = Z ft , - X ft , - vV - X„ , ft i + Zr. i n i 
n n x, 1 n z, 1 1 P, 1, n z, 1 P, 1, n x, 1 


+ r i2“l,n + r 22 a 2,n + r 32 a 3 > n 


(35) 


and 


Z = X ft 


n 


Y ft 


uw H A uw 

n y, 1 n x, 1 


w. 


ft 


+ X, 


ft 


P 4 ww. ^ ~r ■« _ “W -4 

, l,n x, 1 P, l,n y, 1 


+ r i3 a l,n +r 23 a 2,n +r 33 a 3,n 


where n = 1, 2, . . . , N + 1 and where 


(36) 


a 1 - U’’ + Z-j n ft O “ ^T3 9 9 

n 2 Pj 2, n y ? 2 2, n. z ? 2 


(37) 
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(38) 


a 


2, n 


= v" + X. 


£2 


- Z, 


£2 


P,2,n z,2 P, 2, n x, 2 


and 


a 


3, n 


= w" + Y 


£2 


- X, 


£2 


P, 2, n x, 2 P, 2, n y, 2 


(39) 


Integration of equations (34), (35), and (36) yields the components of the relative dis- 
placement vector from point P^ n to point Pg n - 


Force Equations 

The total force acting on each body is made up of forces caused by the elongation 
of the interconnecting cables and forces caused by external sinusoidal forcing functions 
acting on the respective body. The cables are considered to be perfectly elastic tension 
members. The cable spring force is given by 


FK = CK 
n 




Pi P 0 
l,n 2, n 


CL 


») 


(40) 


where n = 2, ..., N + 1. When P. P 0 is less than CL , cable is slack and 
’ l,n 2, n n’ n 

FK n is set equal to zero by the program. Energy absorption per cycle because of 

damping may be approximated by an equivalent viscous damping term provided in the 
equations of cable force. This damping force is given by 


FD = CD 
n n 


XX +YY +ZZ 
n n n n n n 


Pi Po 

1, n 2, n 


(41) 


where n = 2, . . . , N + 1. 

Note that the equilibrium position of a given rotating system may be determined 
by inputting nonzero values for CD n and zeros for all the forcing function amplitudes 

and then allowing the program to run until all structural oscillations damp out. The 
forces FD n and FK n are directed along cable n and are signed in the arbitrary axes 

of body 1. The components of force, caused by cable n acting on body 1 at P^ n , are 
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then given by 


and 


F 1 C = (FK + FD \ 
x, 1 n ^ n n J 

F .C = fFK + FD ^ 
y> 1 n \ n nj 

F -C = fFK + FD N 
z, 1 n \ n n J 



Pi P«i 
l,n 2, n 



P< P« 

1 j 11 6^ II 



Pi Po 
l,n 2, n 


( 42 ) 


(43) 


(44) 


where n = 2, . . . , N + 1. The cable force acting on body 2 will be equal and opposite 

to the cable force acting on body 1. The components of force (in the arbitrary axis sys 

tem of body 2), caused by cable acting on body 2 at P 0 , are given by 

n z. n 


(f 0 c ^ 


(f t c ^ 

x, 2 n 


x, 1 n 

< F y,2 C n 

ii 

i 

S 

F iC > 
y, 1 n 

F „C 


F 

y z,2 n) 


[ z,l n) 


(45) 


where n = 2, . . . , N + 1. The components of torque (in the arbitrary axis system of 
body 1) acting through c.g. caused by cable acting at P 1 , are given by 

X I* X* 11 


G 1 - Y„ - (F ,C ) - Z- 1 (F ,C \ 

x, 1, n P, l,n\ z, 1 nj P, 1, n\ y, 1 nj 

G 1 =Z„ , (F ,C ]-!, , (f ) 

y, l,n P, l,nV x, 1 nj P, l,n\ z, 1 n ) 


(46) 

(47) 
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and 


where n = 2, 


and 


G 1 = X- , (F iC ) - Y d ! (f 1 C \ 

z, l,n P, l,n^ y, 1 n ) P, l,n^ x, 1 n) 


N + 1. Similarly, for body 2 


G n = Y n „ i 

'f 0 c 1 


r F „C ' 

x, 2,n P, 2, n\ 

( Z,2 n) 

P, 2,n( 

^ y,2 n 


f F n C \ 

- „ / 

f F 0 C 

y, 2,n P, 2,n( 

^ x,2 nj 

P, 2, n\ 

^ z ’ 2 n 

G n = - X-r, « j 

(f „c \ 

- Y_, 0 I 

'f 0 C 

z,2,n P,2 ? n ! 

( Y, 2 n) 

P, 2, n' 

^ x, 2 n 


(48) 

(49) 

(50) 

(51) 


where n = 2, . . . , N + 1. The components of torque (in the arbitrary axis system of 
body 1) acting through c. g. caused by all of the cables, are given by 


N+l 

CG , = V* G. - 
x, 1 x, 1, n 

n=2 


(52) 


N+l 

CG y,i = E G y,i,n 

n=2 


(53) 


and 


N+l 

CG . = V 1 g: - 
z, 1 z, 1, n 

n=2 


(54) 
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Similarly, for body 2 


N+l 


CG x, 2 = L G x, 2, n 


(55) 


n=2 


N+l 


CG „ = VG 0 

y,2 y, 2,n 


(56) 


n=2 


and 


N+l 

CG z,2 = E%2,n 

n=2 


(57) 


The components of torque acting through c. g. , caused by all of the cables, may now 
be transformed to the principal axes as follows 


(bo, 1 

x f , n 
< CG y,n 

II 

i2_j 

CG 'I 

x, n 

CG > 

y, n ' 

CG , 

V Z >J 


CG 

l Z ’ n J 


(58) 


where n = 1, 2. The components of external sinusoidal force acting on body n along 
the arbitrary body axes are given by 


F = AF sin^„ t\ 

x, n x, n V F, n / 

F = AF sin(u>_ t) 

y, n y, n \ F,n/ 


(59) 

(60) 


and 


F = AF sin/w„ t) 
z, n z, n \ F,n J 


( 61 ) 
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where n = 1, 2. The components of external sinusoidal torque acting on body n about 
the principal body axes are given by 


G = AG 
x, n x, n 


G = AG 
y»n y,n 


s ^ n ( w T,n t ) 
sln ( w T, n‘) 


(62) 

(63) 


and 


i = AG sin 
z, n z, n 


(“T,n‘) 


(64) 


where n = 1, 2. Any component of external force or torque may be zeroed out by in- 
putting a zero for the amplitude of that component. The components of the total force 
acting on body n along the arbitrary body axes may now be obtained from 


N+l 

TF = F + V' F C 
x, n x, n / y x, n m 

m=2 


(65) 


TF 


y>n 


= F 


y>n 


N+l 

+ Z 


m=2 


F C 
y,n m 


( 66 ) 


and 


TF = F 
z, n z,n 


N+l 

+ 

m=2 



F 

z, 


C 


n m 


(67) 


where n = 1, 2. The components of the total moment acting on body n about the prin- 
cipal body axes are given by 


TG = CG , + G 

x, n x , n x, n 


TG = CG , + G 

y,n y , n y,n 


(68) 

(69) 
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where n = 1, 2. 


CONCLUDING REMARKS 


This paper has presented the six-degree-of -freedom rigid body equations of mo- 
tion for each of two bodies connected by massless cables. A basic computer program 
was presented for determining the dynamic response of the complete configuration 
subject to external sinusoidal forces and torques on both bodies. The program was 
written in subroutine form to facilitate the addition of equations representing other 
perturbations and/or control systems to the basic configuration. 


Manned Spacecraft Center 

National Aeronautics and Space Administration 
Houston, Texas, May 31, 1968 
908-40-01-03-72 
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Figure 5. - Order of rotation for structural Euler angles. 
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APPENDIX A 


COMPUTER PROGRAM LISTING 


The FORTRAN IV source program presented in this appendix presently requires 
a Stromberg- Carlson (S-C) 4020 high-speed microfilm recorder for some of the output. 

The S-C 4020 output section* of the listing runs from SUBROUTINE FILM through SUB- 
ROUTINE RCLOK. The program may be modified for use on a system which does not 
have an S-C 4020 recorder by removing those sections of the program indicated by the 
applicable comment cards. If the program is modified in this manner, SUBROUTINE 
OUTAID should also be modified to print out those output variables which are now output 
in graphical form. 


This section was programed by P. H. 
Section, Manned Spacecraft Center. 


Thornton, Landing and Docking Analysis 
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c r 


C f HE a. CARDS A RE SYSTEM SET-UP CARDS FOR THE MSC 7094 OPERATING SYSTEM 
C AND ARE NOT PART OF THE BASIC PROGRAM. 

S JOB t02T THOMAS 03591 ET2S6 

3. Id JOB ThOMAS GO 


SIdFIC Mm In 

COMMON/Fl_M/OM»ON»TwMf TWN . THM . THN» FOM » FON r F 1M r FIN » SIM > SIN» SEM> SEN, 
lElM»tlN»XNIM#XNlN»TEM» TEN, ELM, ELN, TWEM, T WEN, TIEMPO, THIM» THIN, KUT » 

2 ITmBOP , MOVIE, IGROP, PHI , THETA ,NUM1 , NUM2,NUM3,K WHICH 
COMMON VaR , KNT , KFST » L 
U I MENS I ON VAK(b8U0) 

CALL STAkT 
10 CONTINUE 

C REMOVE Th£ NtXT 4 CARDS IF SC-4020 NOT AVAILABLE 

CAlL FILM ( OM , ON , TWM» TWN, THM , THN » FOM , FON , FIM » FIN» SlM » SIN , SFM , 

1SEn,EIM,lIN,XNIM,XNIN,TEM, TEN , ELM , ELN, TWEM, TWEN, TIEMPO,THIM» THIN. K 
2NT , ITHbOr,MOVIE» IGROP , PHI , THETA , NUMl , NUM2 » NUM3 , K WHICH ) 

CALL CLEmN 
CALL PK 
00 TO 10 
ENU 


C REMOVE THE FOLLOWING -SOkIGIN ACE- CARD IF SC-4020 NOT AVAILABLE 
SORIolN ACE 

4IBFTC SkCO 

SUBROUTINE START 

DIMENSION DYDX(IOO) , VAR l 6000 ) 

Common var 

EQUIVALENCE (VAR(lOl) »DYDX(1) ) 

C ZERO CORE AT INITIAL LOADING 

DO 20 J=1 * d800 
20 VAR(J) = 0.0 

SET DERIVATIVE OF INDEPENDENT VARIABLE WR/T ITSELF EQUAL 
TO ONE 

DYDX(l) = 1.0 
CALL RK 

RETURN 

END 


SIBFTC SRC1 

SUBROUTINE RK 

DIMENSION Y(100) ,DYDX(100) , G ( 100 ) , D ( 1 00 ) *P(6200) , NTEGER ( 1 
1 00) * VAR(6800) 

COMMON VAR, KNT , KFST 

EQUIVALENCE ( VAR ( 1 ) , Y ( 1 ) ) , ( V AR ( 101 ) »DYDX ( 1 ) ) , ( VAR ( 201 > , Q ( 
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1 

2 


2u 


C 

40 

50 

C 


90 

C 


12U 

C 


160 

C 


190 

C 


230 

C 


26U 

C 


300 

C 

330 


1) ) » ( VAR (401) t NTEGEK ( 1 ) ) * (VAR (501) »D(1) ) » (VAR (601) »P(1) ) 
( NTEGER ( 6 ) »N) 

LOAD INPUT DATA INTO COMPUTER 
CALL INPUT 

REWIND 9 
REWIND 11 
REWIND 13 
KNT=0 
KFST=0 

P ( 6964 ) =P ( 5966 ) * 0 . b 
P ( 5965 ) = P ( 5967 ) *0 • b 
CALL UYUXS 
CALL OUTPUT 

IF ( Y ( 1 ) ~P (2) ) 40 » 330 t 330 

calculate the delta yij) at yid = o 

DO 50 J - 1 f N 
D ( J ) = UYDX ( J ) *P ( 1 ) 

CALCULATE THE Y(J) AT T = 0 
DO 90 J = 1*N 
R = . 5* ( U ( J ) - Q( J) ) 

Y ( J) = Y(u) + R 
Q ( J) = Q(b) + 3 • 0*R - . 5*D ( J) 

CALCULATE DELTA Y(J) AT Y ( 1 ) = HALF STEP 
CALL UYUXS 
DO 120 J — 1 » N 
D ( J) = UYuX(J)*P(l) 

CALCULATE THE Y(J) AT Y ( 1 ) = HALF STEP 
DO 160 o = 1 » N 
R = .292893219* (D( J) - o(J)) 

Y(J)=Y(J)+R 

QIJ) = O(J) + 3 • 0*R - • 2928932 19*D ( J ) 

CALCULATE THE DELTA Y(J) AT Y(l) = HALF STEP (AGAIN) 

CALL UYUXS 

DO 190 J = 1 r N 

D ( J) = UYDX ( J ) *P ( 1 ) 

CALCULATE ThE Y(J) AT Y ( 1 ) = HALF STEP (AGAIN) 

DO 230 J = IfN 
R = 1*70710 678* (D(J) - Q(J)) 

Y ( J ) = Y(J) + R 

Q ( J ) = y(J) + 3 . 0*R - 1 • 70710678*0 (J) 

CALCULATE THE DELTA Y(J) AT Y(l) = FULL STEP 
CALL UYUXS 
DO 260 O = 1 r N 
D( J) = UYDX ( J ) *P ( 1 ) 

CALCULATE THE Y(J) AT Y(l) = FULL STFP 
DO 500 J = 1 r N 

R = . 166666666E+00* (D ( J) - 2.0*Q(J)) 

Y ( J ) = Y ( «J ) + R 
Q ( J ) = Q(J) + 3. 0*R - • 5*U ( J) 

PROCEED TO THE NEXT INTEGRATION STEP 

GO TO 20 

RETURN 

END 



SIBFTC SRC2 


C 


C 


L 


C 


1 

1 

2 


30 

£>00 


SuBKOUTINt INPUT 

DIMENSION Y(100) *Q(100) * FIRSTY (100) » P ( 6200 )> NTEGER ( 100 ) »V 
AR(b800) 

COMMON VAR 

EQUIVALENCE ( V AR ( 1 ) » Y ( 1 ) ) * (VAR (201) #0(1) ) * ( VAR (301 ) » 
FIRSTY ( 1 ) ) f (VAR (401) * NTEGER ( 1) ) r ( VAR ( 601 ) » P ( 1 ) ) 

* (NTEGER(b) »N) * (NTEGER(2) *NP) » (NTEGER(44) »NPAGE ) 

SET PAGE NO. OF FIRST PAGE 
NPAGE = 1 
CALL PAGEHD 

READ CONTROL INTEGERS INTO PROBLEM 
READ (5*30) (NTEGER( J) * J=l*8) 

FORMAT (815) 

WRITE (6»500) (NTEGER( J) » J=l»8) 

FORMAT ( lhOtilb) 

NTEGER ( 2i ) =NTEGER ( 4 ) + 1 
NTEGER ( 24 )=NT EGER (6) - 1 
NTEGER ( 2b ) =N7 EGER ( 7 ) 


NTEGER ( 2b ) =NTEGER ( 6 ) 

NTEGER ( 6 ) INTEGER (5) 

check for individual floating point data entry 

IF(NP) 380 * 380 *110 
110 DO 140 J = 1 »NP 

READ (5*130)1* (P( I ) ) 

130 FORMAT (I5»E15.0) 

WRITE (6* 150 ) I * P ( I ) 

150 FORMAT( Ib»E20.6) 

140 CONTINUE 

P ( 1201 ) = 1 • 0E36 
380 CALL INAID 

ZERO THE Q AND SET IN IC 
DO 420 J = 1 * N 
Q ( J ) = 0.0 
Y ( J ) = FIRSTY ( J) 

A20 CONI INUE 

RETURN 
END 


SltSFTC SRC 3 


SUBROUTINE PAGEHD 

DIMENSION NTEGER(IOQ) *VAR(6800) 

COMMON VAR 

EQUIVALENCE ( VAR ( 401 ) r NTEGER ( 1 ) ) » (NTEGER(44) » NPAGE) * 
1 ( NTEGER (l)r IDENT ) 

WRITE (6*20) IDENT * NPAGE 

20 FORMAT ( 17H1 FORFUN OPTION 15* 56H 

1 PAGE NO 15 ) 
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•p u tv 


RETURN 

END 


i lot- i C SkC4 

SUBROUTINE IMA ID 

DIMENSION FIRSTY(IOO) fP(620P) ► NTEGER ( 1 00 ) » VAR ( 6800 ) 
COMMON VAK 

Equivalence (vak( 301) r firsty ( i ) ) » (var(4oi) »nteger(i> ) > (v a 

1 R (601 ) >P(1) ) > (NTEGER (41 ) *LPRINT) » (NTEGER(42) »NLINE) * (NT£G 

2 EK ( 43 ) * NSn IP ) * ( NTEGER ( 24 ) »NTSKIP) 

EQUIVALENCE ( P ( 10 ) » DL1 1 1 ) > ( P ( 1 1 ) » DL1 1 2 ) r ( P < 12 ) * DL1 13 ) » ( P ( 
1 13) »PL121) DPI 14) »DL1 22) » (P(lb) »DL123) » (P(16) »DL131 ) » (P(l 

7) » UL132 ) * ( P ( IB) t UL133 ) » (P(26) rDL.211) r (P(27) »DL212) r(P(28 
) • DL213 ) » (P(29 ) f DL221 ) f l P ( 30 ) »DL222) » ( P ( 31 ) » PL223) » (P(32) 
* DL231 ) t (P(33) t DL232) t (P(34) rUI-233) 

P ( 107 ) = uLi 1 1*P ( 970 ) +0L121*P ( 971 ) +DL131*P ( 972 ) 

P ( 1 UP ) =UL1 1 2*P ( 970 ) +DL122*P ( 97 1 ) +0L132*P ( 972 ) 

P( 109)= ULl 13*P(970)+DL123*P( 97 1 ) +DL133*P ( 97? ) 

P ( IIP ) — uL2l 1*P ( 98U ) +DL221*P ( 9bl ) +DL231*P ( 982 ) 
Plli7)=uL2l2*P(98U)+DL222*P(9til ) +UL232+P ( 982 ) 

P ( IIP ) -uL2l3*P ( 980 ) +DL223+P ( 9bl ) +UL233*P ( 982 > 

C SbT IN INITIAL CONDITIONS 

FiRbTY(2) = P(107)/57. 295779b 
FiRSTY ( 5 ) = P(106)/57.29b779b 

F IRSTY ( 4 ) = P(109)/57. 295779b 

FiRSTY(b) = P(110) 

FIRSTY (fa) = P(lll) 

FIRST r (7) = P ( 112) 

FIRSTY (b) = P(113)/57. 295779b 

F iRSTY ( 9 ) = P(ll4)/57. 295779b 

FIRSTY (10) = P(115)/57. 295779b 
FIRSTY (11) = P(llb)/57. 2957795 
FIRSTY ( 12) = P(117)/57. 295779b 
FiRSTY ( 13 ) = P(llb)/57. 295779b 
F IRSTY ( 14 ) - P ( 1 19 ) 

F IRSTY (15) = P ( 120 ) 

FIRSTY ( 16 ) = P ( 121 ) 

FIRST Y ( 17 ) = P(122)/57. 295779b 
F IRSTY ( 18 ) = P(123)/57. 2957795 
FIRSTY ( 19 ) = P(124)/57. 295779b 
P ( 55 ) =SlN (FIRSTY ( 8 ) ) 

P(36)=C0S(FIRSTY(8) ) 

P(38)=SlN(FlKSTY(17) ) 

P ( 59 ) =COS ( FIRSTY ( 17 ) ) 

P(44)=SIN(FIRSTY(9) ) 

P(45)=C0S(FIKSTY(9) ) 

P(46)— Si N (FIRSTY (lb) ) 

P ( 4-7 ) =C05 ( FIRSTY ( lb ) ) 

P ( 5976) =SIN( FIRSTY (19) ) 

P ( 5977 ) -5 IN (FIRSTY(lU)) 

P ( 5978 )=C0S(F IRSTY ( 19 ) ) 
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W IN) H 


c 


450 

C 


P ( 5979) =COS(F IRSTY 1 1 0 ) ) 

A2=P(39)*P(5978) 

B2=P(39)*P(5976) 

C2=-P(38) 

U2=P ( 46 ) *P ( 38 ) *P ( 5978 ) -P ( 5976 ) *P ( 47) 

E2=P ( 47 ) *P ( 5978 ) +P ( 4b ) *P ( 38 ) *P ( 5976 ) 

F2=P(46)*P(39) 

t>2=P ( 5976 ) *P ( 46 ) +P ( 47 ) *P ( 38 ) *P ( 5978 ) 

H2=P ( 47 ) *P ( 38 ) *P ( 5976 ) -P ( 46 ) *P ( 5978 ) 

AI2=P(47)*P(39) 

Al=P(36)*P(5979) 

bl=P(36)*P(5977) 

Ci=-p ( 35) 

U 1=P ( 44 ) *P ( 35 ) *P ( 5979 ) -P ( 5977 ) *P ( 45 ) 

E 1=P ( 45 ) *P ( 5979 ) +P ( 44 ) *P ( 35 ) *P ( 5977 ) 

Fl=p(44)*P(36) 

61=P ( 5977 ) *P l 44 ) +P ( 45) *P ( 35) *P ( 5979) 

H1=P ( 45 ) *P ( 35 ) *P ( 5977 ) -P ( 44 ) *P ( 5979 ) 

AI1=P(45)*P(36) 

XXPR=A2*0l+B2*Hl+C2*All 
AKKPR= D2*Gl+F2*Hl+F2*AIl 
ZZPR=G2*G1+H2*H1+AI2*AI1 

FIRSTY(20)=ATAN2( (-XXPR) » ( SORT l AKKPR**2+ZZPR**2 ) ) ) 

F IRSTY ( 21 ) =ATAN2 ( AKKPR » ZZPR ) 

F IRSTY ( 22 )=ATAN2( ( A2*D1+B2*E1 ) » ( A2*Al+b2*Bl+C2*Cl ) ) 


FlKSTYt23)=P(5992) 
F IRSTY ( 24 ) —P ( 5993 ) 
F IRSTY ( 25 ) -P l 5994 ) 
FIRSTY(2b)=P(5995) 
F IRSTY ( 27 ) — P ( 5996 ) 
F IRSTY ( 26 )=P( 5997) 


SET IN CABLE INITIAL CONDITIONS 
NCA6LE = NTEGER (21 ) 

NRESRV = NTEGER13) 

NDO = 3+NCABLE 
DO 450 J = If NDO 
NPUT = J + 22 + NRESRV 
F1RSTY (NPUT) = P(J+139) 

SET IN CONTROL NUMBERS FOR PRINTING 

NLINE = 55 

NSKIP = NTSKIP 

LPRINT = 0 

RETURN 

END 


lldFTC SRC5 

SUBROUTINE DYDXS 

DIMENSION Y(100) »DYDX( 100) »P (6200) » NTEGER (100) » VAR (6800) » 
X1PK20) »Y1P1(20) »Z1P1(20) »X2P2(20) »Y2P2(20) »Z2P2(20) r 
Al(20) »A2(20) rA3(20) fFXH(20) »FY1I(20) »FZ1I(20) »FX2I(20) » 
FY2K20) »FZ2I (20) >6X1(20) »GY1 (20) »GZ1 (20) *GX2(20) »GY2(20> 
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cs^io'oi-pojrot-* \CCD^jO''cr-£'OJi\ 5 i-* vc^'JO'cr-POjrvt— ic r >i ^ tr p w wh -.0 a> ->i O' cr -p c*i ro i-* tr-p 


*bZ2(20) » XBH ( 20 ) »YRR(20) »ZBk(20) >XRRD(20) »Y6RD(20) »ZBRD( 
2U ) 

COMMON vak 

EGUi VALENCF (\/Ak( 1) » Y ( 1 ) ) d VAK ( 10 1 ) » DYDX ( 1 ) ) * (VAR (401) » 
NlEbFK(l) ) » (VAR (601) #P(l))dY(2) rOMXlP) r (DYDX (2) »0MX1PP) » 
(Y(3) fOMYlP) f (DYDX (3) »0|v,YlPu) » ( Y (4) »0MZ1P ) dY(b) »U1PP) r 
( uYuX ( 5 ) » U 1 PPD ) , (Y(6) »V1PP) > (DYDX (6) »V1PPD) t ( Y ( 7 ) #W1PP) » 
(DYDX (7) » V» It J PD ) d Y ( 8) »ThT1) » (DYDX(R) t THT1D) • ( Y ( 9) » PHI 1 ) » 
(uYuX (9) fPHUU) » (Y(10) »PSI1) » (DYDX (in) ,PSI1D) r ( Y ( 1 1 ) »0MX2 
P) f (DYDX (11) t OMX2PD ) dY(12) »0MY2P) » ( nYDX ( 12 ) » OMY2PD) » 

( Y(13) » 0MZ2P ) » (DYDX (13) »0MZ2PD) dY(14) »U2PP) » (DYDX (14) » 
U2PP0) dY(lb) f V2PP ) » (DYuX(lb) »V2PPD) dY(16) » W2PP ) * (DYDX 
(16) t W2PPU ) f (Y (17) • THT2 ) * (DYDX (17) » THT2D ) dY(18) »PHI2) 
EoUl VALENCE (DYDXdB) » PH 1 2D) dY( 19) »PSI2) * (DYDX (19) »PSI2D 
) r (Y(20) » THTBR) » (DYDX (20) »THTbRD ) * (Y(21) »PHIRR) , (DYDX (21) 
rPHXBKD ) * (Y(22) »PSlBR) * (DYDX (22) *PSIRRD) dP(3) »CIXX1) » 
(P(4)»CIYY1)»(P(5) *CIZZ1) dP(6) # CM1 ) dP( 10) »DL111)» (POl) 
• uLl 12 ) » (P( 12) * DL1 13 ) dp(13) »uL121) » ( P ( 14 ) t DL 122 ) r (P(15) » 
DH23) » (P(lb> * DlI 31 ) dP(17) ,DL132) dp(18) »0L133) » (P(19) , 
CiXx?) » (P(2u) fClYY2) t (P(21) *C1ZZ2) » (P(22) »CM2) r (P(26 ) t 
Dl21 1 ) » ( P ( 27 ) t DL212 ) ► (P(28) #DL213) * (P(29) »DL221) » (P(30) , 
DL222 ) dP(3l) t DL223 ) dP(32) >PL231) dP(33) »DL232) dP(34) » 
Du233) dP(3b) rSTHTl) dP( 36) t CTHT 1 ) » (P(37) » TTHT 1 ) 

EutJi VALENCE (P(38) * STHT2 ) dP(39) r CTHT2 ) t (P(40) * 1THT2) » 

( P ( 41 ) » STHTiiR ) » (P(42) » CTHTRR ) dP<43) dTHTBR) * (P(44) rSPHIl 
) r (P(45) t CPhI 1 ) , (P(46) * SPHI2 ) dP(47) »CPHI2) » (PUP) »SPHIBR 
) » (P(49) rCPHlbP) d P ( 5U ) > GX1P > d P ( bl ) »GYlP) dP(52) »GZ1P) * ( 
P ( 53 ) f GX2P ) dP(b4) » GY2P) > (P(5b) »GZ2P) * (P(56) »OMxl) dP(57) 
fUMYl) » (P(5H) >OMZl) » (P(b9) »0MX2) » (P(60) »0MY2) r (P(fel) »OMZ2 
) » (P(62)»U1) » (P(63) »V1) r (P(b4) »W1) » (P(65) »U2) » (P(66) »V2) » 
( P ( 67 ) » W2 ) > ( P ( 68 ) »GAMttll) » (P(fo9) »GAMR12) » (P(70) t GAMB13) » 

( P ( 7 1 ) t GAMB2J ) » ( P (72 ) * GAMb22 ) t (P(73) * GAMB23 ) » (P(74) » GAMB 
31) » ( P ( 75 ) • GAMB32 ) » (P(7b) t GAMB33) » (P(77) »SPSIBR) 

EOUi VALENCE (P(78 ) »CPSIbR) » ( P ( 79 ) » GX1T) » ( P ( 80 ) »GY1T) » (P(8 
1) f GZ1T ) dp (82) r GX2T ) f (P(83) » GY2T ) » (P(84) »GZ2T) » (P(85) » 
FaIIT) » ( P ( 86 ) » F Y1 IT ) * ( P l 87 ) » FZ1 1 T ) » ( P ( 88) > FX2IT) dP(89) t 
FY21T) • ( P ( 90 ) » FZ2IT ) dP(91) rGXIPT) » (P(92) rGYlPT) » (P(93) » 
GZlPT ) » (P(94) t GX2PT ) » (P(95) »GY2PT) » (P(96) »GZ2PT) » (P(97) » 
S1NP.H2) » (P(98) » COSPH2 ) » (P(99) »COSTH2) • (P(100) rSlNPHl) » (P( 
lUl) fCOSPHl) r l P ( 102 ) »C0STH1) » (P(221) rXIPKl) ) » (P(241) , 
YlPl(l) ) , ( P ( 261 ) »Z1P1 (1) ) » ( P ( 281 ) t X2P2 ( 1 ) ) f (P(301 ) »Y2P2(1 
) ) »(P (321 ) * Z2P2 ( 1 ) ) t ( P ( 341 ) »A1(1) ) d P ( 361 ) » A2 ( 1 ) ) t (P(3P1) 
»A3(1) ) » ( P ( 461 ) fFXll(l) ) r (P(4bl) rFYH(l) ) 

EOUI VALENCE (P(501) »FZ1I (1) ) dP(561) »FZ2I(1) ) » (P(701) »GX1 
( 1) ) » (P (721 ) » GY1 ( 1 ) ) » (P(741 ) r GZ1 ( 1 ) ) » (P (761 ) * 6X2 ( 1 ) ) * 

( P ( 7 8 1 ) »GY2( 1) ) dP(801) t GZ2 ( 1 ) ) , (P(621) rXBR(l) ) » (P(841) , 
YbR(l) ) » (P(86l) t ZBR ( 1 ) ) » ( P (881 ) t XBRD ( 1 ) ) dP(901) * YBRD ( 1 ) ) 
» ( P ( 921 ) t ZBRD ( 1 ) ) » ( P ( 128) » AGX1PT) * ( P ( 129) > A6Y1PT ) » (P(130) 
r AGZ1PT ) » ( P ( 131 ) t AGX2PT ) dp(132) t AGY2PT) » (P ( 133) » AGZ2PT) » 
(P(134) f AFX1IT) t (P(135) » AFY1IT) » ( P ( 136) » AFZ1 IT) dP(137) » 
APX2IT) t ( P ( 1 38 ) f AFY2IT ) » (P(139) t AFZ2IT ) r (DYDX (4) » 

OMZ1PU) » ( P ( 521 ) » FX2 1(1) ) » (P(541) *FY2I ( 1) ) 

EGUI VALENCE (P ( 5976 ) »SPSI2E) » (P( 5977) »SPSI1E) » (P(5978) »CPSI2E) r 



1 (P(5979) fCPSHE) 

St-T UP THE RELATIVE CONSTANT 
NCABLE = NTEGER (21 ) 

NRESRV = NTEGER (3) 

SET IN X» Y ► Z BAR VALUES WHICH RESULT FROM THE 

INTEGRATION 

DO 12 J = lrNCAttlaC^. 

JUMP1 = 3*J + 20 f NRESRV 
XBRlJ) = Y(JUMPl) 

YbRCJ) = Y ( JUMP1 + 1) 

ZBR(J) = YtJUMPl + 2) 

CALCULATE TRIGNOMETRIC FUNCTIONS 
STHT1 = SIN(THTl) 

CTHT1 = COS(THTl) 

TTHT1 = STHT1/CTHT1 
STHT2 = SIN ( THT2 ) 

CTHT2 5 C0S1THT2) 

TTHT2 = STHT2/CTHT2 
STHTBR = SIN(THTBR) 

CTHTBR = COS(THTBR) 


CPSI1E=C0SIPSI1) 

SPSI1E=SINIPSI1) 

CPSI2E=CuS(PSI2> 

SPSI2E=SIN(PSI2) 



TTH1BR = STHTbR/CTHTRR 
SPHI1 '= SIN ( PHI 1 ) 

CPHI1 = CoS(Phll) 

SPHI2 = SIN(PHl2) 

CPHI2 = C0SIPHI2) 

SPHIBR = SlN(PHIBR) 

CPHIBR = COS(PHIBR) 

SPSIBR = SIN(PSIBR) 

CPSIBR = COS(PSIBR) 

CALCULATE GAMMA BAR VALUES FROM TRIG FUNCTIONS 
GAMBll = CTHTBR+CPS1BR 

GAMbl2 = CTHTBR* SPSIBR 

GAMB13 = -STHTBR 

GAMB21 = -CPHIBR*SPSIBR + SPHIBR*STHTBR*CPSIBR 

GAMB22 = CPHIBR*CPSIBR + SPHIBR*STHTBR*SPSIBR 

GAMB23 = SPHIBR*CTHTBR 

GAMB31 = SPHIBR*SPSIBR + CPHIBR*STHTBR*CPSIBR 

GAMB32 = -SPHIBR*CPSIBR + CPHIBR*STHTBR*SPSIBR 

GAMB33 = CPHIBR*CTHTBR 

TRANSFORM PRINCIPAL AXIS ANGULAR VELOCITIES INTO SYMMETRY 
AXIS COMPONENTS 

OMX1 = DL111+OMX1P + DL1 12*0MY IP + DL113*0MZ1P 

OMY1 = UL121*0MX1P + DL122*0MY1P + DL123*0MZ1P - 

OMZI = DL131*0MX1P + DL132*0MY1P t DL133*0MZ1P 

0MX2 = DL211*0MX2P + DL212*0MY2P + DL213*0MZ2P 

0MY2 = DL221+0MX2P ♦ DL222*0MY2P + DL223*0MZ2P 

0MZ2 = DL231*OMX2P + DL232*0MY2P ♦ DL233*0MZ2P 

U1 = U1PH 
VI = V1PP 

wl = W1PP ;v -V , ■ 
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U? = U2PP 
^2 = V2PP 
w2 = W2PP 

pMXlPP=Of»,Xl*GAMbll+OHYl*GAMR12-OMZl*STHTBR 
OMY1PP=0 nX1*GAMR21+OmY1*GAMB22+0mZ1*GAMR23 
0M21PP=0mX l*GAMHb 1 +0 mY1*GAMB32+0MZ1*GAMB33 
C CALCULATE THETA# PHI# PSI BAR DERIVATIVES 

IHIBRn-CPHiBK* ( Oi 1 Y 2—OMY 1 PP ) “SPHltoR* ( 0M^.2^0MZ1PP ) 

PH1BRU— 0i''iX2-oMXlPP+TTHTRR*SPHIBP* ( OMY 2— 0MY1PP)+TTHTRR*CPHIBR* 
K0MZ2-OM41PP) 

PSiBHD-SKHiBH*(0i'':Y2--of-'YlPP)/CTHTbR+CPHlBK*(0MZ2-0MZlPP)/CTHTbR 
O CALCULATE TUT# PHI# PSI DERIVATIVES 

ThTlD = CPH I l*OMY 1 - SPhIl*0M2t 
PHliD = OMX1 + TTHT1*(SPHI1*0MY1 + CPHI 1+OMZI ) 

PS I ID = (GPHI1*0MY1 + CPHI l*OhZ 1 ) /CTHT 1 
THT2D = CPHI2*0MY2 - SPHI2*0M2? 

Ptil2H = OMX2 + TTHT2*(SPHI2*D|V|Y? + CPHl2*0MZ?) 

Pbl 2 D = (SP,Hl2*0MY2 + CPHI2*0|v,Z2 ) /CTHT2 
C NOW CALCULATE THE ALPHA VALUFS 

DO t)2U o — lfNCABLE 

Al(^) = 112 + Z2P2 ( J ) *OMY 2 - Y2P2 ( J ) *OMZ2 

A2(j) = V2 + X2P2(J)*0M22 - Z2P2(U)*OMX2 

A3(u) = W2 + Y2P2(J)*0MX2 - X2P2(U)*0MY2 

TnEN CALCULATE X# Y# Z bAR UFKIVATIVES FOR EACH 
ATTACHMENT point 

XbRU(u) = YbR(J)*OMZl - ZbR ( J ) *OM Y 1 - U1 

1 - Z1P1(U)*OmY 1 + Y 1P1 ( J ) *OMZl 

2 +oAmR11*Ai( J) +GAMB21*A2( J) +GAMB31*A3 (0) 

YuRo(J) = Zt‘R( J)*OMXl - XbR ( J) *OMZl - V 1 

1 - XiPl(u)*0'lZl + Z1P1 ( J) *OMXl 

2 + GAMbl2*Al(J) + GAMB22*A2(U) + GAMb32*A3(J) 

o2U ZbRU(J) = XbR ( J) *OMY 1 - YBR(J)*OMXl - W1 

1 - Y1P1(j)»OMV 1 + Xl PI ( J ) *OMY 1 

2 + GAM»13*Al ( J) + GAMB23*A2(U) + GAMB33*A3(J) 

I RaNSR tP TO THE CAbLE FORCE SUBROUTINF 
CAl_L C Abh OR 

TRANSFORM SYMMETRY AXIS FORCFS IN BODY 1 INTO SYMMETRY 
AXIS FORCES IN BODY 2 
DO 790 o = 1#NCABLE 

FX2KJ) = - GAMBl 1*FX1 1 ( J ) - GAMB12*FY1 I ( 0) 

1 - GAMB13+FZ1I ( J) 

FY2KJ) = - GAMb21*FXlI ( J) - gAMR22*FY 1 1 ( J ) 

1 - GAMB23*FZ1I ( J) 

F22KJ) = - GAMB31+FX1I (J) - GAMR32*EY 1 I ( J ) 

1 - GAMB33+FZ1I ( J) 

c calculate symmetry axis moments on both bodifs 

GX1(J) = Y1P1 (J)*FZ1I(J) - Z1P1 ( J) *F Y 1 1 ( J ) 

GYi(U)= 2 IPI ( J ) *FXl I ( J ) - X1P1 ( J) *F 71 1 ( J ) 

GZ1(J)= A IP 1 ( J ) *FY 1 1 ( J ) - Y1P1 l J ) *F X 1 1 ( J ) 

0X2(0)= Y2P2 ( J) *FZ2 I ( O) - Z2P2 tJ)*FY2l(J) 

GY 2 ( J ) = <£. 2P2 (J)*FX2I(J) - X2P2 (J>*f-72I(J) 

790 uZ2(0) = X2P2 ( J) *F Y2I (<J) - Y2P? (J)*FX2I(J) 

C NOW SUM THE SYMMETRY AXIS COMPONENTS OF MOMENT 

Gxir = 0.0 
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GY 11 = O.U 

GZ11 = 0.0 

GX2T = 0.0 

GY2T = 0.0 

GZ2T = 0.0 

DO V20 J = 1 » NCABLE 

GX1T = GXiT + GXKJ) 

GY1T = bYlT + GYl(J) 

GZ1T = GZiT + GZ1(J) 

GX2T = GX2T + GX2(J) 

GY2T = GY2T + GY2(J) 

V20 GZ2T = GZ2T + GZ2(J> 

C NEXT » SUM THE SYMMETRY AXIS COMPONENTS OF FORCE 

FX1IT = 0.0 
FYliT = 0.0 
FZliT = 0.0 
FX2.LT = 0.0 
FY2IT = 0.0 
FZ2LT = 0.0 
Du lObO J = 1 f NCABLE 
FxliT = FX1IT + FX1KJ) 

FYliT = FY1IT + FYII(J) 

FZHT = FzlIT + FZII(J) 

FX2LT = FX2IT + FX2KJ) 

FY2LT = FY2IT + FY2I(J) 
lObu FZ2iT = FZ21T + FZ2KJ) 

TRANSFORM SYMMETRY AXIS COMPONENTS OF TOTAL MOMENT INTO 
PRINCIPAL AXIS COMPONENTS 

GX1PT = DL111*GX1T + UL121*GY1T + DL131*GZ1T 

GY1PT = DL112*GX1T + DLl22*bYlT + DL132*GZ1T 

GZ1PT = DL113+GX1T + DL123*GY1T + DL133*GZ1T 

GX2PT = DL21 1*GX2T + DL221*GY2T + DL231+GZ2T 

GY2PT = DL212*GX2T + DL222*GY2T + DL232*GZ2T 

GZ2PT = DL213*GX2T + DL223*GY2T + DL233*GZ2T 

C CALL FORCING FUNCTION SUBROUTINE 

CALL FORPUN 

C CALCULATE THE ANGULAR VELOCITY DERIVATIVES 

OMX1PO = (GX1PT + 0MY1P*0MZ1P*(CIYY1-CIZZ1) )/ClXXl 
1 + AGX1PT/CIXX1 

OMY1PD = (GY1PT + 0MX1P*0MZ1P*(CIZZ1-CIXX1) )/ClYYl 
1 + AGY1PT/CIYY1 

0MZ1PD = (GZ1PT + 0MX1P*0MY1P*(CIXX1-CIYY1) )/ClZZl 
1 + AOZ1PT/CIZZ1 

0MX2PU = (GX2PT + 0MY2P*0MZ2P* ( Cl YY2-CIZZ2 ) ) /CIXX2 
1 + AOX2PT/CIXX2 

0MY2PD = (GY2PT + OMX2P*OMZ2P* ( C IZZ2-C I XX2 > ) /CIYY2 
1 + AOY2PT/CIYY2 

0MZ2PU = ( GZ2PT + 0MX2P*0MY2P* ( CIXX2-C I YY2 ) ) /C IZZ2 
1 + AOZ2PT/CIZZ2 

C CALCULATE BODY AXIS VELOCITY KATES 

U1PPD = - 0MY1*W1PP + OMZ1+V 1PP + FX1IT/CM1 

1 + AFX1IT/CM1 

V1PPD = - 0MZ1*U1PP + 0MX1*W1PP + FY1IT/CM1 
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+ AFY1IT/CM1 

W1PPD = - 0mX 1*V1PP + 0MY1*U1PP + FZlIT/CMl 
1 + APZ1I1/CM1 

U2PPD = - 0MY2+W2PP + 0MZ2*V2PP + FX2IT/CM2 
1 + AFX2IT/CM2 

V2PP0 = - 0mZ2*U2PP + 0MX2*W2PP + FY2IT/CM2 
1 + AFY2IT/LM2 

W2PF0 = - OmX2*V2PP + OMY2*U2PP + FZ2IT/CM2 
1 + AFZ2IT/CM2 

L StT IN KATES OF CHANGE OF COOKDIMATES AS DYOXS 

DO 1280 J = 1 » NCABLE 
JUMP2 = 3*J + 20 + NREShV 
DYDX ( JUMP2 ) = XBRD(J) 

D YDX ( JUMP2 + 1) = YBRD(J) 

128U DYDX ( JUMP2 + 2) = ZbRD(J) 

UYuX(23)-CTHl 1 *CPSI 1 E*UlPP+ ( SPHI 1*STHT i*CPSl 1E-SPSI 1E*CPHI 1 ) *V1P 
lP+(SPSIlt*SPhIl+CPHll*STHTl*CPSIlE)*WlPP 
uYuX ( 24 ) -CTH T 1 *SPSi 1E*U1PP+ ( CPHI 1*CPSI 1E+SPHI 1*STHT ] *SPSI 1 E ) * V1PP 
l+(CPHll*STHTi*SPSIlE-SPHIl*CPSIlE)*WlPP 
UYUX(2b)=SPHil*CTHTl*VlPP-STHTl*UlPP+CPHIl*CTHTl*wlPP 
UYuX(2b)-CTm?*CPS12E*U2PP+(SPHl2*STHT2*CPSl2F-SpSI2E*CPHT2)*V2PP 
1+(SPSI2E*SPH12+CPHI2*STHT2*CPSI2E) *W2PP 
0YuX(27)-CrHT2*SPSI2E*U2PP+(CPHl2*CPSl2E+SPHI2*STHT2*SPSl2E)*V2PP 
l+(LPtll2*bTHT2*SP5l2E-SPHI2*CPSI2E)*W2PP 
uYUX 128)=SPHI2*CTHT2*V2PP-STHT2*U2PP+CPHI2*CTHT2*W2PP 
RETURlM 
END 
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SUBROUTINE CAbFOR 

DIMENSION P ( 6200 ) »NTE6EK<100) »VAR(6800) >FX1I(20) >FY1I(20) 
»FZ1I (20) t XbR ( 20 ) » ZbR ( 20 ) t XbRU ( 20 ) t YRRD ( 20 ) »ZBRD(20) f 
CABLE (2u ) *SPRK(20) »CDAMP(20) »KP1P2(20) r FORS (20 ) t COEE ( 20 ) 
»rBK(20) 

COMMON VAN 

EQUIVALENCE ( VAR ( 40 1 ) » NTEGER ( 1 ) ) » (VAR(60l) »P( 1 ) ) » (P (461 > * 
FX1K1) ) r (P(481) rFYlim >r<P(501)»FZlI(l) ) r ( P ( 821 ) * XBR ( 1 ) 
) r ( P ( 841 ) r YbR(l) ) » ( P ( 861 ) »ZBR(1) ) » (P(881) t XBRD ( 1 ) ) » (P(901 
) »YbRL)(l> ) » ( P ( 921 ) * ZBRD ( 1 ) ) r (P(1221) rSPRK(l) ) » <P(1241> » 
CDAiviP(l) ) , ( P ( 1261 ) »RP1P2(1) ) » (P ( 1301 ) rFORS ( 1 ) ) r (P(1345) » 
COEE(l) ) * ( P ( 1201 ) » CABLE ( 1 ) ) 

NCAbLE = NTEGER (21) 

COMPUTE INSTANTANEOUS CABLE LENGTH 
DO 180 J = 1 > NCABLE 

RP1P2 ( J ) = SORT (XbR ( J) **2 + YbR(J)**2 + ZBR(J)**2) 

CHECK FOR CARLE SLACK AND ACT ACCORDINGLY 

IF ( RP1P2 ( <J ) - CARLE( J) ) 100»50»50 

COEE(J) = SPRK ( J ) * ( 1 » 0 - CABLE ( J) /RP1P2 (J) ) 

+ CDAMP( J)*(XbP( J)*XBRD( J) + YBR ( J) *YBRD ( «J ) 

+ ZbR ( J ) *ZRKD ( J ) ) / ( RP1P2 ( J ) **2 ) 

FXli(J) = COEE ( J ) *XBR ( J ) 
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FY1KJ) = 
FZ1KJ) = 

COEE(J) *YBR( J) 
COEE ( J) *ZBR ( J) 

100 

GO TO 130 

IF THE CAbLE IS SLACK SET FORCES TO ZERO 
FXII(J) =0.0 

130 

FYll(J) = 
FZli(J) = 
FORS(J) = 

0.0 

0.0 

SQRKFX1I ( J)**2 + FY 1 1 ( vi) **2 + FZ1I(J)**2) 

160 

CONTINUE 



P(1342)=kPlP2(l) 

FCAB(viX=0.0 
UO 300 U=1»NCABLE 
IF (FoRS ( o ) --FCABMX )300»300»301 
bOl FCaBmX=FuRS( J) 

AAa=U 

30 U CONTINUE 

P( 1344)=FCABMX 
P(d99B)=aAA 
P(1340) - KP1P2(2) 
p(134l) - RP1P2(3) 

P ( 1343) 3 KP1P2 ( 4 ) 

RETURN 

END 


4IBFTC SKC/ 

SUBROUTINE FORFUN 

DIMENSION Y ( 100 ) »P(fa200) »NTEGER C100) » VAR (6800 ) 

Common var 

EQUIVALENCE (VAR(l) »Y(1) ) » (VAR(401) »NTEGER ( 1) ) * (VAR(60l). 

1 » P ( 1 ) ) 

EQUIVALENCE (P(128) »AGX1PT) » (P(129) » AGYlPT) r ( P ( 130 ) » AGZ1P 

1 T ) * (P ( 131 ) * AGX2PT )»(P(132)r AGY2PT ) » (P ( 133) t AGZ2PT ) » ( P < 13 

2 4) » AFX1 IT ) » (P(138) t AFY2IT ) r (PC 139) » AFZ2IT) 

3 » (P(135) »AFY1IT) f (P(136) rAFZlIT) » (P(137) »AFX2IT) 
i=NTEGER(l) 

GO T0(10i»102»103) * I 

101 AGX1PT=P(595B)*SIN(P(5999)*Y(1) ) 

AGY1PT=P( 5957) *SIN(P( 5999) *Y( 1) ) 

AGZ1PT=P ( 5956 ) *SIN ( P ( 5999 ) *Y < 1 ) ) 

«FX1IT=P (5946) *SIN(P (5945) *Y( 1) ) 

AF Y 1 1T=P l 5947 ) *SIN ( P ( 5945 ) *Y ( 1 ) ) 
aFZ 1IT=P( 5946) *SIN(P( 5945) *Y( 1) ) 
bO TO 105 

102 HGX2PT=PC595a)*SlN(P(5999)*Y(l) ) 

AGY2PT=P ( 5957 ) *SIN ( P ( 5999 ) *Y ( 1 ) ) 

AGZ2PT=P ( 5956 ) *SIN ( P ( 5999 ) *Y ( 1 ) ) 

AFX2 IT=P ( 5946 ) *S IN < P ( 5945 ) * Y ( 1 ) ) 

AFY2IT=P 15947) *SIN(P(5945)*Y(1) ) 

AFZ2IT=P( 5946) *SIN(P( 5945) *Y( 1) ) 

103 CONTINUE 
RETURN 
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SIbFTC SKCtf 

SUBROUTINE OUTPUT 

DIMENSION Y( 100) *DYDX(100) *P(b2u0) * INTEGER (100) * VAR (6800) * 

1 RtC (14) 

COMMON/Fi_M/OM*ON*TwM*TwN*THM*THN*FOM*FON*FIM*FIN*SIM*SIN*SEM*SEM* 
1EIM»EIN» aNIM*XNIi\i» TEM*TEN»ELM*ELN*TWEM*TWEN*TIEMPO*THIM*THIN*KUT# 
2iTHB0P * MOV IE * IGRoP * Phi * THETA * NUMI * NUM2 * NUM3 * KWHICH 
COMMON VAN * KNT * KFST 

EQUIVALENCE (VAR(l) * Y ( 1 ) ) * (VAK(lOl) * OYDX (1) ) ► (VAR (A 01) * 

1 NTEGER(l) ) * ( VAR ( 60 1 ) » P ( 1 ) ) * (NTEGER (24) »NTSKIP) * (NTEGER (41 

2 ) * LPR INT ) * (NTEGER(42) *NLlNE) » (NTEGER (43) *NSKIP) * (NTEGER ( 

3 44 ) * NPAGE ) 

4 * ( P ( b991 ) * YCG) * (P(b99G) »XCG) * (P(b989) *ZC6) * (P(5988) *XBR2CG) * (P(59 
5a7) * YBR20G) * (P(598d) »ZBR2CG) * (P(b965) *PSICAP) * (P(5984) rPHlCAP) 

EQU I VALENCE (P(S970 ) *EX2) » (P(S969) *EY2) * (P(S9b8) *EZ2) 

SOrt XI- (Y (26)— Y (23) ) **2+ ( Y ( 28 ) - Y ( 25 ) ) **2 
UR 12 -SGR I (SQAX1+(Y(2 n-Y(24) )**2) 

CG=(P(22)*bRl2)/(P(6)+P<22) ) 

C r COORD. OF BODY ICG IN INERTIAL PRIME SYSTEM 

YCG=C6*( t (24)-Y(27) J/BR12 
rtLNAX=SQr<T(SQAXl) 

ALAXI-SQkT ( A bS ( Cb**2-YCG**2 ) ) 

C X COORD. OF BODY 1 CG IN INERTIAL PRIME SYSTEM 

XCb=ALAXi*(Y(23)-Y(26) J/ALNAX 
CALL DVChK (KOOOFX) 

GO TO ( 200 * 201 ) *KGO0FX 
20 U XCG=0.0 

C l COORD. OF BODY I CG IN INERTIAL PRIME SYSTEM 

201 2CG=ALAXi*(Y(2b)-Y(2h) J/ALNAX 
CALL DVChK (KOOOFX) 

bO TO (202*203) *K000FX 

202 2Cb = 0 • 0 

C X* Y*Z*COoRDS.OF RIGID BODY CG IN INERTIAL SYSTEM 

200 XBK2CG=Y (23) -XCG 
YBK2CG=Y(24)“YCG 
ZBk2CG=Y(26)-ZCG 
eX2~Y ( 26 ) “XBK2CG 
EY2=Y(27)“YBR2CG 
EZ2=Y (28)-2BK2CG 
AUXCA1=SoRT(EX2**2+EY2**2) 

HSiCAP=AhTN(EY2*EX2) 

IF (PSICrtP)204*205*20S 

204 PS I CAP— 6 ,28318+PSICAP 

20 5 P(b973)=(-1.0)*EZ2/AUXCA1 
SHIPT=P(b973) 

PHI CAP- AT An ( SPIPT ) 

CPHlCP=CbS(PHlCAP) 

CPSICP=CoS(PSICAP) 

SPSICP=SiNE(PSICAP) 

SPHlCP=SiNE(PHICrtP) 
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P ( 5963 ) =UYDX ( 25) + ( UYDX ( 28 ) -DYDX ( 25 ) ) *CG/BR12 
P ( o9b2 ) =U YDX ( 24 ) + ( D YDX ( 27 ) -D YDX ( 24 ) ) *CG/BR12 
P ( 59bl ) =uYDX ( 23) + ( DYDX ( 2b) -DYDX ( 23 ) ) *CG/BR12 
YDZRCG=DYDX ( 27 ) -P ( £>962 ) 

XDzRCG=:DYDX(26)-P(5961) 

ZDZRCG=0 YDX ( eLS ) -P ( 5963 ) 

aLuAPE=( YD2RCG*SPS1CP+XD2RCG*CPSICP)*SPHICP+ZD2RCG*CPHICP 
HFbAHF=YL2HCG*CPSICP-XD2RCG*SPSlCP 
P ( 3960 ) =kESAPE/AUXCA 1 
H(5959J=bLUAPF 
1 / ( bG-RRl2 ) 
bL.UAPM=-bLUAPF 
P(5950.)=ARTN(RLUAPM»RESAPE) 

P ( 39*+9) =ARCOS ( CPHlCP*COS (P ( 5950 ) ) ) 

P ( 5951 ) =bQRT (BLU aPM**2+RESAPE**2 ) / ( BR12-CG) 

5FuA=CPHiCP*CPSICP 
bPLB=-bPSICP 
bMt-C=SPHiCP*CPSlCP 
bMLD=CPHiCP*bPSICP 
5 HlF=SPH i CP*SPS I CP 
bKi-G--bPhlCP 
m 2=P l 39 ) *P ( 5978 ) 
l>2— P ( 39 ) *P 1 5976) 

L?--P ( 38 i 

b2=P 1 4b ) *P t 3b ) *P ( 5978) -P ( 5976 ) *P ( 47 ) 

1 2=P ( 47 ) *P ( 5978 ) +P C 46) *P ( 38 ) *P ( 5976 ) 

62=P(46)*P<39) 

b2=P ( 597b ) *P ( 46 ) +P ( 47) *P < 36 ) *P ( 5978 ) 

H2-P ( 47 ) *P ( 36 ) *P ( 5976 ) -P l 46 ) *P ( 5978 ) 

AI2=P(47)*P(39) 

A1=P ( 36 ) *P ( 5979) 
dl=P(3b)*P(5977) 

C1-”P ( 35 ) 

ul -P ( 44 )*P(35)*P( 5979 ) -P ( 5977 ) *P ( 45 ) 

1 1 =P < 45 ) *P ( 5979 ) +P ( 44 ) *P ( 35 ) *P ( 5977 ) 

P 1 = P l 44 ) *P ( 3b ) 

b 1=P ( 5977 ) *P ( 44 ) +P ( 45 ) *P ( 35 ) *P ( 5979 ) 

H 1=P ( 45 ) *P ( 35 ) *P ( 5977 ) -p ( 44 ) *P ( 5979 ) 

All=P(45)*P(36) 

p(39o2)=hRTN( (A2*SMLB+B2*CPSICP) t ( A2*SMLA+B2*5MLU+C?*SMLG) > 

P(5963)=mRTN( (A 1*SMLB+B1*CPSICP) r ( A1*SMLA+B1*SMLD+C1*SMLG) ) 

COMMll=Gi*SMLC+Hl*SMLF+AIl*CPHICP 

C0MSM2=D2*bMLC+E2*SMLF+F2*CPHICP 

C.OMSMl=Di*SMLC+El*SMLF+Fl*CPHICP 

C0MM22=G2*SMLC+H2*SMLF+AX2*CPHICP 

P(5975)=ARTN( C0MSM2 » C0MM22 ) 

P(5974)=ARTN( COMSM1 * COMM 11 ) 

P(3980)=mRTN( ( (- 1.0)*(A2*SMLC+B2*SMLF+C2*CPHICP) ) »SQRT( 
lbM2**2+CuMM22**2) ) 

r(596l)=MRTN( ( ( -1 . 0 ) * ( A1*SMLC+B1*SMLF+C 1*CPHICP) ) t SORT ( C0M$M1 
1+C0MM11**2) ) 

M0VlE=NTc_GER(26) 

CVT=57.2956 


COM 

**2 
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IF(KFSI) 1 » 3 f 1 

3 0M=P(598i) 

UN=0M 

TWM=P(b9/4) 

TWN-TWM 

IHw=P(59to3) 

THIM-1HM 

i-OM=P(b9bO) 
b ON— (- OM 
FlM=P(b9/5) 
t- Ih=F IM 
bIIW=P(b9b2) 
blN— blM 
bt|vi=Y (20) 
bE i\i=Y (20) 
tlH-P ( b9b4 ) 
t IN=tI,M 
xMiM=P(5^8b) 

XNIN=XNIm 
T F M— Y ( 2 3 J 
ie:n=y<23) 

tLM=Y (24) 
c.Lin=Y (24) 

I ViitM=Y (2b) 
rwtN=Y(2b) 

I FUNi— P 1 1342 ) 

!HlN=P(lo42) 

NFbT=l 
bO TO 66 

1 IF* (P ( 598 i ) -O n ) b » b » 4 

4 uM-P ( 598 J. ) 

b IF (P ( 6981 ) -ON) 6» 7 r 7 

6 UN— P ( b98i ) 

7 IF(P(S974)-TWM)9»9»8 
a rv,'M=P(b974) 

9 iF(P(b974)-TrtK) 11H2» 12 

11 TWN=Plb974> 

12 IF ( P ( 5983 ) — THN" ) 14 » 14 » 13 

13 rHM=P(59b3) 

14 1F(P(5983)-THM)15»16»16 
lb lHn=P(59b3) 

lb IF ( P l 598u ) — FOM ) Id r 18 » 17 
17 FOM=P(59bO) 
lb IF ( P ( b98u ) —FOP) 19 » 21 * 21 
19 F0N=P(b9sG) 

21 IF ( P ( 597 b ) —FlN* )23»23»22 

22 FlM=P(b975> 

23 IF(P(b97b)-FIf024»25r2b 

24 FlN=P(b975) 

2b IF ( P ( 5982 ) -SIP ) 2.7 1 27 1 2b 
2b SlM=P(59b2) 

27 IF ( P ( b982 ) “SIN ) 28 » 29 r 29 
2b SlN=P(b9o2) 



29 IF l Y ( 2G ) -StM ) 32 » 32 * 31 

31 SEM-Y ( 20 ) 

32 IF ( Y (20) -SEN) 33 » 34 » 34 
35 >Ew=Y(20) 

34 1F(P(59H4)-E1M)36»56» 3b 
3b tlM=P(59b4) 

3b IF IF ( 5984 ) -E1N> 37 » 38» 38 
37 c Iim-F ( b9b4 ) 

30 IF l P 1 598b )~X nIM)41»41* 39 
39 aNXiV-P ( 598b) 

41 IF ( P ( 598b ) -XNIN) 42 » 45 » 43 
4 c. aM1N=P(598o) 

43 lF(Y(23)-TEM)4b»45»44 

44 IEm=Y<23) 

4b IF (Y ( 23) —TEN) 46 » 47,* 47 
4b I t'N—Y ( 23 ) 

47 iF( Y(24)-ELM)49*49*48 
4b ELM = Y ( 24 i 

49 IF l Y 1 24 ) -ELN) 51*52* 52 

51 ELN=Y(24) 

52 iF(Y(25)-TwEM)54*54*53 

53 lWEM=Y(2b) 

54 IF(Y(2b)-TwEN)b5*56*56 
5b (WEN=Y(2b) 

5b IF(P(1542)-THIM)58*58*57 
57 IH1M=PU042) 

5b IF ( P ( 134c! ) -ThlN ) 59 * 66 * b6 
59 IH1N=P(1342) 
bb COimT INUE 
KNl=KNT+'i 
KEC(1)=YU) 

KE\.(2)=CVT*P(598i) 

KEC<3)=CvT*P(5974) 

kEC(4)=C\/T*P(5983) 

KEt (b)=C(/T*P(5980> 
riFt (b)=CVT*P(597b) 

Kfct (7)=CvT*P(59B2> 

HEC(d)=CVT»Y(20) 
kEC(9)=CVT*P(5984) 
kEC ( 10 ) =eVT*P ( 5985) 

HEC(il)=T (23) 

KEC( 12)=Y (24) 
kEC(13)=Y(25) 
hEC ( 14 ) =P ( 1342 ) 

WRITE (9) (NEC(I) *1=1*14) 
iF(POVIE)3U0»]0r50U 
30U x I 1=XCG+P ( 56 ) *P ( 5979 ) *P ( 5967 ) 
YU=YCb+P(56)*P(5977)*P(5967) 

Zll=2C6-P(55)*P(b9b7) 

a J i=XC6+ ( P ( 44 ) *P ( 35 ) *P ( 5979 ) -P ( 5977 ) *P ( 45 ) ) *P ( 5967 ) 
t Jl=YC6+(P(4b) *P( 5979) +P(44)*P( 3b) *P( 5977) )*P(59&7) 
/L Ji=2C6+P ( 44 ) *P ( 36 ) *P ( 5967 ) 

XK1=XCG+ ( P ( 5977 ) *p ( 44 ) +P ( 4b ) *P ( 3b ) *P ( 5979 ) ) *P ( 5967 ) 
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YK 1= YC6+ ( P ( 4b ) *P ( 3b ) *P l 5977 ) -P { 44 ) *P ( 5979 ) > *P ( 5967 ) 
ZKl=ZC6+P(45)*P(56)*P(5 c 'b7) 

Xl2=tX2+P ( 59) *P ( 5978 ) *P ( 59b7 ) 

Y I2=EY2+P ( 59 ) *P ( 5976 ) *P ( 5967 ) 

2 l2=tZ2-P ( 38 ) *P < 59o7 ) 

xJ2=EX2+(P(4b)*P(38)*P(597b)-P(5976)*P(47) )*P(5967) 
YJ2=EY2+(P(47)*P(5978)+P(46)*P(36)*P(5976) )*P(59b7) 

2 J2=tZ2+P ( 46 ) *P ( 39) *P ( 5967 ) 

xk2=EX2+ l P ( 5976 ) *P ( 46 ) +P ( 47 ) *P ( 3b ) *P ( 5978 ) ) *P ( 5967 ) 
YK2=EY2+(P(47)*P(38)*P(5976)-P(46)*P(5978) )*P(59b7) 

2Kz=tZ2+P ( 47 ) *P ( 39 ) *P ( 5967 ) 

WRITE ( lx ) P ( 3990 ) » P ( 5991 ) rP(5989) * XI 1 » Y 1 1 » 21 1 * XJ1 » Y J1 » ZJ1 » XK1 » YK 1 * 
1ZK1 1 p ( 59 / 0 ) t P ( 5969 ) »P(59&8) * X 12 * Y 12 » 2 12 » X 02 * Y J2 » Z J2* XK2 * YK2 » ZK2» P ( 
25 970) r P ( 5969 ) t Y(l) 

XlU=xC6+Ml*P(5967) 

Y1u=YCg+d1*P(5967) 

Z1u=ZCg+c1*P(5967) 

Xl x=aC6-m 1*P(5967)-P(5965)*(01+ 6l) 

Yll=YCG-bl*P(5967)-P(5965)*(tl+ Hi) 

21 l=ZC6-Ul *P ( 5967) -P( 5965) *(F1+Ail) 

Xl2-XC6”Ml*Pl 5967 ) — P ( 5965 ) * (61- Ul) 

Y12=YC6-d 1*P ( 5967 )-P( 5965)* (HI- El) 

Zl2=ZC6-6l*P ( 5967 ) -P (6966 ) * ( AI 1-F 1 ) 

Xl3=xC6+Al*P ( 5967) -P ( 5965 ) * ( 61 -Ul ) 

Yl3=YC6+bl*P ( 5967 ) -P ( 5965) * ( HI -El) 

Zl3=ZC6+6 1*P ( 5967 ) -P ( 5965 ) * ( AI 1-6 1 ) 

X14-XC6+m 1*P ( 5967 ) -P ( 5965 ) * ( 61 +Ul ) 

Y14=fC6+ol*P(5967)-P(5965)*(Hl +E1) 

Zl4=ZC6+U 1*P ( 5967 ) -P ( 5965 ) * ( A 1 1+6 1 ) 
a15-xC 6+A1*P ( 5967 ) — P ( 5965 ) * ( Ul— 61) 

Y 15=YCG+bl*P ( 5967 ) -P ( 5965)* (£1- Hi) 

2 1 5=ZC6+6 1 *P ( 5967 ) -P ( 5965 ) * ( 6 1 -A 1 1 ) 

Xlb— xC6“Al*P l 5967) -P ( 5965) * (Ul- Gl ) 

Ylb=YC6-ul*P( 5967) -P( 5965)* (El- Hi) 

2 lb=2CG-b 1*P ( 5967 ) -P ( 5965 ) * ( F 1-A I 1 ) 

X17=XCG-M1*P( 5967) +P( 5965) *(D1+ 61) 

Y1 /=YCG-bl*P(5967)+P(5965)*(El+ HI) 

Z 1 7-ZCG— 6l*P ( 5967 ) +P ( 5965 ) * ( F 1+A 1 1 ) 

Xlb=XC6+Al*P( 5967) +P (5965) *(01+ 61) 

Ylb=YCG+bl*P ( 5967) +P (5965) *(E1+ HI) 

Zlb=zC6+6 1*P ( 5967 ) +P ( 5965 ) * ( 6 1+A 11 ) 

X2U=EX2-h2*P(5966) 

Y2U=EY2-d2*P(596b) 

Z2u=EZ2“U2*P(596b) 

X21=EX2+A2*P ( 5966 ) +P ( 5964 ) * ( 02- 62 ) 
f 21=EY2+b2*P ( 5966 ) +P ( 5964 ) * ( E2- H2 ) 

Z21=EZ2+C2*P ( 5966 ) +P ( 5969 ) * < F2-A 12 ) 

X22=EX2+A2*P ( 5966 ) -P ( 5964 ) * < 02+ 62 ) 

Y22=EY2+b2*P ( 5966 ) -P ( 5964 ) * ( E2+ H2 ) 

Z22=EZ2+t 2*P ( 5966 ) -P ( 5964 ) * ( F2+ A 12 ) 

X23=EX2-A2*P ( 5966 ) -P ( 5964 ) * ( 02+ 62 ) 
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Y 23=EY2~b2*P ( 5966 ) -P ( 5964 ) * ( E2+ H2 ) 

Z2 3=EZ2-C2*P ( 5966 > -P ( 5964 ) * ( F2+AI 2) 

X24=EX2-A2*P ( 596b ) -P ( 6964 ) * ( G2- D2 ) 
t 24=EY2-b2*P ( 596b ) -P ( 5964) * ( H2- E2 ) 

ZP4=EZ2-C2*P l 5966 > -P l 5964 ) * ( A I2-F2) 
aP5=EX2-a2*P (5966) +P ( 5964) *(02+ 62) 

Y25-EY 2~b2*P 1 5966 ) +P ( 5964 ) * ( E2+ H2) 

Z2d=EZ 2-C2*P ( 5966) +P ( 5964) * (F2+A12) 
x2b=EX2+A2*P ( 5966 ) +P ( 5964 ) * ( 02+ 62) 

Y 2b=EY2+d2*P ( 5966 ) +P ( 6964 ) * ( E2+ H2) 

Z2b=£Z2+c2*P ( 5966) +P ( 5964) * (F2+A12) 
a 27=EX2+A2*P ( 5966 ) +P ( 5964 ) * < G2- 02) 
rP7=tY2+u2*P(5966)+P(5964)*(H2- E2) 

Z27=tZP+c2*P ( 5966) +P ( 5964 ) * ( AI2-F2) 

X2b=tXP-A2*P(596b)+P(5964)*(G2 -02) 

Y 2b=EY2-D2*P ( 5966 ) +P ( 5964 ) * ( H2 -E2 ) 

ZP6=EZ2-t2*P l 5966 ) +P ( 5964 ) * ( AI2-F2) 

wHiTt (13)XlU*X13*Xl4*X15»Xl8*X17*X12*Xll*Xl6*X20*X23*X24*X25*XP8* 
1x27* X2P *A2l*X26*Y10*Y13*Y14*Y15*Y18*Y17*Y12*Yll*Y16»Y20*Y23* Y24* 

2Y25*Y28*)27»Y22*Y21*Y2b*Z10*Z13*Z14*Zl5*Z18*Z17*Z12*Zll*Z16* Z20* 

3Z23* Z24 r 425 *Z28*Z27* Z22* Z21 * Z26 * P ( 5970 ) * P ( 5969) » Y ( 1 ) 

10 I F ( Y ( 1 ) - P<2) )20*50*50 

C DETERMINE WHETHER OR NOT PRINT ON THIS INTEGRATION STEP 

PU IF ( NSKIP - NTSKIP)30*50*50 

3u NSKIP = NSKIP + 1 

GO TO 150 

C DETERMINE IF A NEW PAGE IS REQUIRED FOR PRINTING RESULTS 

' 5u IFlNLINE - 52)90*60*61 

6u CALL LASOUT 

61 NPAbF = NPAGF + 1 

CALL PAbEHD 
NLlNE = 0 

9U CALL OUT AID 

NLINF = NL.INE + LPRIMT + 1 
NSKIP = 0 

IF ( Y ( 1 ) - P(2) ) 150* 130* 130 
c COMPLETE PRINTOUT AND GO TO FILM SUBROUTINES 

130 CALL LASOUT 

TIlMR0=Y (1) 

PHI=P(59/2) 

THtTA=P(b971) 

1GK0P=NTEGER t2b) 

ITHB0P=NTEGER(27) 

NUM1=600 
NUM3=600 
KWHlCH=l 
ENU FILE 11 

ENU FILE 9 

lNl) FILE 13 

150 RETURN 

END 


ilttFTC SRCV 
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SUBROUTINE OUTAID 

DIMENSION Y (100) *DYDX(1U0) >P(6200) t NTEGER ( 100 ) * VAR (6800) » 

1 XBR120) f YbR(20) r ZBR ( 20 ) » OUTP1 ( 13 ) * CARLE ( 20 ) * RP1P2 ( 20 ) » 

2 0UTP2 ( 57» 13) * 0UTP4 ( 57 r 13) 

0 I MENS I ON MGUT(2u) 

common vak 

EQUIVALENCE ( VAR ( 1 ) »Y(1) ) r (VAK (101) »DYDX(1) ) » (VAR (601) » 

1 P ( 1 ) ) t ( VAR ( 401 ) * NTEGER ( 1 ) ) » 

2 (NTEGER ( 42 ) . NLINE ) # (Y(8) t THT1 ) » (Y(9) »PHI1) t ( Y (18) »PHI2) t 

3 ( Y ( 17 ) » THT2 ) » (Y(20) » THTbR ) » (DYDX(20) rTHTBRD) » (Y(21) r PHlb 

4 R ) t ( DYDX ( 21 ) * PHIBRD ) » (Y(22) »PSIBR) * (DYDX(22) »PS1RRD) » 

5 (P(5B) » OMZ1 ) » (P(59) » 0MX2 ) » ( P < 60 ) t OMY2 ) » (P(61) »0MZ2) » 

6 (P(b6) f OMX1 ) * ( P ( 57) # OMY 1 ) » (P(fa41 ) t YbR ( 1 ) ) » (P ( 821 ) r XBR ( 1 ) ) 

7 » (P ( 861 ) * ZbR ( 1 ) ) » (P ( 941 ) » OUTP1 (1))»(P( 1201 ) » CABLE ( 1 ) ) » 

8 ( P ( 1261 ) t RP1P2 ( 1 ) ) » (P(2151) »0UTP4(1) ) * (P(103) »K) » (P(1401 

9 ) »OUTP2 ( 1 ) ) » ( NTEGER ( 21 ) » NCABLE) 

C CALCULATE QUANTITIES TO BE PRINTED ON PAGE 1 

C TO PRINT OUT THE GRAPH VAR I ABLES » MAKE THE HOLLOWING CHANGES IN THE OUT 
C Pl( ) CAKDS-CHANGE THT 1 TO P(5981)»PHll TO P(5974)»THT2 TO P(5980)»PHl 
C 2 TO P ( 597b ) t XBR ( 1 ) TO P ( 5983 ) *57 . 29B78 • YPk ( 1 ) TO P ( S982 ) *57.29578* ZBR 
C (1) TO P ( 598 a ) *57 . 2957b 

0UTPK1) = Y ( 1 ) 

OUTP1 (2) = THT1 *57.2957795 

OUTP1 ( 3 ) =P ( 5960 ) *57 • 2957795 

OUTP1 (4) = Phil *57 . 29b779b 

0UTPK5) '= THT 2 *57. 295779b 
OUTPl (6) =P( 5959) *57.2957795 

OUTP1 (7)=PH12 *57.2957795 

OUTP1 (8) = THTBR+57. 2957795 

0UTPK9J = PSIBR*57. 2957795 
OUTPl(lO) = PHlbR*57. 2957795 
OUTPK11) = XBR ( 1 ) 

0UTPK12) = YBR(l) 

OUTP1 ( 13) = ZBR ( 1 ) 

C WRITE COLUMN HEADINGS IF NEW PAGE 

IF (NLINE) 15u> 150 » 170 
15U WRITE ( 6 * 16(J ) 

DO 611 L-l » 57 
DO 612 1=1*13 
OUT P2 ( L * 1 ) =0 . 0 
bl2 CONTINUE 
611 CONTINUE 

DO 613 L=l*57 
uO 614 1=1*13 
0UTP4 (L*l)=Q.O 
ol4 CONTINUE 
bl3 CONTINUE 


160 FORMAT ( // 127h 

TIME 

THETA 1 

PS1RBD 

PHI 1 

THETA2 

1 THETRBD PH 12 

THETAB 

PSIB 

PHIB 

XBR 1 1 ) 

YBR ( 1 ) 

2 ZBR ( 1 ) / 124H 

SEC 

DEG 

deg/sec 

DEG 

DEG 

3DEG/SEC DEG 

DEG 

DEG 

deg 

IN 

IN 


4 IN//) 
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C 


WHITE OUT DATA ON FIRST PAGE 
170 WHITE (6»180) (OUTPl(I) »I=1,13> 

18U FORMAT(3(F11.4,F9.4»F10.4) » Fll . 4 , 3F9. 3) 

MOUTPR=0 

UO 500 I=2,NC'ABLE 
IF ( RP1P2 ( I ) -CABLE ( I ) ) 502 , 502 » 501 
5U1 MGUT(I)=0 
bO TO 500 
502 MGUT(I)=1 
MOUTHR=l 
500 CONTINUE 

IF l MQUTP k ) 504,504,50b 
505 WRITE (6,503) 

505 FORMAT HbH SLACK CABLES ARE > 

UO 506 1=2 , NCABLE 

IF l MGDT(I) )506,506,507 

507 bO TO ( 511, 512, 513, 514, 515 » 516, 517, 518» 519 ,520, 521 » 522 »523» 524 » 525 
1 ,526,527,528,529,530) » I 

511 WRITE ( 6 , 531 ) I 

531 FORMAT ( 1H+19X » 1 1) 

GO TO 50b 

512 WRITE ( 6 , 532 ) I 

532 FORMAT ( 1F1+21X, II ) 
bO TO 50b 

513 WRITE (6,553)1 

533 FORMAT Un+-23X» ID 
bO TO 50b 

514 WRITE (6,534)1 

534 FORMAT (1B+25X, ID 
GO TO 50b 

515 WRITE (6,535)1 

535 FORMAT (lh+27X» ID 
bO TO 50b 

51b wRITE (6,536)1 
53b FORMAT (lri+29X,Il) 
oO TO 50b 

517 WRITE (6,537)1 
537 FORMAT (1N+31X, ID 
bO TO 50b 

51b WRITE (6,538)1 
53b FORMAT (1H+33X, ID 
bO TO 50b 
519 wRITE (6,539)1 

539 FORMAT ( 1F1+35X, 1 1 ) 

GO TO 50b 

52U WRITE (6,540)1 

540 FORMAT ( lFf+37X , 12 ) 

GO TO 50b 

521 wRITE (6,541)1 

541 FORMAT (lh+40X»I2) 

GO TO 50b 

522 WRITE (6,542)1 


46 



542 FORMAT ( lh+43X . 12 ) 
bO TO 50b 
b2o WRITE ( 6 » 543) I 
545 FORMAT ( 1H+46X » 12 ) 
bO TO bOb 
b24 wRiTE ( 6 . 544 ) I 

544 t-0RMAT(lH+49X» 12) 
bO TO 50b 

b2b wRiTE (6.545)1 

545 FORMAT! 1F1+D2X. 12) 
<30 To bOb 

52b WRITE (6.546)1 
54b FORMAT ( 1FI+55X. X2) 
GO To 50t> 

527 WRITE (6.547) I 
547 t-ORMAT ( ln+58X» 12) 
bO TO bOo 
52b wRiTt (6.548)1 
54b FORMAT ( ln+ol/. » 12) 
jO To bOu 

529 WRITE (6.549)1 
549 t-ORMAT ( lh+o4X » 12 ) 


bO TO bOb 
b3U wRiTE (6.550)1 
bbu FORMAT ( lh+b7X. 12) 
bOb CONTINUE 

b04 IF ( MOUTPr ) 509 » 509 » 51 U 
51U NL 1 NE-NL 1 NE+ 1 

CALCULATE DATA FOR SECOND PAGE 
509 K = NLInE + 1 

0UTF2 ( K . 1 ) = Y(l) 

U(J I P2 ( K . b ) — P ( 5950 ) *57 . 2957795 
OU I P2(KO)=P (5949) *5 7.2957795 

0UTP2 ( K . 4 ) = P( 5951) *57. 2957795 

0UTF2(K.5)=THTBRD* 57.2957795 
0UTP2(K.6)=PSIBRD* 57.2957795 
0UTH2(K.7)=PHIBRD* 57.2957795 
0UTP2 ( K . 8 ) =P(5990) 

0UTP2 ( K . 9 ) = P ( 5991 ) 

0UTP2 (K . 1U ) =P<5989) 


C 

C TO 
C 4 ( 


PRINT OUT 


0UTP2(K. 11) 
0UTP2(K. 12) 
0UTP2 (K . 15) 
CALCULATE 
THE GRAPH 


) CARUS-CHANGt PQ34U) 
0UTP4(K. 1) 
0UTP4 ( K » 2 ) 
0UTP4 ( K » 3 ) 
0UTP4 ( K » 4 ) 
OUTP4 (K . 5) 
QUTP4 ( K . 6 ) 
0UTP4 ( K . 7 ) 


=P(5988) 
=P(5987) 
=P(5986) 

DATA FOR THIRD 
VARIABLES. MAKE 


PAGE 

the following 


TO P(5985)*57. 29578 
: Y ( 1 ) 

: OMX 1*57 .2957795 
: 0MY1*57. 2957795 
: 0MZ1+57. 2957795 
: 0MX2*57. 2957795 
: 0MY2*57. 2957795 
: 0MZ2*57. 2957795 


change in the outp 
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L 


4 lU 

VII 


0UTP4 ( K » H ) 
0UTP4 ( K » 9 ) 
0UTP4 ( K * 1 0 ) 
0UTP4 ( K » 1 x ) 
0UTP4(K# l^) 
0UTP4(K» lo)=P(S996) 


P(1340) 
P ( 1341 ) 
P ( 1342 ) 
P( 1343) 
P ( 1344 ) 


WHITE OUT PAGtS 2 AND 


IMfMl.iNt - b? ) 420 » 410 

Call LAbOUT 

RETURN 


EnO 


3 IF at end of page \ 
410 


x 1 : ! r It ShC i 0 


SUBROUTINE LAbOUT 

DIMLNbluN P ( 620 0 ) » NTEGEH ( lOU ) t OIJTP2 (57*13) » 0uTP4(57*l3 
1 ) t VmI- (6bOU ) 

CuMjv.ON V AH 

EwU i V aLENoE (VAR ( 40 1 ) * NTeGER C 1 ) ) » (VAR(60l) *P<1) ) * 

1 (NTlGLR ( 44) rNPAGE) » (P( 1401) *OuTP2( 1 ) ) * (P(2151) *0UTP4( 1) ) * 

2 ( P ( 1 n 3 ) »K) 

'JkIif our Page 2 
NPAbF = NPAbF + 1 
CALL PAbEhU 
WKllt- (fa» 44 U ) 

44 u r OKiVaT ( // 123u TIME GAMMA ALPHA I A VB2 THET ABD P 

lblbU Pi-iIl,!' XCb YCG ZCG XdR2CG YbR2C6 ZB 

2k2uG/ 1 23n SEC uEG DlG DEG/SEC DEG/SEC DFG/SE 

3C UlG/SeC INCHES INCHES INChFS INCHES INCHES INCHES 
4 //) 

Ob 4AU L = 1 » K 

W k I IE (b*47u) ( OUTP2 ( L * I ) * 1 = 1 * 13) 

47 u pOHMaTIFi1.3*F9.4*P9.1*F11.3»2F9.4*F9.3*F11.4»F10.4»FB.2*3F9.2) 

4 rt u CONTINUE 

WHITE OUT PaGE 3 
NPAbF = NPAOE + 1 
CaLL PAbEllD 
WHITE (b*b2u) 

b2u K'HPaTF// 127h TIME OMEGAXl OMEGA Y 1 OMEGAZ1 

lot* c.GaY 2 OMEbA/2 RP1P2 (2) HP1P2(3) HP1P2(1) RP1P2(4) 

2 ChHLE/112H SEC OEG/SEC UEG/SEC DEG/SEC 

3t_G/SbC uEb/SFC In IN IN IN 

DO b6U L = i*K 

WHI IF ( b ► b50 ) (0UTP4(L»I) *1=1*13) 
bbu FORMAT (Eli. 3» 2 (Fll. 3»F9.4»F9.4) * 5F1 U . 3 * F6. 0 ) 

bbu CONIINUE 

RtTUPN 
EnQ 


0MEGAX2 
FCABLEMAX 
DEG/SEC D 
Lb//) 


bloFTC SHCil 
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FUNCTION SINE(X) 

SInE=S1N(X) 

RETURN 

END 


SlBFlC SrC12 

FUNCTION AKTim ( X» Y ) 
ARTN=ATAim21X# Y) 
kETUKN 
END 


C REMOVE THE Rt-ST uF THE SOURCE DECK IF S-C4020 NOT AVAILABLE 
iOKlblw ACE 

VliiFTC SrFU 

SUdROU I InE Fll. M ( OM » ON * TWM » TWN » THni » THN » F DM » FON »FlM»FIN»SlM»SIN> SFm * 
1 SEN * t IM » t-lN » XMIM t XNIN * TEM t TEN » ELM » ELN » 1 WEM » TW^N * T IEMPO *THlM»THINrK 
2 NT * I THbOF » MO V IE » 1GR0P» PHI * THETA » NUM1 » NUM2* NUM3 »K WHICH) 

U I MENS i On vAMbBOO) 

uImEnSIOn BEL ( 26 ) fX(iB) » Y(18) »Z(18) 

COMMON VAR t KNT t KFST » L 

Lsn-'si 

I\=1 

L=U 

REWIND li 
REWIND lo 
IF ( I CROP J4lrll»41 
41 REWIND 9 
R=i 
L=L+i 

GO TO (l»2»3»4»5»6f7rB»9rl0rll»l2) »L 

1 XL=0.0 
XR^TIEMPu 
YB— ON*CVI 
YT=OM*CV I 

CALL CAR) (XL»XR»YBrYI ) 

CAuL APRimTV (u»— 12»— lb» 16HTHET ASTI DEGREES * 12 » 650 ) 

M=L+1 
GO TO 33 

2 YB=TwN*CVT 
YT=TrtM*CvT 

CALL CAR 1 (XLf XR> YBf Y i ) 

CALL APRNTV(U»-12»-I4f 14HPHIST1 DEGREES * 12 » 650 ) 

M=L+1 
GO TO 33 

3 YB=THN*CVT 
YT=THM*CVT 

CALL CAR) (XL>XR» YB>YT) 

CAUL APRNTV(U»-I2r-14»14HPSIST1 DEGREES t 12 # 650 ) 

M=L+1 
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uO Tu J>3 
h YFs-F ON*C v r 
If T-FOM*C V'T 

cAlL CAR] (xL»VR» YR» Y I ) 

CALL AHRijT V (ur-ler»-l»:» l6HTHETASTk DEbPfcES » 12 » AbO ) 
to— L+i 
uO TU o3 
s YB-Fju.*CvT 
i T-F I '+CvT 

C AlL C aR i ( xL f V R» Y B t Y t ) 

CALL ApR H Tv(l/»-l«:r-li < ri.4HPHlST2 bEGPtFS * 12 1 6 S 0 ) 

i-Fl+i 

CO Tu 53 

0 j F— S i f !*C m T 
rT-=S* *C\j T 

CALL CAR i .( aL» XK» YH » YT ) 

CALL APRimTv/ (U f-I2f-14»14HPSIST2 uE6Rfc.ES 1 1 2 1 6b0 ) 

to=L+i 

uU Tu 

1 iF l 1 1 1 luOp ) yo » 41 » so 
uu COImT lNUE 

YR=SfcN*CvT 
r T-St' 1 *CvT 

call CAR l (XL» XK» YB» YT ) 

CALL APRimTV ( U t —lid t -1 i+ f 14FTTHET4B DEGREES » 12 # 6bH ) 

to-L+± 

uO TU 

a yh-t i >*CvT 

f t — 1 1 f“' * c v r 

CAlL Car I ( XL» XKHHi Y } ) 

CAlL APRimTv ( U » - 12 » -14 » 14FITFIE1 Rb uEbRtFS * 12 » 6S0 ) 
1*1—1—+ i 

uO Tu o3 
9 \ P — X !M 1 l\i*CV I 
Y T = XlM I |V|*CV J 

c A LL CAR I (XL»XR»YH»Y1 ) 

call APRimTY/ (u#“l2»—l4> 14Fi PSIRw UEGRtF S » 12 * 6S0 ) 

to-L + 1 

GO TO 53 
lu YF-TnlW 
YT-TnlM 

CALL CAR] ( XL r XK » YB » YT ) 

CALL APRi' J TY(Ur-l^»-lt)»15HRPlP2-l- INCHfcSr 12» 690 ) 

i»i— L+4 

50 COliTlNUE 

RFaD 19T (KEC(l) »I=lrl4) 

1X1=NXV(kEC(1) ) 

1 Y i=.ft Y V ( kEC ( M ) ) 

4u pLaQ (9i (KEC(I) r I=lf 14) 
lXP=iMXv(pEC(i) ) 

1 Yc-imY V (nEC (M) ) 

IF 1 1 r 2—1 1 1 ) 34 t 56 1 3b 
54 il <P — i Yi—i Y2 



bO TO 37 
3b IDR=lY2“iYl 
37 IF t IuF-3 ) 36 » 38 » 3b 
3b IF l ( IX2“iXl ) -3 ) 39 * 38 » 3b 
3b CALL LINEV (IX1»IY1»IX2»IY2) 
iXl— 1X2 
1Y1=1Y2 
39 CONTINUE 
K-K+ i 

IF (KNT—K ) 41 » 41 »4u 

11 iF (MOVIE ) 5U d? * 5U 
5U xTRM-0 . 0 

uO 6U U = i*KNT 

KFAL) (ll) (REC(I) * i = D26) »T 
uO hi 1 = 1>26 

IF lXT«M-MBb(KECl 1 ) ) )b2*61»fel 
o2 aTkM— A tiS (REC ( 1 ) ) 
bl CONTINUE 
hU CONTINUE 
kEwINIJ Ijl 

aL— “XTRM— X’ i R m/20 . 0 
XR-X TRH+aTkM/PO.U 
YR-XL 
l T=XK 

l AlL FROimTINUMD 

READ (ll) (REC( I) » i=l»26 ) f I 

UO 1 u 4 KlU6=1#NUM3 

CALL GKlUlV(l»XLfXKr YB r YT » 0 . 0 » 0 . 0 » 0 » 0 » U r 0 » 0 » U ) 

CALL AXlb(XL»XK» PHI t THET A » 0 ) 

CALL DKAm ( KEC ( 1 ) » Phi # THET Arl) 

HEAD (15) (X(D d=ld8) dYd) d=ld8) » (Z( I) d=ld8) >RE 

lb X r RlFY t I 
Rt w I nD lo 

CALL GRAr(X( 1) »Y(1) »Z(1) » PHI » THE! A » RtFX » REF Y ) 
oALL IJNDwRT (KWHICH) 

CAlL RCLuK(170»9u0d00d80d0.0d ) 

104 CONTINUE 
REWIND ll 
REkInD 13 

52 CAlL GRlulV(l#XLrXRf YBf YT»0.0»0.Uf 0»U»0»0r0>0) 

CALL AXlb(XLdRfPHldHETAd) 

IE ( K'n'T—K > 12 r 31 * 51 
51 kEaD (11) (RtC(I) » i=l»26) fT 
CALL DRAjfii(RECd) r PHI » THETA f 0) 

CALL RCLoKd70*9u0» lOOflSOr 10. Or T) 

K=K+i 
GO To 52 

12 CONTINUE 
201 CONTINUE 

IF ( MOV IE ) Io2d02du0 
100 CONTINUE 

UO lu5 KRUU— 1 * 50 
CALL RESET V 
10b CONTINUE 

UO lul KkUU=1»KNT 
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C AlL GKlUl V ( 1 » XL t XK » YB » YT t 0 . 0 # 0 . U » 0 » U r U » U » U t 0 ) 

CAl_L AXIS(XL»XK» PHI t THETA * 1 ) 

HEAD l lb) (XII ) » i=l»lA) » ( Y ( I ) » I=l#18) » IZ(I ) »I=1»18) t HF 

1FX » REF Y » I 

CALL GHAi'i ( X ( 1 ) »Y(1) fZ(l) » PHI » THETA » REFX » REFY ) 

CALL RCL0K( 170»9u0» 100* 180*10.0*T) 

11)1 CONTINUE 

uO 106 Ki-dJU^l t SO 
CALL RLStTV 
lUb CONTINUE 
102 CONTINUE 
kE TURN 
lMu 


alrtFTC ShFI 

SlJbROUT InE CART(XL>XR» YB*YT) 

COMMON VAR * KNT t KFST » L 
DIMENSION VAK(bftUO) 

CALL DXDYV ( 1*XL*XR* UX*N* I»NX* 16 • 0 » IEHR ) 

CALL DXDtV( 2> Yb* YT»OY»M* J»NY» 16.0* IERR ) 

CALL EUG(XL»XR» Yb» YT*DX*OY*N*M*0. 000001) 

IF l L— 1 0 ) 6 * 2 * 6 
2 NX=NX+2 

wY— ny+2 

O i-iX~NX+l 
imY-NY+1 

CALL GRlulV ( i*XL* XR* Yb* YT*UX*DY*N»M*-N*-M*NX»NY> 
CALL PRINTV(— 12* 12HT I ME SECONDS * 400 * 12 ) 

HE TURN 

t-NU 


4IHFTC SKF 2 

SUbROUTlNE EUG(XL*XR*Yb* YT»DX*DY* I*J*C) 
UX=AbS(Dx) 

UY=AbS(DY) 

IF(DX)60*45»60 
6U IF l DY ) 61 * 4S» 61 
61 LONTINUE 
F I — I 
F J-J 

IF I XL ) 28 *8*1 

1 iJ=UX*FI 
S=U 

2 AG-S-XL 

IF l Ab)4*b*d 
0 IF(C-AG) 7*6*6 
7 XL-S-0 
( 5 0 TO b 
b XL-S 
bO TO S 

4 iFlC + AG)27*5*5. 
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27 b-b + u 
oO TU id 
b l)=UX*FI 
b b=U 
9 MG— XK — S 

iF l Ab ) 1 1 » 12 » 1 « 

ID iF ( C— Ab ) i3 » 12 • 12 

11 aH-S 

12 iF l Yi • ) 22 1 lb » 1" 
lb b = b+u 

bO Tu 9 
LH J=UY*Fu 
CK-D 

7b iF lYn/CK-luOU.U) 70*71*71 
71 oK=Ct\+10u0.0*( 
oO To 7b 
7u b=CK 
lb at-’— b-Yb 

iF l Mo ) 19 * lb * 1 7 
17 iF IC-Ab) aH* ] t>> lb 
lb 7P=S-D 
GO To lb 

19 il- lC + Ab)c;l*2u*2U 
2i b— o+u 

uO TO Lb 

20 tH=S 

oO To lb 
2c; U=uY*FJ 
b=-U 

2b MG-Yn-b 

IF l Ao ) 24 1 lb» 25 
24 iF l C+AG ) ab * lb * 15 
2b YH— S 

oO TO 15 
2b b = b~U 

nO TO 23 
2o U=UX*F1 
b=”U 

29 mG— XL” b 

iF ( A b ) 3 0 * 3 3 * 3 ? 

30 IF l C+ Ab ) b 1 * 32 * b2 

31 b=b-U 

bO TO 29 

32 aL-S 

33 IF ( Xh ) 34 * 1<£ * b 

34 b=-0 

3b AG-XK-b 

IF l Ab ) 38 * lb ► 38 
3b 1F(C+Ab)b9rlb»15 
39 S=b-U 

bO To 3b 

3o iF(C-Ab)j7*lb*ib 

37 AR=b+D 

lb IF l Y f ) 4b * 4b* 40 



4u u=uY*Fu 
CK-D 

74 iF IY r/CK-luOU.U)75»7b»76 
7b CK-CK+10u0.0*D 
uO TO 74 
7b b-CK 

41 aG-YT-S 
lFlAG)44»4b»42 

42 1 F i C—A g ) 43 » 4b t 45 
4b b— b+U 

oG TO 4l 
44 YT-S 

*oG TO 45 
4b D=uY*FJ 
b=-D 

47 aG=YT-S 

IF ( Au ) 48 r 4b » bl 
4b IF ( C+AG ) SO » 49 # 49 
5U S=b-U 

oO TO 47 
49 YT-S 

GO TO 45 

51 iF l C— AG ) b2 * 4b • 45 
be: YT-S+D 
4b CONTINUE 
kE 1 UKN 
END 


SluFIC SKFb 


SUBROUTINE FRONT (NUM1) 
uIlviElMSIOu v Ak ( 68U0 ) » H ( o6UO ) 


COMMON VaR 

EGUI VALENCE (VAR (GUI) »P(1) ) 
U = 3Bo»4*P(b) 

V=38b.4*P(22) 


up lu l = l*NU(Vtl 
CALL REStTV 

CALL RITt.2V(150» 169 » 1023 »180»2»36»-lf 38HCABLE CONNECTED SPACE STAT 
1 ION UYNA(»,ICS»NLAST) 

uaIl RiTt2 V ( 300 >250 » 1023» 180 » 1 r 29»-l » 29H WEIGHT OF BODY 1 = 

1LB. rNLASI ) 

CALL R1Tl 2V(435»25U»1023»180#1»29*— lr 29H WEIGHT OF BODY 2 = 

IlB • • NLAST ) 

CALL RITc.2V ( S70 » 178 » 1023* 180 » 1 » 37» -1 » 37HINITI AL SPIN SPEED = 

1 uEG./StC. »NLAST) 


CALL RITt2V 
luY C.G. S= 
CALL RITL2V 
1 PS 1 f NLaST 
CALL R 1Tl 2 V 
CAuL R1TL2V 


(705»88»1023»180» 1 > 48 » -1 r 48HINITI AL DISTANCE BETWEEN BO 
FT. * NLAST) 

(640»232*1023»180»1*31»-1* 31HCABLE ELASTICITY = 

) 

( 690 *700»1023»180»1»1»—1»1H*» NLAST ) 

(975»214f 1023 »180»1»33»-1» 33HT0T AL CABLE AREA = 
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IS (3. IN.MmLAST) 

CALL HEAu(3»3» 18 > 26 » 300 r 592 » 1023»180* l»LX»LYr0»0» 2»U) 

CALL HE.Au ( 3 » 3 * 18 » 2b > 435 » 592 » 1023»l80»l*LX»LY»0»0i2»V) 

CALL HEAu(3»3r 18 , 26 » 570 » 556 * 1023 » 180 , 1 r LX > LY > 2 r 0 , 2 » P ( 5955 ) ) 
CALL HEAD<3*3#18*2o»705f772fl023»180rl»LX»LY»l,0*2»P(5952) ) 
CALL HE AU (3»3»18» 2b * 840 » 574 * 1023 rl80»l*LX# LY »0»0»2»P( 5954) ) 
CALL HEAu (3»3»1B»26»975»556»1023»180»1»LX»LY»4»0?2»P( 5953 ) ) 
10 CONTINUE 
HETUKN 
ENU 


SlbFTC SKI-4 

SUBROUTINE HtAU(LVw»LVH» ispace> irow#mx*my»limit»k*int»lx»ly* irt» 

UVARf IFlKf VAK) 

DIMENSION KEC ( 400 ) 

CALL CAMhAV ( 35 ) 

CALL CHSiZV(LVW»LVH) 

EXTERNAL TABL1V 

CALL RITSTV<ISPACE»IR0W*TA6L1V) 

l\iDS=U 

LX=MX 

LY— M Y 

1FIIVAR.EQ.0)G0 TO 10 
nth=i 

ImTOT AL— 6 

lF<IFIR.I\iE.0)GO TO 52 
READ 1 5 » 5u 1 ) NU 
READ(5»5uO) (KFC(J) »J=1»NU) 

IFIR=1 
52 CONTINUE 
501 FORMAT ( I b ) 

SOU FORMAT ( lb At> ) 

UO 60 L=1»NU 
BCUTXT=REC(L) 

CALL RITE2V(LXrLY»LlMlT»Kf INT»NTOTAL»MTH»BCDTXT»NLAST) 

CALL RlTxYV (LX»LY) 

IF(K.EQ.O) GO TO 61 
1FIK.EQ.90) GO TO 90 
IF (K.EQ.180) GO TO 180 
IF ( K • EQ • 270 ) GO TO 270 

61 IF((LY-4*LVW ) .LT. LIMIT) GO TO 62 
GO TO 60 

90 IF( (LX+4*LVW ) .GT. LIMIT) GO TO 63 
GO TO 60 

180 IF l (LY+4*LVW ) .GT.LIMIT)GO TO 64 
60 TO 60 

270 IFULX-4*LVW ) .LT. LIMIT) GO TO 65 
GO TO 60 

62 LY=MY 
LX=LX-IR0W 
GO TO 60 

63 LX=MX 
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LY=LY-1R0W 
60 TO 60 
64 LY=MY 

LX=LX+IRuW 
60 To 60 
6b LX=MX 

LY=LY+1RuW 
6U CONTINUE 
60 TO 20 

10 IF(VAR.Gt.O.U) GO TO 43 

CALL R ITL2V (LXrLY»LlMlT»K» I NT » 1 » -1 » 1H- 1 NLA ST ) 

CALL RITxYV (LX»LY) 

43 CALL BNBCDV ( VAR»CON»NDS) 

IF ( NUS) 44 * 45 » 46 

44 NDS=IABS(NDS) 

CALL RITt2V(LX»LY»LlMIT»KrINTf2*-lf2H0.» NLAST ) 
IFINUS.GL.IRT) GO TO 47 
UO 46 I = 1 * NDS 
CALL RITxYV (LX»LY ) 

CALL RITL2V (LX »LY» LIMIT ► K » INT » 1 # -1 » 1H0 * NLAST ) 

40 CONTINUE 

NToTAL=IKT-NUS 

IF ( NTOTAl • GT . 6 ) NTOT AL = 6 

CALL RITXYV ( LX r LY ) 

CALL RITL2V ( LX *LY» LIMIT »K * INT » NTOTAL » 1 » CON » NLAST) 
GO TO 20 
47 UO 49 I=i»IRl 

CALL RITxYV (LX.LY) 

CALL RITL2V ( LX > LY » LlMl T » K » INT » 1 » -1 f 1H0 r NLAST ) 

49 CONTINUE 
60 TO 20 

4b CALL RITL2V(LXrLYrLlMlT»K» INT*2»-lr2H0. .NLAST) 
CALL RITxYV (LX»LY) 

IF t IRT • GT . 6 ) 60 TO 73 
CALL RITXYV (LX»LY) 

CALL RITt2V (LX * LY r LIMIT » K » INT» IRT. 1»C0N» NLAST) 

60 TO 20 
73 KUT=6 

CALL R1TXYV(LX»LY) 

CALL RITL2V ( LX »LY» LIMIT » K» INT»KUT»1#C0N» NLAST) 

60 TO 20 

4e> lF(NuS.Lt»6) GO TO 70 
MOKE— NUS—6 

CALL RITl 2V(LX»LY»LIMIT»K» INT » 6 » 1 » CON r NLAST ) 

DO 71 KLU=1 f MORE 
CALL RITXYV(LX»LY) 

CALL RITt2V (LX »LY» LIMIT »,K » INT » 1 » -1 » 1N0 » NLAST ) 

71 CONTINUE 

IF ( IRT . Ec . U ) 60 TO 20 
CALL RITXYV(LX'LY) 

CALL RITt.2V(LX»LY»LlMlT»K»INTrl»-l»lH.r NLAST) 

DO 7b KLU=1»IRT 
CALL RITXYV(LX»LY) 
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CALL R1TE2V ( LX #LY» LIMIT# K» INT# 1 » -1 » 1H0 # NLAST) 

7b CONTINUE 
GO TO 20 

70 CALL R1Te2V(LX#LY»LImIT»K#INT» NDS# 1 # CON# NLAST) 

IF ( IKT • Eg • 0 ) 60 TO 20 

IF ( (NDS+iRT) .LE»6) GO TO 72 

KUI=o-NOb 

NTH=NOS+l 

CALL RITxYV (LX #LY ) 

CALL R ITt2V ( LX # LY » LIMIT #K» INT# KUT# NTH# CON# NLAST) 
GO TO 20 
72 NTH=NDS+1 

CALL RITXYV ( LX ► LY ) 

CALL RITL2V (LX »LY# LIMIT » K » INT # 1 # -1 » 1H. » NLAST ) 
CALL RI (XYV(LX»LY) 

CALL R1Tc2V(LX»LY»LImIT»K» INT# IRT»NTH» CON# NLAST) 
20 CALL RITXYV(LX»LY) 

retukn 

LNU 


SIBFTC SRFb 

SUBROUTINE AXIS(XL»XR#PHI # THETA »LBL) 
bPH— blN(PHI ) 

CPH=COS(PHi) 

STH=S IN (THETA) 

CTh=COS(lHETA) 

CALL POlNTV (U.0#0»0#0) 

XP=XR*CPh*STh 

YP=-XR*SPH 

1X1=NXV(XP) 

IY1=NYV( YP) 

IX2=NXV(-XP) 

IY2=NYV(-YP) 

CALL LINLV(IX1»IY1#IX2»IY2) 

CALL LlNtV(IXl#IYl#IX2#lY2) 

IX=IXH-12 

IY=IY1-12 

CALL RITe2V( IX#IY#1023#180#2»1»-1»1HX»N) 

XP=XR*SPH*STH 

YP=XR*CPh 

IX1=NXV(XP) 

1Y1=NYV( YP) 

IX2=NXV(-XP) 

IY2=NYV(-YP) 

CALL LINEV (IX1»IY1# 1X2# IY2) 

CALL LINEV(IX1#IY1#IX2#IY2) 

IX=IX1+12 

IY=IY1+12 

CALL RITE2V(IX»IY»1023#180»2»1#-1»1HY»N) 
XP=-XR*CTH 



YP=0.0 
1 X1=imX V ( aP ) 

IY1=NYV ( KP) 

1X2=NXV ( u . U ) 
lYk=NY\/(u.O) 

CALL LlNt-V IIX1»IY1»IX2»IY2) 

CALL LIMtV(IXl»lYl»IX2»IY2) 

IX-IXl-lt: 

iY=IYl 

CALL RITL?.\/{iXr IY»1023»180»2»1»-1»1H^»N) 

kEI UkN 

ENU 


3>IbRlC SkF o 

SUBROUTINE PhAw (kEC»PH1 t THLT A t LRL ) 
U I MENS I ON KEC ( 26 ) 
bPfi— s In ( eh i ) 

CPH=COS(pHi ) 

STH=SIN( IHLTA) 

CTM=COS( I HuTm ) 

X=kEC ( i ) 

Y=KEC(2) 

Z=KEC(3) 

ARbl=X*CpH*Sl B+Y*SPH*STH-Z*CTH 
ARb2— —X*SPr1 + lf *CPH 
iXl=NXV(MRbl) 

1Y1—NYV ( mRl2 ) 

UO lu I=i»H 

j=5*l+l 

a=KEC < J) 

Y=KEC( J+i) 

Z=KEC( J+^) 

ARbl-X*CpH*ST H+Y*SPH*STH-Z*CTH 
ARb2— ”X*bPn+ Y*CPh 
1X2 =NXV (nROI) 
iY2=NYV(ARto2) 

CALL LINbV ( 1X1 r IY1» 1X2 » I Y 2 ) 

J.F ( LbL ) lo » 19 » 1 8 
18 CALL IJKU»IX2»IY2) 
iy CONTINUE 
1U CONTINUE 
iXi=iX2 
IYi=IY2 
UO 11 
J=3*L+1 
X— KEC ( U ) 

Y=KEC(U+i) 

L=HEC ( J+c) 

ARbl=X*CbH*STH+Y*SPH*STH-Z*CTH 

ARb2=”X*SPH+Y*CPii 

IX2=NXV(ARbl) 
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iY2=NYV<ARb2) 

CALL LINtVlIxl >IYl»lX2rIY2) 

IF l LbL )2u»21»2U 
2U CALL IUKII» IX2» IY2) 

21 CONTINUE 
11 CONTINUE 

IF(LbL)2b>26r2b 
2 d CONTINUE 
X=KEC<25) 

Y-REC ( 261 

CALL REFr(X» Y »SPti»CPH»STH»CTH) 
2b CONTINUE 
RETURN 
ENU 


3>IbFTC SKF7 

bUbROUTIlME IJK ( I > IX» 1Y) 

EXTERNAL TABL1V 
EXTERNAL TABlbV 
CALL CHSiZV ( 3» 3) 

ILX=lX+9 

!LY=IY-6 

CALL VCHaRV ( 180»1*ILX*ILY» 42 » T ABL1V ) 
bO TO ( 1 » 2 > 3 » 9 » b * fa » 7 » 9 ) » I 

1 IX1=ILX+I2 
IY1=1LY-12 

CALL VCHARV ( 1BU » 1 » 1X1 » I Yl » 25 » TAB15V ) 
fa CALL CHbIZV(2*2) 

IXl=IXl+o 

IY1=IY1+13 

CALL VCHaRV ( 1R0» 1 » 1X1 » I Yl » 1 » TABL1V ) 
bO TO 9 

2 IX1=ILX 
1Y1=ILY+I2 

CALL VCHARV ( 180 » 1 » 1X1 r I Y1 » 33» TAB15V ) 
bO TO 8 

3 IX1=ILX-12 
IY1=ILY 

CALL VCHARV ( 180 » 1 r IXlr I Yl > 34 # TABlbV > 
GO TO 8 
b IX1=ILX+12 
IY1=ILY-I2 

CALL VCHARV ( 180 » 1 1 1X1 » I Yl » 25 » TAB15V ) 
10 CALL CHSiZV(2»2) 

IXl=IXl+o 

IY1=IY1+13 

CALL VCHARV (180»1»IX1» IYl » 2» TABL1V ) 
bO TO 9 
fa 1X1=ILX 
lYl=ILY+2 

CALL VCHaRV ( 180»1» 1X1 » IYl »33»TAB15V) 
bO TO 10 
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7 1X1=ILX-12 
1Y1=ILY 

CAlL VCHmRV ( 180*1* 1X1* I Yl * 34 * TAR 15V ) 
CO TO 10 
G CONTINUE 
KFTUKii 
lNU 


bldFTC SKFtt 

SUdROUTInE RtFR( a*Y *SPH*CPH*STH*CTH) 

x 1=0.0 

Y 1=0 . 0 

SLOPE=Y/X 

bLiNV=X/Y 

UX=X/20.U 

UY=Y/20.U 

X=U.U 

r=u . u 

IF ( AdS l SLOPE) -l.U) 13* 13* 14 
10 UO lb 1=1*20 
Y=bLOPt*x 

xT=X*CPH»STH+Y*SPH*STH 

YT=-X*SPh+Y*CPH 

CALL POINT V ( XT » YT* 0 ) 

X=X+UX 
lb CONTINUE 
\jO TO 16 
14 uO 17 1=1*20 
X=bLiNV*Y 

XT=X*CPH*S I H+Y*SHH*STH 

YT=-X*bPN+Y*CPH 

CALL PoIuTv (XT* YT* 0) 

If = Y +l)Y 
17 CONTINUE 
lb CONTINUE 
KEIUKN 
tNU 


bldFTC SKFG 


10 


bUbR0UTIi\iE 
U I WENS ION 
bPlt=SIMpH 
CPn=COb(HH 
bTh=SlN ( TH 
CTH=COb(TH 
CALL PK0 o( 
UO 1U 1=2* 
CALL PROu ( 
CONTINUE 


X 

) 

) 

) 

) 

X 

y 

x 


GRAN ( X * Y * Z» PH* TH* REFX » REF Y ) 

(18) * Y ( 16 ) *2(18) 

(1)*Y(1)*Z(1) *PH*TH*XCG1* YCGl) 

(I) *Y ( I) *Z( I) *PH*TH*X(I-1) *Y ( 1-1) > 
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CALL OUAu(XCclrYCGl» 1001) 

I GU= 1 QU 1 

CALL PROu ( X ( 10 ) » Y ( 10 ) » 2 ( 10 ) » PH » Th » XCG2 » YCG2 ) 

KlJAL— 1 

CALL QUAu(XCG?» YCG2* IQD2) 

2 U CALL CASl(X( 1) »Y(l)»X(2)»Y(2)rX(3)rY(3)»X(4)»Y<4)»X<5)rY(5)»X(6)r 
lY(b)»X(7)»Y(7)»Xlti)»Y(b) > XCG1 » YCgI » XCG2 » YCG2 » IQD t ICASE ) 

CALL BOX ( X ( 1 ) » Y(l) * ICASE » A » B t C * D r E » F » G » H » 0 t P » 0 * R » S t T r U» V t W , A 
ltif AC) 

IF l KUAL—i )2S»21»22 

21 Xl=X(l) 

Y1=Y(1) 

X2=X(2) 

Y2=Y(2) 

a3=X(3) 

Y 3-Y ( 3 ) 

X4— X ( 4 ) 

Y4=Y(4) 

UO li I=ilrlb 

CALL PROu ( X ( 1 ) »Y(I) »Z(I) r PH * TH » X ( 1-10 ) »Y(I-10) ) 

11 CONTINUE 
KUAL=2 
iQU=IQU2 
ICASe1=IcAsE 
GO TO 20 

22 1CASe2=ICAsE 

CALL KABtL(XCGl» YCGl»XCG2»YCG2rXl»Yl.X2»Y2»X3»Y3fX4»Y4»X(l) »Y(1) » 
1X(2) »Y(2) »X(3) # Y(3) fX(4) #Y(4) »ICASElr ICASE2) 

2b CONTINUE 

CALL REFh(REFX*REFYfSPH>CPH>STH»CTH) 

return 

END 


SlbFTC SKF10 

SUbROUT INE PRO J ( X » Y f Z » PHI t THETA t T X t TY ) 

TX=X*COS ( PmI ) ♦SIN ( 1 HETA ) +Y*SIN(PnI)*SIN( THETA )-Z*COS( THETA) 
rY=-X*SlN(PHl)+**COS(PHI> 

RETURN 

ENU 


SIbFIC SKF11 

SUBROUTINE QUAD l XRST » YRST t IQD) 
If (XRST>200»201»201 

201 IF ( YRST >203*202*202 
200 IF ( YRST >204*204*205 

202 IQD=1 

GO TO 206 

203 IQD=2 
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204 


GO TO 206 
IUD=3 
GO TO 2u6 
20b I0D=4 

20b CoNTINUt 
RtTUHN 
END 


SldFIC SkF 1 ? 

SUBROUTINE CASE(X1>Y1»x2»Y2»X3»Y3»X4»Y4»X5»Y5»X6*Y6»X7»Y7»X8»Y8» 
1XCg1»YCG1» XCG2 » YCG2 1 I QD t IC ASE ) 

Kl=Xl**2+Yl**2 

K2=X2**2+Y2**2 

k3=X3**2+Y3**2 

k4=X4**2+Y4**2 

KM=R1 

ITH=i 

iF (R2-RMJ6U0»603f6U3 
bOU HM=R2 
ITH =2 

b03 IF ( R3— RM )6U1»604»6U4 

601 HM-R o 
ITh=3 

b04 IF ( R 4 ~KM ) 6u2 t 605 t 60 S 

602 RM=R4 
ITH=4 

bOb CONTINUE 

GO To ( 82i * 822*823*824 ) * ITH 

021 CALL S1ML(XCG1*YCG1* XCG2 * YCG2 »X1*Y1*X2*Y2* XL1 * YL1 ) 
RL1=(X1-a 2)**2+( Yl-Y2)**2 
ULI= l XI- aLI ) **2+ ( Y1-YL1 ) **2 
EL1= ( X2-XL1 ) **2+ ( Y2-YL1 ) **2 
IF 1 (UL1.gT.RlD .OR. (EL1.GT.RL1) )GO TO 702 

701 1_INE=1 

GO TO 621 

702 LINE=2 

GO TO 622 

824 CALL SIML ( XCG1 ► YCGl * XCG2 * YCG2 * XI * Y1 * X4 » Y4 » XL2 » YL2 ) 

RL2= ( X1-X4 ) **2+ ( Y1-Y4 ) **2 
UL2= ( X1-AL2) **2+ ( Y1-YL2) **2 
EL2= ( X4-XL2 ) **2+ ( Y4- YL2 ) **2 
IFC (UL2.GT.RL2) .OR. (EL2.GT.RL2) ) GO TO 704 
GO TO 702 
704 LINE=4 

GO TO 624 

822 CALL SI ML ( XCG1 ► YCG1 * XCG2 * YCG2 » X2 » Y2 * X3* Y3* XL3* YL3) 
KL3=(X2-X3)**2+(Y2-Y3)**2 
UL3= ( X2-XL3 ) **2+ ( Y2-YL3 ) **2 
EL3= ( X3-XL3 ) **2+ ( Y3-YL3 ) **2 
IF( (DL3.GT.RL3) .OR. (EL3.GT.RL3) )GO TO 701 

703 LINE=3 

GO TO 623 
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b23 CALL SIML(XCg 1 » YCGl#XCG2# YCG2>X3»Y3*A4»Y4»XL4»YL4) 
KL4= ( X3-X4 ) **2+ ( Y3-Y4 ) **2 
UL4= ( X3-XL4 ) **2+ ( Y3-YL4 ) **2 
lL4= ( X4-XL4 ) **2+ ( Y4-YL4 ) **2 
1F( (UL4.GT.RL4) .OR. (tL4.GT.RL4) >G0 TO 703 
CO TO 704 

621 GO Tu ( b3i t 631 » 633 » 633 ) # IQD 

622 GO TO (641 *642 *642 >641 ) rlQD 

623 GO To ( 651 * 652 * 652 » 651 ) ♦ IQD 

624 GO TO ( 661 * 661 * b63* 663) * IGD 

631 S2b=(Y8-Y2)/(Xb-X2) 

d 2 b=Y 2 -Sc:a*X 2 
d3=Y3-S2b*X3 
IF (B28.L I .63) GO TO 671 
GO TO 67 (j 

b33 Slb=(Y5-Yl)/(Xb-Xl) 
dlb=Yl-Sl5*Xl 
H6=Y6-Slb*X6 
lF(B15.Li.ri6) GO TO 673 
GO TO 672 

641 S4b— ( X6 ~a 4 ) / ( Yb—Y 4 ) 

A4b=X4-S48*Y4 

A3=X3-S4o*Y3 

IF ( A48.L I • a 3) GO TO 67b 
GO TO 67 h 

642 b4b=(Xb“A4)/(Y8“Y4) 

A4b=X4-S48*Y4 

A3=X3-S4o*Y3 

IF ( A48.L I .A3) GO TO 677 
GO TO 67 b 

651 b3b=(Xb-X3)/(Yb-Y3) 

A3b=X3-S55*Y3 
A4=X4-S3b*Y4 
IF (A35.L I . A4) GO TO 678 
GO TO b7G 

b52 S2o = ( X6~X2 ) / ( Y6-Y2 ) 

A26=X2-S26*Y2 

Al=Xl-S2o*Yl 

IF( A26.L I .Al) GO TO 680 

GO TO 681 

663 S46=(Y6-Y4)/(X6“X4) 
d46=Y4-S46*X4 
Bl=Yl-S4b*Xl 
IF (B46.LT .til ) GO TO 684 
GO TO 68b 

661 S37=(Y7-Y3)/(X7-X3) 

B37=Y3-So7*X3 
b2=Y2-S37*X2 
IF (B37.LI .B2) GO TO 682 
GO TO 68,} 

670 ICASE=4 
GO TO 80U 

671 ICASt=12 



60 TO BOU 
672 iCASE=16 
00 TO BOu 
675 1CASE=8 
60 TO BOu 

674 JLCASE=10 
bO TO BOu 

675 1CASE=2 
60 TO BOu 

b76 1CaSE = 6 
60 TO BOu 
677 ICASE=14 
60 TO BOu 
o76 iCASt = l 
60 TO BOu 
t>79 1CA5E=9 
60 TO BOU 
6bU ICASE=i3 
60 To BOu 
681 ICASE=5 
60 TO BOU 
t>82 ICaSE= 3 
bO TO BOu 
b83 1CASE=11 
60 TO BOU 
t>84 1 CaSE=15 
60 TO BOU 
685 1CASE=7 
BOO CONTINUE 
RETURN 
ENU 


SIbFTC SHF 13 

SUBROUTINE SI ML ( XC61 * YCG1 r XCG2 * YCG2 »XP1»YP1» XP2» YP2» X » Y ) 

0XPl2=XPi-XP2 

uYP12=YPi-YP2 

UXCG=XCG1“XCG2 

UYCG=YCGi-YCG2 

X=IDXP12*DYCG*XCG1-DXCG*DYP12*XP1+DXCG*DXP12*YP1-DXCG*DXP12*YCG1)/ 
1 (DXP12*DYCG-DYP12*UXCG) 

Y=(DYCG/UXCG)*X-(0YCG/0XCG)*XCG1+YCG1 

REV URN 

ENU 


SIBFTC SRF14 

SUBROUTINE BOX(X»Y» ICASE>ArB»C»DrE»F»,G»H»0»P»0»R»S»TrU»V»W»AB»AC) 
uImENSIOn X ( 18) » Y ( 18) 

60 TO(ir2’3’4t3’6’3»lr9rlO»lOrl’13r4flO»4 ) * I CASE 
1 IX1=NXV(X(1) ) 
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IY1=NYV< Y <1> > 

UO 20 I=2»b 
IX2=NXVU(I) ) 

IY2=NYV(Y (I) ) 

CALL LINtV(IXl>lYl»lX2,lY2> 

IX1=IX2 

IY1=IY2 

20 CONTINUE 
iXl=NXV(x(l) ) 

1Y1=NYV(Y (1) ) 

LALL LINLV<IX1'IY1'IX2'IY2) 
IX2=NXV(X(4 ) ) 
i Y2=NYV(Y (4) ) 

CALL LINtV<IXl>IYl»IX2»IY2) 
IXi = NXV ( X ( ) ) 

1Y1=NYV<Y (3) ) 

IX2=NXV(x(b) ) 

iy2=nyv(y ( a) ) 

CALL LINLV 1 1X1 1 IY1 1 1X2 r IY2) 
1X1=NXV ( X ( 5) ) 

IY1=|\|YV<Y (b) ) 

CALL LINtV(IXl»lYl»IX2»lY2) 
GO TO 30 

2 IXi=NXV(X(4) ) 

IY1=nYV(Y (4) ) 

DO 21 I = b»8 
1X2=NXV(X(I) ) 

IY2=NYV(Y(1)) 

CALL LINtVUXl»IYl»lX2»lY2) 

lXl=IX2 

IY1=IY2 

21 CONTINUE 
IX2=NXV(X(5) ) 

1Y2=NYV(Y (b) ) 

CALL LINLV(IX1»IY1»IX2#IY2) 
IX1=NXV(X(7) ) 

IY1=NYV(Y(7) ) 

I X2=NXV(X(2) ) 

IY2=NYV( Y (2) ) 

CALL LINEV(IXlrIYl»IX2»lY2> 

IX1=IX2 

1Y1=IY2 

IX2=NXV(X(1) ) 

1Y2=NYV(Y(1) ) 

CALL LINtV<IXl»IYl»IX2»IY2> 

IX1=IX2 

IY1=1Y2 

1X2=NXV(X(4) ) 

IY2=NYV(Y(4) ) 

CALL LINLV(IX1»IY1»IX2»IY2) 
IX2=NXV(X(6) } 

IY2=NYV(Y(6) ) 

CALL LINLV(IX1»IY1» IX2»IY2) 
GO TO 30 

3 IX1=NXV(X(2) ) 
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iYi=wYv ( Y ( 2 } i 
00 22 1=0*6 
1X2=NXV(X(I) ) 

I Y2=NY V ( Y ( 1 ) ) 

CALL LINlV (IX1»IY1#Ix 2»IY2) 
IX1=1X2 
iYl=IY2 
2 l CONTINUE 

1X2=NXV (x(b) ) 

1Y2=nYV( Y ( b) ) 

CALL LlNtV(IXl»IYl* 1x2* IY2) 
lX2=NXV(x(b) ) 

1Y2=NYV( Y(b) > 

(-ALL LINLV (Ixl *IY1*IX2*IY2) 
1 X i=NXV ( X ( 7) ) 
i Y i=l\iXV ( Y ( 7 ) ) 

IX2=NXV(X<2) ) 

IY2=NYV(Y <2> ) 

CALL L INlV l 1X1 » I Y1 » 1X2* IY2 ) 
GO TO bO 
4 IX1=NXV ( A(l) ) 
i Y i=nY v ( r ( i ) ) 
uO 2b I=^*b 
iX2 =NXV (a( i ) ) 

IY2=NYV( Y (1) ) 

CALL LINlV 11X1* IY1* 1X2* IY2) 
IXl=iX2 
IY1=XY2 
2b CONTINUE 

IX2=NXV (X(6) ) 

IY2=NYV( Y(ti) ) 

CALL LINLV<IX1*IY1*IX2*IY2) 

IX1=1X2 

IY1=IY2 

IX 2 =NXV (x<7) ) 

iY2=NYV(Y(7 ) ) 

CALL LINlV (IX 1»IY1*IX2*IY2) 
IX2=NXV ( X ( b ) ) 

I Y2=NYV ( Y ( b) ) 

CALL LINLV(IX1*IY1*IX2*IY2) 
IXi=NXV ( X ( 7 ) ) 

1Y1=NYV(Y(7) ) 

IX2=NXV(X(2) ) 

I Y<;=NYV ( Y ( 2 ) ) 

CALL LINlV (IX1*IY1»IX2»IY2) 
IX1=NXV (X ( 1 ) ) 

I Y i=NXV ( Y ( 1 ) ) 

IX2=NXV(X(4 ) ) 
i Y2=NYV ( Y (4) ) 

CALL LINLV(IX1*IY1*IX2*IY2) 
GO Tu bO 
6 IX1=NXV ( X ( 1 ) ) 

IY1=NYV( Y (1) ) 

00 24 l=e; * 4 
IX2=NXV ( X ( I ) -) 



IY2=nYV( Y ( 1) ) 

CALL LINtV(IXl»IYl»IX2»IY2) 
1X1=IX2 
IY1=IY2 
24 CONTINUE 

lX2=NXV(x(l) ) 

IY2=NYV(Y (1) ) 

CALL LINLV( 1X1 » IY1» IX2» IY2) 

1X1= 1X2 

IYl=iY2 

UO 2b I=b » 8 

1X2=NXV ( X( I ) ) 

IY2=NYV(Y <I) ) 

CALL LINtV (lXl»lYlflX2»lY2) 
IX 1=1X2 
I Yl=IY2 
2b CONTINUE 

IX2=NXV(X<3 ) > 

1 Y2=NYV ( Y ( 5 ) ) 

CALL LINtVUXl » IY1» 1x2* IY2) 
I X1=NXV ( X ( 2 ) ) 

I Y1=imYV ( Y ( 2 ) ) 

I X2=NXV ( X ( 7 ) ) 

IY2=NYV(Y (7) ) 

CALL LINLV (IX1»IY1»IX2»IY2) 
GO TO 30 
9 1 X1=NXV ( X ( 1 ) ) 

1Y1=NYV(Y (1) ) 

00 26 1=2 » 7 
IX2=NXV ( X ( I ) ) 

IY2=NYV<Y(i) ) 

CALL LINtVl IX1» IY1» IX2»lY2) 
IX1=IX2 
IYI=IY2 
2b CONTINUE 

I X2=NXV ( X ( 2 ) ) 

IY2=NXV( Y (2) ) 

CALL LINtV( IX1» IY1» IX2» IY2) 
IX1=NXV ( x (6) ) 
lYl=NY\/ ( Y (fa) ) 

IX2=NXV ( X ( 1 ) ) 

I Y2=NYV ( Y ( 1 ) ) 

CALL LINtV(IXl»lYl»IX2flY2) 
IX1=NXV ( X ( 4 ) ) 

IY1=NYV(Y(4) ) 

CALL LINtVilXlr IYl»IX2rlY2> 
bO TO 30 
10 IXl=NXV ( X ( 3) ) 

IY1=NYV(Y(3) ) 

UO 27 1=4 » 8 
IX2=NXV (X( I ) ) 

IY2=NYV(Y(I) ) 

CALL LINEVdXldYlf IX2»IY2) 
IX1=IX2 



m=iY^ 

27 CONTINUE 

IX2=NXV ( A< 3) ) 

IY2=NYV(Y(3) ) 

C Ai_L LINtV< Ixir lYldX2»lY2) 
i X2— NX V ( X ( b) ) 
i Y cL — H'i'v ( Y (5) ) 

C Ai_L LlNfc.VtIXl»lYl»IX2»lY2) 
IXi=NXv(A(l) ) 
iYi=wYV( t (i) ) 

1X2=NXV ( a (4) ) 

IY2=NYV( l (4) ) 

CALL LINtV(IXl»lYl»IX2»lY2) 
1 Xl = NXV ( A ( to ) ) 
iY^=NYV(Y(t)) ) 

CAlL LINcVdXl* IYI* Ix2» IY2) 

00 TO 30 

lb IX1 =NXV (a( 1) ) 

1Y1=nYV<Y (1) ) 
uo 2a i=^o 

iX2=NXV(A(i) ) 

1 Y 2 = NYV (til)) 

CAlL L1NlV( I xl » IY1» IX2» IY2) 
iXl=iX2 
IY1=IY2 
2b CONTINUE 

IX2=NXV ( A (a) ) 

iy^=nyv( y ( a) ) 

CAlL L1NlV 1IX1»IY1»IX2»IY2) 
lXi=NXV(x(b) ) 

IY1=NYV( Y (b) ) 

CALL LlNtV<IXl » IY1» Ix2» IY2) 
UO 2S< I=o»a 
IX2=NXV(A(I) > 

1Y2=NYV (Y ( I) ) 

CALL LINtV(Ixl> IYI* IX2dY2> 
1X1=1X2 
IYI=IY2 
2V CONTINUE 

iXi=NXV(x(i) ) 
lYi=NYV(Y (1) ) 
iX<£ = NXV ( A ( b ) ) 

IY2=NYV( Y (b) ) 

CALL UNLV ( lx 1 » IY1» I X2 » I Y2 ) 
1X1— NXV ( A ( 2 ) ) 

IY1=NYV( Y (2) ) 

IX2=NXV(x(7) ) 

IY2=NYV( Y < 7 ) ) 

30 CONTINUE 
KE 1 Ukn 

fc NU 


SIbFIC SKFI5 
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SUBROUTINE KAHEL ( XCG1 » YCG1 » XCG2» YCG2» XI »YlrX2»Y2»X3»Y3»X4»Y4»X5»Y5 
1» Xfo»Y6f XY# Y7»X8» Y8» ICASEl * ICASE2 ) 
iCASE=lCASEl 
KUAL=1 
XCG=XCG1 
YCG=YCG1 

50 GO TO ( 1 * 2» 3» 1 f 5» 1 • 3 * 1 * 1 1 2 » 3» 1 1 5> 1 * 3» 1 ) r ICASE 

1 IXC1=NXV(XCG) 

IYC1=NYV(YCG) 

GO TO 20 

2 CALL SIML(XCGl»YCGl#xCG2»YCG2fXl»Yl»X4fY4»X»Y) 

30 IXC1=NXV(X) 

IYC1=NYV(Y) 

GO TO 20 

3 CALL SIML(XCGlfYCGl»XCG2» YCG2fX3» Y3»X4»Y4»X»Y) 

GO TO 30 

b CALL SIML(XCGlfYCGlfXCG2»YCG2»X2»Y2fX3fY3»X»Y) 

GO TO 30 
20 CONTINUE 

IF IKUAL— i ) 42 r 41 » 42 

41 lXC2=IXCi 
1YC2=IYC1 
XCG=XCG2 
TCG=YCG2 
KUAL=2 
1CASL=1 CaSL 2 
Xl=Xb 
X2=Xfa 
X3=X7 
X4=Xb 
Yl=Yb 

Y2=Y b 
Y3=Y7 
Y4=Yb 
GO TO 50 

42 CALL LINLV( IXClUYClt IXC2» IYC2) 

RETURN 

END 


SIBFTC SRF16 

SUBROUTINE UNDWRT UWHICH) 
DIMENSION VAR(b800) »P( 6800 ) 
COMMON VAR 

EQUI VALENCE ( VAR (faOl) »P(1I) 

EXTERNAL TAB15V 

EXTERNAL TAbLIV 

Vl=Pt595ti)/12.0 

V2=P(5957)/12.0 

V3=P(595G)/12.0 

CALL CHSIZV (3»3) 

CALL RITSTVt 18» 2b » TABLIV ) 
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CALL 

RiTE2V(623» 97 »lt)23»18U»l» 45 » -1 » 45hT ORQUE 

ABOUT 

r - 

SI 

lN< 

T) FT.-Lb. »NLAST) 




CALL 

RITL2V(873»97f 1 023 » 180 » 1 » 45» -1 » 45HT0RGUE 

ABOUT 

— 

SI 

1N( 

1) FT.-LB. »NLAST) 




CALL 

RITt2V(923r97»l023fl80» 1 » 45 » “1 » 45HT0RQUE 

ABOUT 

— 

SI 


IN ( T) FT. -Lb. * NLAST) 

CALL VCHaRV ( 180»l»b32»325»25 
cAlL VCHmR V ( 18U » 1 * b82 » 325 » 33 
CALL VCH«RV( iflU*lf932»325»34 
CALL CHSiZV ( 2 » 2 ) 

CAlL VCHmRV( 1B0»1»B38» 338 ► 1 » 
CALL VCHMRV(180»l»888f338»lr 
CALL VCHARV(180rl»938»338»l» 
CALL H£AU(3r3»18f2o»823»421r 
CALL HEAU(3»3f 18 » 2b > 873 » 421 » 
CALL MEAu ( 3 > 3 » 18 » 2b * 923 » 421 » 
CALL HE Ac (3»3»18»2b»823»bl9» 
CAlL HEAu ( 3 * 3» 18 t 2b t 873 * 619» 
CALL HEAu(3»3»18»2b»923»619» 
RETURN 
ENU 


» TAB15V ) 
» TAB15V ) 
» TAB15V ) 


TABL1V) 

TABL1V) 

TABL1V) 

1023»180»1»LX»LY#0»0»2#V1) 

1023» 180»1#LX»LY»0»0>2»V2) 

1023 »180fl*LX» LY »0»0»2rV3) 

1023 » 180 » 1 » LX »LY»2»0r2»P( 5999) ) 
1023 » 180 rlrLX»LY*2r0»2rP( 5999) ) 
1023 » 180 » 1 » LX »LY » 2» 0 » 2» P( 5999) ) 


SloFTC SKF 17 

SUBROUTINE RCLOK ( IX » 1Y » IK#K» FACTOR » T) 

EXIERNAL TABL1V 

R=1R 

PI=3. 1415927 
C=2.U*PI*R 
lX=IX 
LY=IY 

IF t IK • LT • ISO ) OO TO 13 
CALL CHSiZV ( 4 » 3 ) 
k1=IR-22 
GO TO 14 

13 CALL CHSiZV (2*2) 

Kl=lK“16 

14 CALL RITbTV(23»26»T ABL1V) 
rt2=0.9*Rl 
OTh=2.0*Hl/10.0 

|H=0.0 

CX-IX 

cY=IY 

UO 11 1=1»10 

TH=TH+uTh 

lX=CX-R*COS(TH) 

LY=CY+R*SIN(TH) 

bO TO( 1»2»3*4»5»6»7»8»9> 10 ) * I 

1 CALL RITl2V(LX»LY»1023»K * 1 » 1 » -1 » 1H1 * NLAST ) 
GO TO 11 

2 CALL RITe2V(LX»LY»1023»K t 1 ? 1 t -1 t 1H2» NLAST) 
O0 TO 11 
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o CAUL RITL2V(LX»LY*1023»K 1 1 > 1 » -1 r lHbt NLAST ) 
oO TO 11 

4 CALL RlTt.2V(LX»LY» 1023 » K » 1 » 1 * -1 » 1H4» NLAST) 

GO TO 11 

b CALL RITl2\/(LX'LY>1023*K # 1 » 1 ► -1 t IHb » NLAST ) 

GO To 11 

b CALL R ITt.2V ( LX » LY » 1023 * K >l»l»-l*lHb» NLAST ) 

GO TO 11 

7 CALL RITl2V ( LX * LY » 1023 » K »1»1»— 1» 1H7» NLAST ) 

GO TO 11 

8 call R1Tl2V(lX»LY»1023»K »1»1»-1*1H8»NLAST) 

bO TO 11 

G CALL RITt2V(LX»LY*i023*K # 1 » 1 » “1 » 1H9 » NLaST ) 

GO TO 11 

1U CALL RITl2V(LX»LY» 1023»K 1 1 » 1 1 -1 * 1H0 t NLAST ) 

11 CONTINUE 
uTH=2.U*Pl/lOO.O 
fH=0.0 

CALL SLT|ViIV(0»0r0»U) 

CALL GRlulV (2 »0.U» 1023.0 » 0.0 » 1023* 0 » U . U » 0 . 0 » 0 r 0»0»0f 0*0) 

uO 1 2 I=i'10U 

Al=CX-Rl*CoS(TH> 

Yl— CY+Rl*SIN(TH) 

CALL POlNTV(Xl.Yl.u) 
f H-T H+UT h 

12 CONTINUE 
uTh=2.0*pI/1(J.U 
IH-0.0 

uO lb l—l » 1 0 
Xl— Cx—Rl*COS ( TH ) 

Yl=CY+Rl*SINlTH) 

CALL POINT V ( XI »Yl#b) 
iH=TH+UTh 

lb CONTINUE 

ANG = 2 • U*R I*T 

IF (K.NE. 180 )mMG=AIMG+1. 570795 
LX=CX-R1*C0S( ANG) 

LY=CY+Rl*SiN(ANG) 

CALL LINtV(IX*IY»LX*LY> 

CALL LINtV (IX»1Y*LX»LY) 

ARG=ANG/EACT0R 
LX=CX-R2*COS(AkG) 
lY=CY+R2*S1N(AKG) 
call LINLVt IX» lYf LX»LY) 
call LINtVl IX» IY»LX»LY) 
aRG=ANG/10u.U 
LX=Cx-.8*R2*C0S( ARG) 
lY=C Y+ . 8*R2*SIN ( ARG ) 

CALL LlNtV(IX»lYrLX»LY) 

CALL LINlVIIx. 1Y»LX»LY) 

rE IUkN. 

lnu 


71 





PRECEDING PAGE BLANK NOT FILMED. 

APPENDIX B 

COMPUTER PROGRAM DATA 


This appendix presents a description of program input, program output, and a 
sample problem. 


General Input 

This section presents a description of program input, including correct format 
and proper ordering. Format statements are given at the beginning of the particular 
group of data to which the statements apply. 

The following fixed-point data should be punched on a single data card in the 
order given. This card should be the first card in the data deck. 


FORMAT (8 15) 


A code number indicating the desired forcing function option (The integer 1 is 
used to force body 1, the integer 2 is used to force body 2, and the integer 3 is 
used to indicate no forcing function. ) 

Number of floating-point values to be input 

6 (The number of auxiliary differential equations, such as control equations, to 
be integrated in the Runge-Kutta (R-K) subroutine is given in this datum location. 
This number is presently 6. ) 

N (This number must be less than or equal to 19. ) 

Number of differential equations to be integrated in the R-K subroutine (This 
number is equal to 3N + 25 + the number of auxiliary differential equations to be 
integrated in the R-K subroutine. This number is presently equal to 3N + 31. ) 

Number of integration steps desired between output printout intervals 

An S-C 4020 output option code number (The integer 1 calls for nine graphs to 
be output. The integer 0 omits the graphs. ) 

An S-C 4020 output option code number (The integer 1 calls for output in the 
form of motion pictures. The integer 0 omits this output. ) 

Each line of the following data should be punched on a single data card. Data 
having a value of zero may be ignored. The order of these cards in the data deck is : 
unimportant. 
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FORMAT (15, El 5. 7) 


1 Integration step size, sec 

2 Program termination time, sec 


3 

V,! 

4 

V,i 

5 

X k',l 

6 

Mi 

10 

h ' T i 

11 

V H 

12 

V~ k i 

13 

V T i 

14 

VH 

15 

V4 

16 

V T i 

17 

VH 

18 

V*i 

19 

X i«,2 

20 

V, 2 

21 

X k»,2 

22 

m 2 


\ 


Direction cosines 
for body 1 


) 
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Direction cosines 
for body 2 


26 

*2 ' X 2 

27 

V T 2 

28 

i 2 ' k 2 

29 

V T 2 

30 

V T 2 

31 

V4 

32 

V T 2 

33 

4 ' ^2 

34 

V 4 

110 

u'l 

111 

v i 

112 

W’j’ 

113 

0 i 

114 

\ 

115 

. 

119 

U” 

2 

120 

V 2 

121 

w 2 

122 

9 2 

123 
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124 

*2 

970 


971 

Cl m 

y, i 

972 

°z,l 

980 

°x,2 

981 

Cl n 

y,2 

982 

°z,2 

140 

X 1 

141 


142 

Z 1 

143 

*2 

144 

? 2 

145 

1 

1 

140 + 3N 

h 

1 

*N + 1 

141 + 3N 

Vl 

142 + 3N 

? N+1 

222 

i 

I 

?P, 1, 2 
1 

222 + N 

1 

X P, 1, N+l 

242 

i 

i 

Tp, 1, 2 

241 + N 

Y P, 1, N+l 
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I » 


261 + N 

Z P, 1, N+l 

282 

l 

i 

?P, 2, 2 
1 

281 + N 

i 

A P, 2, N+l 

302 

i 

i 

Yp, 2, 2 
1 

301 + N 

Y P, 2, N+l 

322 

i 

i 

fp, 2, 2 

321 + N 

1 

Z P,2,N+1 

1202 

i 

I 

? L 2 

| 

1 

1201 + N 

CL N+1 

1222 

i 

i 

CM 

w 

u — 

i 

1221 + N 

tK N+l 

1242 

I 

l 

?°2 

1 

1241 + N 

™N+1 

5945 

<ti P> n 

59 46 

AF 

z, n 

5947 

AF 

y»n 

5948 

AF 

x, n 
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5952 Initial distance between c.g. ^ and c.g.g, ft* 

5953 Total cable area, sq in. * 

5954 Cable elasticity, lb/in^ 

5955 Initial spin speed, deg/sec* 

5956 AG 

z, n 

5957 AG 

y,n 

5958 AG 

x, n 

5966 One-half of length of box representing body 2 in motion picture 

output, in. * 

5967 Length of arbitrary body-fixed axes i^, j n> and k n , and 

one-half of length of box representing body 1 in motion pic- 



ture output, in. * 


5971 

Angle of rotation of 
t 

picture, rad' 

2’ axis out of projected plane of motion 

5972 

Angle of rotation of 

X’ and Y' axes about the Z' axis, rad^ 

5992 

5 i 


5993 



5994 



5995 

S 2 


5996 

% 



$ ~ 

This input is required for motion picture output only. 

^If these angles are both equal to zero, the Y'Z' plane will be parallel to each 
frame of the motion picture. 
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5997 


5999 




Data for the first of a series of runs to be made must include a card for all non- 
zero floating-point parameters as well as the card of fixed-point data. Data for suc- 
ceeding runs may omit any nonzero floating-point parameters which remain unchanged 
from the respective preceding rim but must include the card of fixed -point data. 


General Output 

Table B-I relates the program output symbology to the symbols section of the 
paper. Time histories of the first 36 dependent variables (THETA1 to CABLE) (ta- 
ble B-I) will be output by the system printer for every run made. In addition to this 
fixed output format, two optional forms of output are available through the use of the 
S-C 4020 high-speed microfilm recorder. This option is controlled by two fixed-point 
numbers. The first of these options consists of nine graphs representing the time his- 
tories of the structural Euler angles, the pseudorigid body Euler angles, and the 
pseudorigid body length. The pseudorigid body Euler angles are subject to the following 
restrictions. 


2 < RB < 2 
0£ *'eb < 2 » 


The restriction on concerns only the output values, that is, when the value of 


^RB reaclies 


2 it , a value of zero will be plotted, causing a discontinuity in the graph. 


If the S-C 4020 output section of the program is removed, as described previously, 
there will be no output record of either the structural Euler angles or the pseudorigid 
body Euler angles. This output may be retained by modifying the SUBROUTINE 
OUTADD, as indicated by the applicable comment cards, and by making the following 
changes in the headings printed by the system printer. 
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Original headings 

Modified headings 

THETA1 

THETASTl 

PHIl 

PHISTl 

THETA2 

THETAST2 

Pffl2 

PfflST2 

XBR(l) 

PSIST1 

YBR(l) 

PSIST2 

ZBR(l) 

THETRB 

RP1P2(2) 

PSIRB 


A second output option consists of an 8 -millimeter motion picture of vehicle mo- 
tion during the run. The motion picture consists of two distinct parts. The first part 
shows pseudorigid body and arbitrary body axes motion during the run. The second 
part represents complete vehicle motion. The vehicle is simulated by two rectangular 
parallelepipeds connected by a single cable tied to the geometrical center of the two 
opposing faces. The length of each body is controlled by data input. Body width and 
height are set by the program to 0. 6 of the body length. Both parts of the movie show a 

dotted reference line in the X’Y' plane. This line is the projection of that part of the 
pseudorigid body between c.g. com p and c.g. 2 and serves to give the viewer a qual- 
itative idea of the magnitudes of and 0 RB - Both parts of the movie also show a 

three-handed clock in the upper right-hand corner. One revolution of the largest hand 
represents 1.0 second of vehicle motion. A revolution of the middle hand represents 
10.0 seconds, and a revolution of the smallest hand represents 100. 0 seconds. The 
S-C 4020 generates one frame for every integration step; therefore, an integration step 
size of 0.04167 second will result in a real-time movie. Motion picture output should 
be used with discretion since it greatly increases computer run time. 


Sample Problem 

The vehicle configuration used for the sample problem was taken from refer- 
ence 1 and is shown in figure B-l. Body 1 is manned, and body 2 is an empty booster 
casing. The launch weight of body 1 was approximately 19 000 pounds, which is within 
the Saturn C-l payload capability. The spinning configuration of body 1 includes a 
small unmanned resupply vehicle and two Gemini capsules for emergency escape. 

Four sample runs were made. A copy of the input data cards for the runs is 
given in figure B-2. The data for each run are headed by the card containing eight 
fixed-point values for that run. The first three runs demonstrate how the program can 
be used to obtain the dominant response characteristics of the vehicle for the arbitrary 
axes of body 1. The forcing function used for run one was 


G y x = 1400 sin (1.5t) (Bl) 
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The pertinent output from this run is shown in figure B-3. This curve was output by 
the S-C 4020 and was used to obtain the point indicated on figure B-4. The forcing 
function used for run two was 


G 1 = 1400 sin (0. It) (B2) 


The pertinent output from this run is shown in figure B-5. This curve was used to ob- 
tain the indicated point on figure B-6. The forcing function used for run three was 

G 1 = 30 sin (0. i5t) (B3) 


The pertinent output from this run is shown in figure B-7. This curve was used to ob- 
tain the indicated point on figure B-8. A complete linear analysis of the uncontrolled 
dynamic response of the vehicle is presented in reference 1. The additional response 
characteristics given in reference 1 can also be determined by the program presented 
here. The final sample run was made to demonstrate all available output formats. The 
forcing functions used for this run were 


G 1 = 4000 sin (0. 25t) 
x, i 


G y x = 1 500 000 sin (0.25t) 


and 


G 1 = 1 500 000 sin (0.25t) 


(B4) 

(B5) 


(B6) 


The first three sheets of printed output are shown in figure B-9. The nine graphs out- 
put by the S-C 4020 are shown in figure B-10. Typical S-C 4020 motion picture output 
is given in figure B-ll. The first two frames shown in figure B-ll appear prior to the 
motion pictures of the actual run and provide run identification information. The last 
two frames in figure B-ll were taken from the two types of movies produced during the 
run. 


The sample runs presented illustrate only one of many possible applications of 
the program . The nonlinear approach and the generality of the subroutine structure 
make the program highly adaptable to any type of motion study desired. 
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TABLE B-I. - OUTPUT SYMBOLS 


Program output 
symbols 

Variables 

represented 

Program output 
symbols 

Variables 

represented 

TIME 

t 

YCG 

Y’ 

c.g. 

THETA1 

e i 

ZCG 

Z f 

c.g. 

PSIRBD 

^RB 

XBR2CG 

I 

c.g. 

PHI1 

*1 

YBR2CG 

Y 

c.g. 

THETA2 


ZBR2CG 

Z 

c.g. 

THETRBD 

^RB 

OMEGAX1 

S2 . 
x, 1 

PHI2 

^2 

OMEGAY1 

0 , 
y> i 

THETAB 

PSIB 

PHIB 

e 

W 

? 

OMEGAZ1 

OMEGAX2 

0 - 
z, 1 

n o 
x, 2 

XBR( 1) 

x i 

OMEGAY2 

0 o 
y,2 

YBR(1) 


OMEGAZ2 

9 

Z, a 

ZBR(1) 


RP1P2(2) 

P P 
1, 2*2, 2 

GAMMA 

y 

RP1P2(3) 

P P 
1, 3*2, 3 

ALPHA 

a 

RP1P2(1) 

P P 

i, r 2 , i 

IAVB2 

A 2 

RP1P2(4) 

p p 

l, 4*2, 4 

THETABD 

PSIBD 

PHIBD 

t 

V 

i 

FCABLEMAX 

CABLE 

F 

c, max 

c r 

f, max 

XCG 

X' 

c.g. 

THETAST1 

6 i 

s, 1 
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TABLE B-I. - OUTPUT SYMBOLS - Concluded 


Program output 
symbols 

Variables 

represented 

Program output 
symbols 

Variables 

represented 

PfflSTl 

V 

THETRB 

0 RB 

PSISTl 

\i 

PSIRB 

^RB 

THETAST2 

6 o 

T 

t 


s, 2 



PH3ST2 

V 

X 

X’ 



Y 

Y’ 

PSIST2 

*s,2 

Z 

I' 
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f igure B-l. - Vehicle configuration for sample problem. 









1 


1 


0 


1 109 

6 8 

1 

5-2 

2 

24+0 

3 

1236+3 

4 

1086+3 

5 

2076+3 

6 

10167-2 

10 

1+0 

14 

1+0 

18 

1+0 

19 

36+4 

20 

876+3 

21 

876+3 

22 

46417-3 

26 

1+0 

30 

1+0 

34 

1+0 

111 

- 207022-3 

120 

453454-3 

972 

2292-2 

982 

2292-2 

140 

1651192-3 

143 

1199992-3 

144 

1536-1 

146 

1199992-3 

147 

- 1536-1 

149 

1199992-3 

150 

- 1536-1 

152 

1199992-3 

153 

1536-1 

155 

1199992-3 

157 

- 1536-1 

158 

1199992-3 

160 

1536-1 

161 

1199992-3 

163 

1536-1 

164 

1199992-3 

166 

- 1536-1 

222 

1884-1 

223 

1884-1 

224 

1884-1 

225 

1884-1 

226 

1884-1 

227 

1884-1 

228 

1884-1 

229 

1884-1 

242 

- 768-1 

243 

768-1 

244 

768-1 

245 

- 768-1 

246 

- 768-1 

247 

- 768-1 

248 

768-1 


Figure B-2. - Input data for sample runs. 
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249 

768-1 

262 

768-1 

263 

768-1 

264 

- 768-1 

265 

- 768-1 

266 

768-1 

267 

- 768-1 

268 

- 768-1 

269 

768-1 

282 

- 2628-1 

283 

- 2628-1 

284 

- 2628-1 

265 

- 2628-1 

286 

- 2628-1 

287 

- 2628-1 

288 

- 2628-1 

289 

- 2628-1 

302 

768-1 

303 

- 768-1 

304 

- 768-1 

305 

768-1 

306 

- 768-1 

307 

- 768-1 

308 

768-1 

309 

768-1 

322 

768-1 

323 

768-1 

324 

- 768-1 

325 

- 768-1 

326 

- 768-1 

327 

768-1 

328 

768-1 

329 

- 768-1 

1202 

120515-2 

1203 

120515-2 

1204 

120515-2 

1205 

120515-2 

1206 

120515-2 

1-207 

120515-2 

1208 

120515-2 

1209 

120515-2 

1222 

228706-3 

1223 

228706-3 

1224 

228706-3 

1225 

228706-3 

1226 

228706-3 

1227 

228706-3 

1228 

228706-3 

1229 

228706-3 

5957 

14+2 

5992 

2+2 

5993 

2+3 

5994 

2+2 


Figure R^2. - Continued. 


86 



4 


55 


1 


5995 

5996 

5997 
5999 

1 

2 

5956 

5957 
5999 

1 

2 

5956 

5958 
5999 

1 

2 

5952 

5953 

5954 

5955 

5956 

5957 

5958 

5966 

5967 

5971 

5972 
5999 


1851192-3 
2+3 
2+2 
15-1 
6 8 

7+1 
14+2 
0+0 
1-1 

6 8 

1+2 
0 + 0 
3+1 


15-2 

13 6 8 

25+0 
1376-1 
1512-4 
14583333+0 
2292-2 
15+5 
15+5 
4+3 
2628-1 
1884-1 
5235-4 
5235-4 
25-2 


55 


55 


1 


1 


Figure B- 2. - Concluded 
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Figure B-3. - Pertinent S-C 4020 output for sample run one. 
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Figure B-5.- Pertinent S-C 4020 output for sample run two. 
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■■■■■■IBWmHM BMmMMMMlWWllimMIlllilil 

■■■■■HmamiiMm MBMMamiiBMiimaigflgilggggg 

■jMiaBMiaa BMiiaiBaiiaMaitaiHn 
« Miiajawaiiaien iMiBBi ma 

^■^^il■■l^^■^■■■^■■^^^^■^^^■g■»^■■gggggg|ggggg^lligggggg 


b ia aa fflBHMBSisap 

ImaniuiiwiiiMaiiBiaiilMl] 


ISSESSSSBSSSSiSBSSSSSSSSl 


■■■■■■■■■■iiumuuglggggggggii! 


BaiaBBeaas aaaaaaMHHaiiWiga 

aaaaaaaaaaanaaaaamaaaaaaaaaaaaaaaaaf gag g ggggggggggggg 

iaaaBBii B:atamB gai;aaMiBB!aaBa 

a aaaiBaMBiai BaaBBiiBBaaBaa 

■■■■■■■■■■■■lMf,BBMW |gBBgBBBBlBBflB!IBMBgiWBBBiBBB 

■■■■■■■■■BBmiaHBlBmBBimggBgglg fjggggiggg|gigggBB|l 

■■■■■MMlBlllBBWBMBBIBBBMBBBHg KggggMgggggBBBglMB 

aaaaaaaaaaaawaaaaaaaa aaaaaaaaaamaaaaaaaaaaajagggggg 

■HMMBHBM^BBMBlIBBBBBBggB ggggggggSgliiBBBBBBBIBBB 

■■■■■■■■■■■■■■■■BHBlIBBIBIIBlIBBBBIIIBIIBiBBIllBlIIIII 


!SShSSSSSSSSSSSES5SSSS5BSS|HHHHHH| 
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Figure B-7. - Pertinent S-C 4020 output for sample run three. 


92 







FUKFUN CPTIUN 1 PAGE NO 



AJ-^-inA-X^O-CAXOA- 

>OCCN^rt\-OONO'XO 

in 

X 

A- 

A- 

X 

4- 

Aj 

-4 

X 

X 

, ^1 

-4 

AJ 

AJ 


— * cr — * x -* o o a~ r -4 — <>#- 

C-'ONO'inOJtT'-tN 

-r 

X 

—4 

X 

4- 

4- 

X 

X 

4“ 

X 


X 

0 

f- 


OO^N>rmNfn^rCif\ 

IT (f -t -4 (MD N 0 in N 

r- 

•4- 

0 

X 

X 

X 

0 

X 

A- 

X 

AJ 

X 

X 

X 

a z 

OOOOO^iMWiPNtMM 

mco<M'OCT'fAN--4inco— • 


A- 

0 

AJ 

4- 

X 

O' 

-* 

X 

X 

O' 

AJ 

X 

O' 

CD *- 
aj 

—4 


in 

X 

X 

X 

X 

X 

X 

A- 

A- 

A- 

A* 

X 

X 

X 


Q* sf H ff 1 a i<l ^ N >0 |U 

Off'fllNM'OO'vOfOO'fCl 


0 

Al 

X 

A- 

A- 

AJ 

X 

O' 

X 

0 


0 

X 


O in KO « O (? O' «1 >f 

—iGfrXA*«d'XOin4'X--4 

•4- 

-4 

—4 

X 

A- 

Aj 

X 

4- 

0 

X 

A- 

A- 


Aj 

— 

oo— '-j-A-xoa'-uxx 

n^oNN^fncoi/'Wfij 

—4 

—4 

—1 

—4 

—4 

AJ 

AJ 

X 

4- 

4- 

X 

X 

A- 

O' 

o' z 


*•••••• •••« 

• 

• 

a 

e 

• 











000000-4Alf\JvJ-XA- 

T— 'st'CTAlXCAXA--^ 

in 

O' 

X 

A- 


X 

O' 

X 

A- 


X 

O' 

X 

A- 

>- 

i l t 111 ll l 11 

1 M-4-4^4fSJMCMfflrt«t 

1 1 11 II III 1 

1 

1 

X 

1 

X 

1 

X 

1 

X 

1 

X 

1 

1 

A- 

1 

X 

1 

X 

1 

X 

1 

O' 

1 

T 


(M N <f <5 «0 O 0 3 ^ W in N 

co in -o 0 ^ai sfoxoo 

AJ 

>t 

•4* 

X 

4- 

X 

O' 

0 

A- 

X 


O' 

-4 

X 


O'O'O'O'O'OOO'O'* £ ro 

C'-j-Ncoxir,4 l nOMr' 

X 

X 

0 

X 

AJ 


X 

X 

X 

0 


X 

X 


— 


OOCNcoN-'Ocnm-40'N- 

in 

AJ 

0 

X 

X 

O' 

4- 

O' 

4- 

O' 

X 

A- 

-4 

X 


—4 -4 —4 1—4 —4 — 4 r—i r—A —4 — 4 —4 —4 

-t^OOOOOOOO'O' 

CN 

O' 

O' 

X 

X 

A- 

A- 

X 

X 

X 

X 

.* 

4- 

X 

<* 

xxxxxxxxxxxx 

in in in in in m m m in sf - ^ 

•4- 

>4- 

4- 

4- 

4- 

4- 

4^ 

4- 

4- 

4^ 

4* 

4" 

4- 

4- 

cc z 

\0 .O X X X \0 X X X X X 

'C'O'Oo-C'OO'O'O-b'C 


0 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

0 

X >— 

—4 —4 —4 —4 —4 ■— 4 -4 —4 —4 —4 —4 —4 

—4 —4 1— 4 —4 -4 —4 —4 —4 H p4 H 

—4 

— • 



—4 


-4 


*-* 

-4 

—4 

H 

—4 

— 1 


OO-'N^nNoh ■* AJ 

^-iffiM'n'O'tO'M'fin 

X 

A- 

4- 

4- 

X 

X 


O' 

X 

X 

X 

X 

X 

-4 


OOOaOO-*A|X,fr x 

c04vjcorgNfno'Nin-f 

>4- 

»4- 

X 

X 

A- 

X 

O' 

X 

A- 

X 


X 

A- 

X 


ooooooooooo 

0-41-4— JCNJAi'-nCI-IsJ-U-'O 

A- 

X 

O' 

0 


AJ 

X 

4- 

X 

X 

A- 

r- 

A* 



ooooooooooo 

ooooooooooo 

0 

0 

0 

-* 

-4 


- 1 


—4 

—4 

—4 




-c 

- -3 
X uj 

66066690606900066600000 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

f\jNJ3COCDM-4 l ninO'N 

t\)443004ifi4^NN 

X 

Al 

0 

—4 

X 

X 

AJ 


X 

AJ 

O' 

AJ 

X 


a 3 


0 — nru'i'C'Otfl'orficotn 

■J-incO^rgrO't'Oi^a) 

vO 

X 

4- 

—4 

X 

X 


AJ 

X 

0 

X 

X 

0 



ooo— 4Aj>*-A-oinoA- 

n-t't'0(Mcl!0'f'4(J l N 

-0 

X 

4- 

X 


O' 


4- 

0 

A- 

AJ 

X 

4 - 

O' 


0000000—4— 4C\jaj 

lA X A- T O AJ 4 'A A- 

it 

—4 

X 

X 

A* 

X 

0 

X 

4- 

X 

A- 

X 

c 

1—4 





• 












• 

— O 
00 UJ 

000000000000 
1 1 1 1 1 1 1 1 1 1 1 

000000-4-4-4-4-4 

1 t 1 1 i 1 1 t 1 1 1 

Y 

AJ 

AJ 

1 

AJ 

1 

AJ 

1 

AJ 

X 

1 

X 

1 

X 

X 

1 

X 

1 

X 

1 

4" 

1 

4- 

1 

o. a 

> 

<a co *4 0 in rg in 0 * -1 co 

X 

'J’ 


X 

4- 

4- 

O' 

X 

X 

X 

X 

4- 

4- 

X 



x vj- x o — 1 -4 in -4 n. x x 

AJ 

X 

A- 

X 

4” 

0 

X 

AJ 

AJ 

X 

X 

0 

AJ 

AJ 


00 — 1 (M m CO m 0 N ON 

0 'mCD'tA 10 !ONinfC |0 

A- 

AJ 

X 

O' 


AJ 


0 - 

X 

X 

X 


O' 

X 

< 

000000-‘"4MfA-J- 

IA N- X O AJ -T LA f— CT -4 X 

>4- 

X 

r- 

X 

0 

-4 

AJ 

X 

X 

4- 

X 

X 

X 

A- 

i- j> 

UJ UJ 

• ••«•••••• •• 

000000000000 

O O O —4 —4 —I —4 —4 —4 A! AJ 

• 

AJ 

AJ 

Al 

AJ 

X 

X 

X 

X 

X 

X 

X 

X 

X 

t 

X 

x a 
>•» 

1 1 1 1 1 1 1 1 1 I 1 1 
OOOOOOHlM'KO-t 

1 1 1 1 1 1 1 1 1 1 1 
AJ >T — <XAJCOAJinXO-4 

1 

<r 

1 

X 

1 

0 

1 

0 

1 

0 

1 

O' 

1 

4* 

X 

1 

4 1 

1 

X 

1 

X 

1 

4- 

1 

O' 

1 


ooooooooooo 

fMninsocno.no'Nin 

X 

X 

X 

X 

X 

AJ 

. AJ 


O' 

X 

AJ 

A- 

0 

X 


ooooooooooo 

0000-4— <—4AJAJX4- 

in 

X 

A- 

X 

O' 

0 

—4 

X 

AJ 

X 

4- 

4- 

X 

X 


ooooooooooo 

ooooooooooo 

0 

0 

0 

0 

0 

—4 

-4 

—4 

—4 

—4 

—4 

—4 

—4 

—4 

AJ 

■ •■••*••*••• 

• • ••••••••• 


• 


• 


• 

t 

• 

• 

• 

• 

• 

• 


- 3 

000000000000 

ooooooooooo 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

X uj 

a. a 

« 1 1 1 t 1 1 1 I 1 

1 1 1 1 1 1 1 1 1 i 1 

1 

1 

1 

1 

I 

1 

1 

1 

1 

1 

1 

1 

1 

1 

Q O 

oo-^moiftco^Hin® 

— <iAOXXOX-44-XA- 

in 

AJ 

A- 

O' 

O' 

X 

O' 

AJ 

X 

O' 

X 

O' 

0 

O' 

X UJ 

OOOOO-HAJvfN-OsT 

0 O >f- N IM |C| >}• N 0 >r CO 

X 

X 

AJ 

X 

0 

4- 

A- 


4* 

A- 


X 

f— 4 

X 

a: i/> 

000000000-4-> 

A1AJX-J-IAXA-COO-4AJ 


X 

r- 

X 

0 

-4 

AJ 

4- 

X 

X 

X 

O' 


AJ 

N- A. 
aj O 

OOOOOOOOOOO 

OOOOOOOO- 4 - 4-4 

• ■••••••••• 

— « 

• 

~* 9 

• 

H 

• 

AJ 

• 

AJ 

• 

Al 

AJ 

Al 

AJ 

AJ 

AJ 

X 

X 

X u 

000000000000 

ooooooooooo 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

>— Q 

oo^ma'-iin-t'OOiN 

0 -O •!■ in H !MN -1 -4 CO 


X 

X 

0 

X 

AJ 

X 

0 

X 

—4 

X 

X 

O' 

X 

AJ 

OOOOO'M'tD't^J) 

msj , f44i4'n-iM^4in 

O' 

O' 

4- 

X 

AJ 

—4 

X 

4- 


O' 

0 

X 

O' 

AJ 

< 

OOOOOOOO— 4AJX 

inNO«lN4 0'4'fl^N 

Al 

A- 

AJ 

X 

O' 

—4 



O' 

X 

X 

X 

AJ 

A- 

1 - 3 

OOOOOOOOOOO 

00-4~4-4C\|AjXX4'4' 

X 

X 

X 

X 

X 

A- 

A- 

A- 

X 

X 

X 

X 

X 

4- 



















r o 

V- 

000000000000 
1 1 1 f 1 f 1 f 1 f 

OOOOOOOOOOO 

I f ( f f f r 1 1 1 1 

0 

1 

0 

1 

0 

1 

0 

( 

0 

1 

0 

1 

0 

1 

0 

1 

0 

1 

0 

1 

0 

1 

0 

1 

0 

1 

0 

1 


OO^AJ^N-X— 4XCTCT 

-0 X A- v* x — 1 AJ — 4 X AJ nJ- 


>4- 

0 

-4 

X 

4- 

0 


—4 

X 

X 

4* 

0 

X 


000000-4 (M .X 4* 0 

CNAJ.O-4XXOX XXX 

A- 

X 

0 

AJ 

4- 

A- 

0 

4" 

X 

AJ 

A- 

X 

0 

A- 


ooooooooooo 

0-4-4AJAJX-4-s}-inXA. 

X 

O' 

—4 

AJ 

X 

4- 

•0 

A- 

X 

0 

—4 

X 

X 

■0 

-u 

ooooooooooo 

OOOOOOOOOOO 

0 

0 

-4 

-4 


•4 

—4 

—4 

—4 

AJ 

AJ 

Al 

AJ 

AJ 






• 

• 

• 

• 

• 

• 

• 






X >JJ 

ooooooooooooooooooooooo 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

O 

a. o 




X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

o 

X X Aj AJ X -Q in S\ O' O' IA 

A- X ■04-OA-4'-4XX4' 

CT 

—4 

X 

AJ 

X 

CO 

A. 

X 

X 

A- 

—4 

O' 

0 



CO X CO CO X O X T O -4 X X 

(J ,| C|DsfH'ONNNO'nj 

in 

0 

X 

AJ 

0 

O' 

0 

AJ 

X 

T 

4* 

X 

4" 

O' 

x oo 

H4444444(MN'NJ1\J 

AixX4*ir\in-oA*xo'-4 

AJ 

•r 

n 


O' 

0 

X 

X 

A* 

O' 

AJ 

4- 


O' 

oc s 

O'Q'Q'O'&'J'O'Q'O'Q'O'O' 

J'O'O'CACTO'TiTCT'CT'O 

0 

0 

0 

0 

0 

-4 

—4 

-4 

—4 

—4 

AJ 

Al 

AJ 

AJ 

-« 3 
oO UJ 

• ••••••••••> 

AJAJAJAJnjr\JC\JAJAJAJAlAI 

"gAJAJAlAJAJAJAJAlAJ.X 

X 

■0 

X 

X 

X 

• 

X 

X 

X 

• 

X 

X 

X 

X 

X 

X 

a. o 

AlAJAJC'dAJvNJ-'JAJAlAJAJ A1 

AJ A| AJ AJ AJ AJ A1 AJ AJ AJ Al 

AJ 

AJ 

AJ 

Al 

AJ 

Al 

AJ 

AJ 

AJ 

Al 

AJ 

Al 

Al 

AJ 




AJ 

Al 

AJ 

AJ 

AJ 

AJ 

Al 

AJ 

Al 

AJ 

AJ 

AJ 

AJ 

AJ 


X X X X X X A| O X AJ O' O A- AJ 
O^OiAO'^ OSNHr« ^0-)-IAQS rf lC07'fO^ — < X O' O (D -ct- AJ AJ X aj .*■ X Al O' 

O O 'OtOHiACONinJ)-^^ vf 4- X X ~t O O' X A- 00 O' — • X -O' 

OOOOOOHHN^-t^'ONO'O-'rfl't'flN^LUTilOUJHijjiMmffl'iJ^uj^aJUlJJ^aJNiJj'O JJ O 'U -4 UU AJ JU 

• • • • • • • • • • • • • • • • • i * • • • c£ « i t a; • o' • g£ • ^ • CC • oC * cC * aj *2 £ • cC • i * cC 

000000000000000-«-*~«-*-4-4-H<r*^<AJ<:'g<AJ<Aj<cg<r\j<Al<Aj<Aj<AJ'3X<X<X< 


3000000 0 0000000 DOOOC 
300000000000000 3 0 0 0' 
"loxoxoxonoxoxoxoxox;; 
4iM^f^(n<t4 , inin^ONNcoBoa'oo< 

• •• • e » » » ••• • • ••••!• 

>0000000000000000-4-4- 


> _J 3 -J O . 
>300 00' 
><X<0- 
' O -T J X I 


10-10-10. 
i o x o x o 3 
io<*n<o« 
i-oo 0 o a- c 

: h ^ ^ ^ h « 

> o o 
; < < 


94 



TIME GAMMA ALPHA IAV82 ' THfcTABO PSIBO PHIBO XCG YCG ZCG XBR2CG Y8R2CG ZBR2CG 

SBC DEG DEG UEG/SEC OEG/SEC DEG/SEC OEG/SEC INCHES INCHES INCHES INCHES INCHES INCHES 


X 

a 

u. 


ooooooooo 

0 0 

0 

0000000000 

0 

0 

0 

0 

0 0 

0 

0 

0 

0 

0 

0 

0 

0 

O 

ooooooooo 

0000000000000 

0 

0 

0 

0 

0 0 

0 

0 

0 

0 

0 

0 

0 

0 

O 

00666660600 

0 

0000060000000000666006000000000060006660 

ooooooooo 

c 0 

0 

0000000000 

0 

0 

0 

0 

0 0 

0 

0 

0 

0 

0 

0 

0 

0 

O 

rgegegcMfMCMCMCMiM 

CM cm 

CM 

CMCyCgCMCMCMCMCMfMCM 

rvi 

04 

CM 

CM 

cm rg 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 


OOOOOOOOOOOOOOOOOOOOOO 

O 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

O 

0 

O O 

60OOOOOOOOOOOOOOOOOOOO 

O 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

O 

0 

O O 

666600660666066660006660606066666666 

.0 0 

606666 

ooooooco 

OOOOOOOOOOOOOOOOOOOOOO 

O 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

O O 

OOOOOOOOOOOOOOOOOOOOOO 

O 

0 

0 

0 

0 

• 0 

0 

0 

■ 0 

0 

0 

0 

0 

O O 

CMfMrMfMCMCMf\JfMCMCMfM<NJCMf\)(MCMCMCMCgCMCMCM 

eg 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM CM 

•O'O'O'O'O'O'O'O'O'O O'O'O'O'C'O'O'O'O'O'O'O 

O 

4>, 

>0 

O 

O 

43 

45 

45 

4) 

4} 

4) 

43 

45 

43 >0 

ininintninthinmminmintnininininmminmm 

in 

in 

in 

m 

in 

in 

IT. 

in 

m 

in 

in 

in 

in 

in m 

— r— r—j — i — f— 

0- 

7 . 

0. 

7 , 

0. 

7 , 

0 

7 . 

0 f- 

6n6p*6non6^6^on 

0 n- 

o e- o n- 

•— H f—4 rH — H —4 r—4 r—i r-4 H —4 r-4 r-H r-4 r-4 ■— 4 —4 —4 r-4 «-4 r-4 --4 

r-H 


r-H 







r-H 

—4 

•—4 

— 1 


N-e-n-N-f'-N-e-N-r'-N-N-N-N'n-N-r'-f'-.e-N-N-N-N- 


r- 

r- 

r- 

r- 

n- 

N- 

r- 

r- 

r- 

r» 

r- 

e- 

r» c- 

0000 0 0 0- HiNN(C|in> 07 -<'t(DNN 

ro 

O' 

43 

CO 


0 

O' 

O' 

0 

-—4 

cn 

in 

CO 

CM 43 

OOOOOOOOOOOOOOO— <^- 4 f\lrM 

ro 

0 

*4- 

m 

4> 

n* 

n- 

00 

0 

x4 

CM 

m 

4* 

4 > r- 

00066606660666006006060060 

0 0 

00000066666-1 

0-4 0 

H 

6 -« 

6 -4 

0—4 0—4 


4 -r^e- 4 - 43 rt 34-0 

OC-CM 45 (MrnfMrnomOCMM- 

-4 

-4 

co 

CO 

43 

4 * 

CM 

O' 

0* 

co 

N- 

CO 

CM 

43 

^4 

04 J-t 0 ^(DN® 

-4 n g- >0 •- i^-rM-Hr-cc'xjsO'n 

41 


4 ) 

43 

c*. 

4 - 

m 

0 

4 - 

eg 

CM 

CM 

—4 

m 

m 

ina '^43 43 incMin 

sOfMinfnr-inconinc 7 ' 4540 co 

CM 

CO 

4 - 


00 

m 


43 

O' 

—4 


O' 

in 

CO 

O' 

m-oocn 430 'CM 4 - 

vOOO'OOOO'cC'OrtO^H 

43 

O' 

CM 

4 - 

4 - 

4 - 

to 

O 


eg 

4 ) 

co 

q 

0 

O' 

OOO— 4-4— «— irsirg 

rgfMcgrncnfnfMrMfMCMCM- 4 -^ 

6 o 6 a* 

o o o ® 

ONO’ObinO't 

O CM 

0 -4 

6 6 

ON O 4 6 40 

r -4 

— <<MO'#-in 4 !r~ao 

0 ' 0 - 4 CMr<*i 4 'insOC-cOO' 0 — • 

CM 

CM 

CO 

4 - 

in 

4 ) 

r*» 

CD 

O' 

0 

0 



ro 

4 - 

1 1 1 l l 1 t 1 

1 -(-rhhh-h-hhhhnm 

1 1 1 1 1 1 1 1 1 1 i i 

CM 

1 

CM 

1 

eg 

1 

CM 

1 

eg 

1 

eg 

l 

CM 

CM 

1 

eg 

1 

to 

1 

CO 

1 

CO 

1 

ro 

I 

CO 

CO 

1 

4 -oOfnomco 4 -r^ 

egmCMCgrM-<N- 0 ' 43 r>-CMC\j 43 

e- 

O' 


CO 

co 

4 - 

O 

O' 

r- 

4 - 

4 " 

O' 

O' 

N- 

O 

43 m rO 43 m r— OO 

mr\i<piHr-(M 4 -'coNON'J' 

CO 

4 - 

e- 

4 - 

N- 

m 



4 ) 

4 - 

4 - 


N» 

CO 

m 

inm'tNONmtj'-t 

( 7 « 4 ,r Df'JinO'(VJ 4 )CMCICONN 

<o 

o 

r- 

in 

m 

CM 

CM 

eg 

4 - 

CO 

4 - 

ro 

4 - 

O' 

CO 

in 4 -* 4 i O' 0 co ^ O' 

- 4 CMOr^- 4 m«tcMoomin 4 )- 4 ‘ 

—4 

43 

CO 

O' 

CO 

in 

O 

to 

4 - 

to 

-4 

r- 

—4 

<0 

4 - 

r>.NNsO(figfAfvj 6 

CM^iniN 0 N'f- 4 N 4 ' 06 cj 6 ffl 6 CO ‘6 aaOcnoOOOCdOCOOCMOOOOO^ 

On o-*o 4 0 

N- 


ooooooa*CT*'»aoo 3 /».r- 

43 

4 > 

in 

in 

4 - 

4 - 

to 

CO 

eg 

eg 


0 

O 

O' 

CO 

minminminmminminininm 4 - 4 - 4 - 4 - 4 ' 4 - 4 - 4 * 
1 1 , 1 1 1 1 t t 1 9 1 1 1 1 1 1 1 > • 1 1 1 

4 * 

1 

M- 

1 

4 - 

1 

r 

4 - 

1 

4 - 

1 

4 - 

1 

4- 

I 

4 " 

1 

4 “ 

1 

4 - 

T 

4 - 

1 

<0 

CO 

1 

oo-HfnKio'iflco 

^Nooo'ON'tcomo' men 

o 

CM 


00 

CM 

CM 

43 

in 

co 

in 

4 ) 

CO 

m 

4 ) 

4 - 

000000 - 4-4 

Nng-^insooO'O'MnKis 

O' 

o 

-4 


CM 

eg 


O 

CO 

43 

ro 

O 

43 

eg 

-4 

0000000, 0 

00000000 - 4 - 4 -^,- 4-4 

_ 

CM 

CM 

CM 

CM 

CM 

eg 

eg 



-« 

-4 

0 

O 

O 

600006000000006000000000006066000 
1 1 , 1 1 1 1 1 1 1 1 1 f 1 1 1 9 t 1 t 1 1 1 1 1 1 1 1 

oooooooooooooooo 
1 1 I I 1 I 1 1 

000 

1 

i ! 

oo®fnom 4 (MO' 

0 'Om>ONin,t«>a 3 rMNN,o 

co 

tD 

r- 

43 

e- 

O' 

co 

43 

CM 

CO 

CO 

m 

CO 

4 -, 

-4 

Orci-H 4343 Oa 50 'O 

HO'-‘(p|M>} , g , M)nin-«>o 

4 - 

o 

CM 

co 

43 

4 > 

O 

m 

O' 

O' 

CM 

m 

-4 

4 “ 

4 * 

O-Hin— iOcnjit, -<o 

o -ma O' ^ ij >f a -< o- N 

CM 

r- 

CO 



4 - 

43 

4 ) 

43 

N- 

O 

4 - 

—4 

O 

CM 

OOOHN|f|^vO» 

0 <MM- 450 '- 4 M'> 00 '- 4 rnM' 4 ) 

N- 


e- 

r» 

r- ■ 

4 ) 

in 

4 - 

tO 

CM 

CM 

—4 

-4 

”* 

—4 

606066666 
I I 1 1 1 I I 1 

.H^rH-H-ipj^pgrjfnnfnnofOofl 
1 .1 1 t l .1 t 1 1 1 1 1 1 1 1 

6 ro o 6 
1 1 

omomomoroomomomorn 
l i I l l i I I 

0 m 0 co 0 co 
l 1 1 


e- 

4 * 

4 ’ 

O' 

m 

m 

in 

~4 

00 

r 4 

•O 

—4 

0 

O* . 

4 * 

4‘aoootococoN-O' 

inco[MOMo- 40 «- 4 'trsB 

m 

O' 

o 

4 - 

CM 

to 

0 

in 

CO 

—4 

O' 

co 

O' 

O 

O 

CMO-HCOcocMOOCeO-C— lO' 4 ‘ 4 Jf 0 Nf-«-«in<O 

4 ) 

4 3 

in 

—4 

r- 

cn 

-4 

r— < 

m 

m 

O' 

O 

N- 


O 

ooegmmcooro 

•ooifl'OaHcnifiONNoin 

m 


O' 


4 " 

eg 

0 

CO 

43 

m 

4 * 

in 

m 

N- 

CO 

6606666-4-4 
II II 1 1 II 

- 4 CMCMCMcgfPim 66 mmm 66 fnofn 6 cM 

t 1 1 1 1 I 1 1 1 1 1 , 1 1 1 1 i 1 

6 eg 6 eg 6 CM 
1 I I 

6 CM O 
1 

1 

0 - 40-4 
1 1 

OHOH 
1 1 

6-46 

l 

-4 O 
1 

-4 

1 

COco6 COCOQDO'O'Oeg.J- 4 ) 0 >t' 0 ' ^ ^ (j 1 N Nco O 

<o 

o 

—4 

43 

m 1 

-4 


eg 

4 1 

e- 

-4 

43 

-4 

43 

eg 


cMcgcMrnmfnsfmtn. 4 >r *-330 

—4 

CM 

4 - 

in 

r- 

O' 

—4 

to 

in 

e- 

O 

CM 

m 

N- 

O 

O' O O' O' O' O' ^ O' 0* O* < 7* O' O'* O' O' o* o* o* 0 

o 

o 

o 

o 

o 

o 

— * 

— H 

-4 

-4 

CM 

eg 

CM 

CM 

ro 

cMrgcMcgcgcgcMrgcM 

CMCMCMCMCMnjpgfMogegrMCMrOOcOOtOOcOOcOOcOOrOOcOO 

cn 

6 6 

0 co 

0606 

O ro 6 CO 6 CO 

CMCMCMCMCMNCMCMCMCMCMCgcMCMCMCMCMCMCMCMcgiM 

CM 

eg 

CM 

<M 

CM 

1 

eg 

CM 

eg 

eg 

CM 

CM 

eg 

eg 

eg 

CM 

O 0 O O 0 a 0 O io O 0 0 ~4 i-H -4 -4 CM CM CM 'CO <0 

co 

* 

4 - 

4 - 

in 

m 

in 

>0 

O 

43 

r- 

f- 

N- 


co 

660006666006660606660666666666 

o 0.0 

6 6 6 6 6 6 06666 

06006666 


O-4<0'»Oc0Oina>r0Or0r0O4’in-H-4egrMr^ 

4- 

-4 

m 

CM 

eg 

4* 

CO 

43 

O' 

-4 

43 

» 

CO 

in 

O' 

0000410 ^ 4 N 4 IN 04 inOON ■ 0 NO 4 ’ 

0 

N- 

fO 

O' 

4 1 

D 

-4 

4* 

0 

O' 

—4 

4- 

O* 

in 

fO 

OOOOOOO-^-*egr0m4333Or0mJ3-4m<D 

eg 

m 

O' 

CM 

43 

O* 

(O 

43 

O' 

eg 

43 

O' 

CM 

43 . 

O 

00000000000000 - 4 - 4 -H— irgcMeg 

10 


to 

4^ 

4- 

4* 

in 

m 

in 

43 

43 

43 

N» 

N- 

'tt 

66666666006660660606666660606000060006000600 

66606666 

OOOOOOOOO OOOOOOOOOOOO 

0 

O 

0 

O 

O 

O 

0 

O: 

0 

O 

O 

O 

O 

O 

O 

momoinoinomo^nomoinomomom 

0 

in 

0 

in 

O 

m 

0 

m 

0 

m 

O 

in 

O 

in 

O 

O-4i4cg(\}ff|cCi4'4’inino4)N*N , c0»0'0'OO 

-4 
• * 

•4 
a • 

eg 

• • 

CM 

• • 

CO 

• • 

CO 

t a 

4- 

• * 

4- 

a a 

m 

a • 

m 

• • 

43 

43 

.at 

N- 

N- 

• • ! 

CO 
• * 


00000000.0 00000000000 -<- 40 - 40 ^ 0 - 40 - 40 *- 40 — * 0 - 40 — < 0 — < 0 — ‘O— 4 0 -* 0 - 4 Q-« 0 -« 


I 


95 


Figure B-9. - Continued. 
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Figure B-9. - Concluded 






Figure B-10. - Pertinent S-C 4020 graphical output for sample run four. 
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Figure B-10. - Continued. 
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Figure B- 10. - Continued. 
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Figure B- 10. - Continued. 
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Figure B-10.- Continued. 
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Figure B-10. - Continued. 
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Figure B-10. - Continued. 
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Figure B- 10. - Continued. 
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Figure B- 10. - Concluded. 
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CABLE CONNECTED SPACE STATION DYNAMICS 

WEIGHT OF BODY 1 = 39285 LB. 

WEIGHT OF BODY 2 = 17935 LB. 

INITIAL SPIN SPEED = 22.92 DEG. /SEC. 

INITIAL DISTANCE BETWEEN BODY C.G.*S= 137. G FT. 
CABLE ELASTICITY = 14583300 PSI 

TOTAL CABLE AREA = 0.1512 SQ. IN. 

Figure B-ll. - Typical S-C 4020 motion picture output for sample run four. 
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TORQUE ABOUT Ij = 333 SINCQ«>25 TD FT. -LB 

TORQUE ABOUT ^ = 125000 SINC0.25 T) FT. -LB 

TORQUE ABOUT = 125000 SINC0.25 TD FT. -LB 

Figure B-ll. - Continued. 
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